Skip to content
Snippets Groups Projects
Commit 2718dcf6 authored by cmaffeo2's avatar cmaffeo2
Browse files

Fixed bug in product potential's apply_dihedral_force

parent 6f782b6d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/TabulatedPotential.h
+ 5
− 5
View file @ 2718dcf6
...@@ -440,11 +440,11 @@ public: ...@@ -440,11 +440,11 @@ public:
f3 = -distbc * crossBCD.rLength2() * crossBCD; f3 = -distbc * crossBCD.rLength2() * crossBCD;
f2 = -(ab.dot(bc) * bc.rLength2()) * f1 - (bc.dot(cd) * bc.rLength2()) * f3; f2 = -(ab.dot(bc) * bc.rLength2()) * f1 - (bc.dot(cd) * bc.rLength2()) * f3;
energy_deriv = (ab.length2()*bc.length2()*crossABC.rLength2() > 100.0f || bc.length2()*cd.length2()*crossBCD.rLength2() > 100.0f) ? 0.0f : energy_deriv; // energy_deriv = (ab.length2()*bc.length2()*crossABC.rLength2() > 100.0f || bc.length2()*cd.length2()*crossBCD.rLength2() > 100.0f) ? 0.0f : energy_deriv;
if ( energy_deriv > 1000.0f ) /* if ( energy_deriv > 1000.0f ) */
energy_deriv = 1000.0f; /* energy_deriv = 1000.0f; */
if ( energy_deriv < -1000.0f ) /* if ( energy_deriv < -1000.0f ) */
energy_deriv = -1000.0f; /* energy_deriv = -1000.0f; */
f1 *= energy_deriv; f1 *= energy_deriv;
f2 *= energy_deriv; f2 *= energy_deriv;
......
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