[gmx-users] Re: A bug for MPI of GROMACS?

Mark Abraham Mark.Abraham at anu.edu.au
Wed Aug 12 05:38:27 CEST 2009


wuxiao wrote:
> Dear Mark and Vitaly,
>   In order to get help, here I describe the more detailed procedure. 
>   I successfully generated the gro and top files using pdb2gmx.
>   pdb2gmx -f -o -p -ter
>   And also the grompp finished successfully.
>   grompp -f -c -p -o 
>   Then I ran the mdrun in parallel way, only to find that the polymer became many pieces.

I doubt that it did. The representation of it that you were viewing may
have made what was one topology look like many pieces.

>   mpirun -np 4 mpimdrun -s -o -c -e -g -v
>   Obviously, it cannot form any bonds between those pieces due to the long distances.
>   With the same tpr file,I ran the mdrun in non-parallel-way, however, the polymer seems to behave properly. 
>   mdrun -s -o -c -e -g -v
>   The partial  files (gro, top, and mdp) are given as follows:
> 
>   gro file:
> Great Red Owns Many ACres of Sand
>  1382
>     1DIA  N1  1  -0.199  0.069  0.078
>     ......
>     125AMI H42  1382 -0.073  -0.681  -1.613
>      3.61500  4.27230  2.77190

Since you are doing simulating a periodic system, are these box
dimensions reasonable?

Mark

>   top file:
>     ......
>     [ molecules ]
>     ; Compounds #moles
>     Protein_D 1
> 
>   mdp file:
>     ......
> define = -DFLEXIBLE
> constraint = none
> integrator = steep
> nsteps = 10000
> ns_type = grid
> nstcgsteep = 100
> nslist = 10
> rlist = 1.0
> coulombtype = PME
> rcoulomb = 1.0
> vdwtype = cut-off
> rvdw = 1.0
> DispCorr = EnerPres
> fourierspacing = 0.36
> pme_order = 4
> ewald_rtol = 1e-5
> emtol = 10
> emstep = 0.01
> 
>  
> 
> _________________________________________________________________
> 上Windows Live 中国首页,下载最新版Messenger!
> http://www.windowslive.cn
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php



More information about the gromacs.org_gmx-users mailing list