Previous Thread
Next Thread
Print Thread
Joined: Jun 2016
Posts: 24
E
Forum Member
OP Offline
Forum Member
E
Joined: Jun 2016
Posts: 24
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

Joined: Sep 2003
Posts: 8,606
Likes: 24
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,606
Likes: 24
The output log file should have clear messages about the OpenMM and GPU usage.


Rick Venable
computational chemist

Joined: Sep 2003
Posts: 4,850
Likes: 10
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,850
Likes: 10
gbsw is not supported with OpenMM


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden

Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.009s Queries: 20 (0.006s) Memory: 0.7375 MB (Peak: 0.7830 MB) Data Comp: Off Server Time: 2022-10-06 17:12:34 UTC
Valid HTML 5 and Valid CSS