[gmx-users] MPI problem------Cannot compile and link MPI code with cc

Gavin Melaugh gmelaugh01 at qub.ac.uk
Mon Mar 29 16:52:07 CEST 2010


Hi guys

I have been running gromacs for a few months now with no real problems
however some of my numbers differ from my supervisor's who uses DL-POLY.
Anyway discussion of the discrpnacies has led me to questioning aspects
of my topology file below (exerpts). We created the topology ourselves.
I always thought that the set up below would calculate 1-4 interactions,
but having read through the archives I am now confused. Do I have to
give a pair list in topology file, or is gen-pairs = yes sufficient. I
never created a pairlist as I thought that gen-pairs would create it
from the parameters in atomtypes.

;Parameter level
[defaults]
; nbfunc    comb-rule     gen-pairs        fudgeLJ     fudgeQQ
     1             3                    yes                   
0.5           0.5

[atomtypes]
;type     mass           charge      ptype     sigma(nm)    
epsilon(kjmol-1)
   CB     12.011000      0.000000       A      0.355000      0.292880
   CA     12.011000     -0.115000       A      0.355000      0.292880
   HC      1.008000      0.115000       A      0.242000      0.125520
   CU     13.019000      0.265000       A      0.350000      0.334720
   NU     14.007000     -0.597000       A      0.325000      0.711280
   CH     13.019000      0.332000       A      0.385000      0.334720
   C3     15.035000      0.000000       A      0.390500      0.732200
   C2     14.027000      0.000000       A      0.390500      0.493712

;Molecular level
[moleculetype]
;       name         nrexcl
       subhexylcage      3

[atoms]
; atomnr  type   resnr  residue   name    cgnr     charge     mass
    1      CA       1    CGE       CA       1     -0.1150     12.0110
    2      CB       1    CGE       CB       1      0.0000     12.0110

[bonds]
;  ai    aj   func   b0(nm)       kb(kjmol-1 nm-2)
    1     2     1     0.1400      392459.2
    1     6     1     0.1400      392459.2

[angles]
;  ai    aj    ak   func  theta(deg)    ktheta(kjmol-1 rad-2)
    2     1     6     1     120.00      711.2800
    2     1     7     1     120.00      292.8800

[dihedrals]
;  ai    aj    ak    al    func       V1           V2            
V3            V4 (all kj mol-1)
    1     2     3     4       5       0.00000      30.33400      
0.00000       0.00000
    1     2     3    10       5       0.00000      30.33400      
0.00000       0.00000

;System level
[system]
;the name of this system
config test of 32 subhexyl cage molecules

[molecules]
; moleculename     number
 subhexylcage        32




More information about the gromacs.org_gmx-users mailing list