Clamp exponent to avoid overflow warnings. [ref #8]

Author: lohedges

src/loch/_sampler.py

@@ -1401,7 +1401,11 @@ def move(self, context: _openmm.Context) -> list[int]:
 
                         # Compute the PME acceptance correction.
                         acc_prob = _np.exp(
-                            -self._beta_openmm * (final_energy - initial_energy - dE_RF)
+                            min(
+                                0.0,
+                                -self._beta_openmm
+                                * (final_energy - initial_energy - dE_RF),
+                            )
                         )
 
                         # Store the PME energy change and acceptance probability.
@@ -1482,7 +1486,11 @@ def move(self, context: _openmm.Context) -> list[int]:
 
                         # Compute the PME acceptance correction.
                         acc_prob = _np.exp(
-                            -self._beta_openmm * (final_energy - initial_energy - dE_RF)
+                            min(
+                                0.0,
+                                -self._beta_openmm
+                                * (final_energy - initial_energy - dE_RF),
+                            )
                         )
 
                         # Store the PME energy change and acceptance probability.