I was wondering whether I am implementing omm incorrectly in my dynamics script, because I am not noticing a speed-up when compared to my CPU run. I'm using the CHARMM22_cmap forcefield.
...
scalar fbeta set 10 select .not. hydrogen end
prnlev 0
stream "/toppar/gbsw/radius_gbsw.str"
prnlev 5 node 0
gbsw sgamma 0.005 nang 50 conc 5e-2
set ctofnb = 12
calc ctonnb = @ctofnb
calc cutnb = @ctofnb + 2
NBOND atom switch cdie vdw vswitch -
ctonnb @ctonnb ctofnb @ctofnb cutnb @cutnb -
mini sd nstep 100
energy
write psf card name prot.psf
write coor pdb name prot_minimized.pdb
omm on
faster on
shake fast bonh
open unit 1 write form name prot_omm.res
open unit 2 write unform name prot_omm.dcd
dynamics leap start timestep 0.002 -
nstep 25000000 nprint 10000 iprfrq 20000 -
firstt 298 finalt 298 twindl -5.0 twindh 5.0 -
ichecw 0 teminc 0 ihtfrq 0 ieqfrq 0 -
iasors 1 iasvel 1 iscvel 0 -
ilbfrq 0 inbfrq -1 echeck -1 -
iunread 1 iunwrite 1 -
iuncrd 2 nsavc 100 nsavv 0 -
omm langevin gamma 5
open unit 1 write form name prot_omm.res
open unit 2 write unform name prot_omm.dcd
dynamics leap restart timestep 0.002 -
nstep 225000000 nprint 10000 iprfrq 20000 -
firstt 298 finalt 298 twindl -5.0 twindh 5.0 -
ichecw 0 teminc 0 ihtfrq 0 ieqfrq 0 -
iasors 1 iasvel 1 iscvel 0 -
ilbfrq 0 inbfrq -1 echeck -1 -
iunread 1 iunwrite 1 -
nsavc 100 nsavv 0 iuncrd 2 -
omm langevin gamma 5
close unit 1
close unit 2
stop
