Commit 6e306f2c authored by peastman's avatar peastman
Browse files

Merge branch 'master' into cleanup

parents eec14ddb ba9f33c7
...@@ -396,7 +396,6 @@ void testVdwAmmoniaCubicMeanHhg() { ...@@ -396,7 +396,6 @@ void testVdwAmmoniaCubicMeanHhg() {
compareForcesEnergy(testName, expectedEnergy, energy, expectedForces, forces, tolerance); compareForcesEnergy(testName, expectedEnergy, energy, expectedForces, forces, tolerance);
} }
// test VDW w/ sigmaRule=CubicMean and epsilonRule=W-H // test VDW w/ sigmaRule=CubicMean and epsilonRule=W-H
void testVdwAmmoniaCubicMeanWH() { void testVdwAmmoniaCubicMeanWH() {
...@@ -435,6 +434,12 @@ void testVdwAlchemical(int power, double alpha, double lambda, AmoebaVdwForce::A ...@@ -435,6 +434,12 @@ void testVdwAlchemical(int power, double alpha, double lambda, AmoebaVdwForce::A
std::string testName = "testVdwAlchemical"; std::string testName = "testVdwAlchemical";
int numberOfParticles = 8;
double boxDimension = -1.0;
double cutoff = 9000000.0;
std::vector<Vec3> forces;
double energy;
setupAndGetForcesEnergyVdwAmmonia2("CUBIC-MEAN", "HHG", cutoff, boxDimension, forces, energy, setupAndGetForcesEnergyVdwAmmonia2("CUBIC-MEAN", "HHG", cutoff, boxDimension, forces, energy,
method, power, alpha, lambda); method, power, alpha, lambda);
std::vector<Vec3> expectedForces(numberOfParticles); std::vector<Vec3> expectedForces(numberOfParticles);
...@@ -459,7 +464,6 @@ void testVdwAlchemical(int power, double alpha, double lambda, AmoebaVdwForce::A ...@@ -459,7 +464,6 @@ void testVdwAlchemical(int power, double alpha, double lambda, AmoebaVdwForce::A
compareForcesEnergy(testName, expectedEnergy, energy, expectedForces, forces, tolerance); compareForcesEnergy(testName, expectedEnergy, energy, expectedForces, forces, tolerance);
} }
// test VDW w/ sigmaRule=Arithmetic and epsilonRule=Arithmetic // test VDW w/ sigmaRule=Arithmetic and epsilonRule=Arithmetic
void testVdwAmmoniaArithmeticArithmetic() { void testVdwAmmoniaArithmeticArithmetic() {
......
...@@ -326,7 +326,7 @@ void testArgonBox() { ...@@ -326,7 +326,7 @@ void testArgonBox() {
integrator.setMaximumStepSize(maxStep); integrator.setMaximumStepSize(maxStep);
for (int i = 0; i < 100; i++) { for (int i = 0; i < 100; i++) {
integrator.step(1); integrator.step(1);
ASSERT(integrator.getStepSize() <= maxStep); ASSERT(integrator.getStepSize() <= maxStep*1.000001);
} }
} }
......
...@@ -305,7 +305,7 @@ void testArgonBox() { ...@@ -305,7 +305,7 @@ void testArgonBox() {
integrator.setMaximumStepSize(maxStep); integrator.setMaximumStepSize(maxStep);
for (int i = 0; i < 100; i++) { for (int i = 0; i < 100; i++) {
integrator.step(1); integrator.step(1);
ASSERT(integrator.getStepSize() <= maxStep); ASSERT(integrator.getStepSize() <= maxStep*1.000001);
} }
} }
......
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