OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add test cases for multiplication and division

Browse Source 
Helped uncover syntax error (misprint) fixed in commit 8105166.
This commit is contained in:
Bård Skaflestad 2013-04-29 14:52:38 +02:00
parent 8105166c26
commit a4cdbdc049
1 changed files with 74 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

74
test_syntax.cpp
View File

@ -97,3 +97,77 @@ BOOST_AUTO_TEST_CASE(Subtraction)
BOOST_CHECK_CLOSE(omb.val(), 1 - b.val(), atol);
BOOST_CHECK_CLOSE(omb.der(), - b.der(), atol);
}
BOOST_AUTO_TEST_CASE(Multiplication)
{
typedef AutoDiff::Forward<double> AdFW;
const double atol = 1.0e-14;
AdFW a(1.0), b(1.0);
AdFW no_a = a * 0;
BOOST_CHECK_CLOSE(no_a.val(), 0.0, atol);
BOOST_CHECK_CLOSE(no_a.der(), 0.0, atol);
AdFW atb = a * b;
BOOST_CHECK_CLOSE(atb.val(), a.val() * b.val(), atol);
BOOST_CHECK_CLOSE(atb.der(), a.der()*b.val() + a.val()*b.der(), atol);
double av = a.val();
double ad = a.der();
a *= b;
BOOST_CHECK_CLOSE(a.val(), av * b.val(), atol);
BOOST_CHECK_CLOSE(a.der(), ad*b.val() + av*b.der(), atol);
av = a.val();
ad = a.der();
a *= 1;
BOOST_CHECK_CLOSE(a.val(), av, atol);
BOOST_CHECK_CLOSE(a.der(), ad, atol);
AdFW bto = b * 1; // b times one
BOOST_CHECK_CLOSE(bto.val(), b.val(), atol);
BOOST_CHECK_CLOSE(bto.der(), b.der(), atol);
AdFW otb = 1 * b; // one times b
BOOST_CHECK_CLOSE(otb.val(), b.val(), atol);
BOOST_CHECK_CLOSE(otb.der(), b.der(), atol);
}
BOOST_AUTO_TEST_CASE(Division)
{
typedef AutoDiff::Forward<double> AdFW;
const double atol = 1.0e-14;
AdFW a(10.0), b(1.0);
AdFW aob = a / b;
BOOST_CHECK_CLOSE(aob.val(), a.val() * b.val(), atol);
const double res = ((a.der()*b.val() - a.val()*b.der()) /
(b.val() * b.val()));
BOOST_CHECK_CLOSE(aob.der(), res, atol);
double av = a.val();
double ad = a.der();
a /= b;
BOOST_CHECK_CLOSE(a.val(), av * b.val(), atol);
BOOST_CHECK_CLOSE(a.der(), res, atol);
av = a.val();
ad = a.der();
a /= 2;
BOOST_CHECK_CLOSE(a.val(), av / 2, atol);
BOOST_CHECK_CLOSE(a.der(), ad / 2, atol);
AdFW bot = b / 2; // b over two
BOOST_CHECK_CLOSE(bot.val(), b.val() / 2, atol);
BOOST_CHECK_CLOSE(bot.der(), b.der() / 2, atol);
AdFW otb = 2 / b; // two over b
BOOST_CHECK_CLOSE(otb.val(), 2 / b.val(), atol);
BOOST_CHECK_CLOSE(otb.der(), -2*b.der() / (b.val() * b.val()), atol);
}