We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 0ecc721 commit 36a4bdfCopy full SHA for 36a4bdf
1 file changed
src/somd2/runner/_repex.py
@@ -1265,9 +1265,6 @@ def _run_block(
1265
),
1266
)
1267
1268
- # Set the state.
1269
- self._dynamics_cache.save_openmm_state(index)
1270
-
1271
# Perform a GCMC move and write ghost water residue indices after
1272
# dynamics so that the ghost state is temporally consistent with
1273
# the saved frame.
@@ -1289,6 +1286,9 @@ def _run_block(
1289
1286
# Save the GCMC state.
1290
1287
self._dynamics_cache.save_gcmc_state(index)
1291
1288
+ # Save the OpenMM state after any GCMC move so the context is consistent.
+ self._dynamics_cache.save_openmm_state(index)
+
1292
# Get the energy at each lambda value.
1293
energies = dynamics._current_energy_array()
1294
0 commit comments