diff --git a/src/TabulatedPotential.h b/src/TabulatedPotential.h index 0aa73a0939163c1cbeac7e83bc6e74dedb00f5bb..bebe922e288c6062fb66d74ee3c4bdee5f123ffc 100644 --- a/src/TabulatedPotential.h +++ b/src/TabulatedPotential.h @@ -440,11 +440,11 @@ public: f3 = -distbc * crossBCD.rLength2() * crossBCD; 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; - if ( energy_deriv > 1000.0f ) - energy_deriv = 1000.0f; - if ( energy_deriv < -1000.0f ) - energy_deriv = -1000.0f; + // 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 ) */ + /* energy_deriv = 1000.0f; */ + /* if ( energy_deriv < -1000.0f ) */ + /* energy_deriv = -1000.0f; */ f1 *= energy_deriv; f2 *= energy_deriv;