Unverified Commit 6bdaa495 authored by David Clark's avatar David Clark Committed by GitHub
Browse files

Increases hydrogen mass for benchmarks (#3208)



* Changes hydrogen mass for benchmarks

* Adds amber case for hydrogen mass
Co-authored-by: default avatarDavid Clark <daclark@nvidia.com>
parent ae2fe2fd
...@@ -88,6 +88,7 @@ def runOneTest(testName, options): ...@@ -88,6 +88,7 @@ def runOneTest(testName, options):
method = app.PME method = app.PME
cutoff = options.cutoff cutoff = options.cutoff
constraints = app.HBonds constraints = app.HBonds
hydrogenMass = 1.5*unit.amu
system = prmtop.createSystem(nonbondedMethod=method, nonbondedCutoff=cutoff, constraints=constraints) system = prmtop.createSystem(nonbondedMethod=method, nonbondedCutoff=cutoff, constraints=constraints)
if options.ensemble == 'NVE': if options.ensemble == 'NVE':
integ = mm.VerletIntegrator(dt) integ = mm.VerletIntegrator(dt)
...@@ -106,7 +107,6 @@ def runOneTest(testName, options): ...@@ -106,7 +107,6 @@ def runOneTest(testName, options):
else: else:
method = app.CutoffPeriodic method = app.CutoffPeriodic
cutoff = 1*unit.nanometers cutoff = 1*unit.nanometers
hydrogenMass = 1.5*unit.amu
elif explicit: elif explicit:
ff = app.ForceField('amber99sb.xml', 'tip3p.xml') ff = app.ForceField('amber99sb.xml', 'tip3p.xml')
pdb = app.PDBFile('5dfr_solv-cube_equil.pdb') pdb = app.PDBFile('5dfr_solv-cube_equil.pdb')
...@@ -132,6 +132,7 @@ def runOneTest(testName, options): ...@@ -132,6 +132,7 @@ def runOneTest(testName, options):
else: else:
dt = 0.004*unit.picoseconds dt = 0.004*unit.picoseconds
constraints = app.HBonds constraints = app.HBonds
hydrogenMass = 1.5*unit.amu
if options.ensemble == 'NVE': if options.ensemble == 'NVE':
integ = mm.VerletIntegrator(dt) integ = mm.VerletIntegrator(dt)
else: else:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment