OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

DenseAD: add copyDerivatives(x) and clearDerivatives() methods to Evaluation

Browse Source 
this makes it possible to implement these functions more
efficiently...
This commit is contained in:
Andreas Lauser
2016-10-27 21:45:40 +02:00
parent a289230202
commit 48d570afc5
2 changed files with 37 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
tests/test_densead.cpp
View File

@@ -66,6 +66,34 @@ void testOperators(const Scalar tolerance)
const Eval xEval = Eval::createVariable(x, 0);
const Eval yEval = Eval::createVariable(y, 1);
Eval xyEval = xEval;
Eval yxEval = yEval;
xyEval.copyDerivatives(yxEval);
yxEval.clearDerivatives();
for (int i = 0; i < numVars; ++i) {
if (i == 0 && xEval.derivative(i) != 1.0)
throw std::logic_error("oops: createVariable");
else if (i == 1 && yEval.derivative(i) != 1.0)
throw std::logic_error("oops: createVariable");
else if (i == 1 && xyEval.derivative(i) != 1.0)
throw std::logic_error("oops: copyDerivatives");
else continue;
// the remaining derivatives must be zero
if (xEval.derivative(i) != 0.0)
throw std::logic_error("oops: createVariable");
if (yEval.derivative(i) != 0.0)
throw std::logic_error("oops: createVariable");
if (cEval.derivative(i) != 0.0)
throw std::logic_error("oops: createConstant");
if (xyEval.derivative(i) != 0.0)
throw std::logic_error("oops: copyDerivatives");
if (xyEval.derivative(i) != 0.0)
throw std::logic_error("oops: clearDerivatives");
}
// test the non-inplace operators
{
Eval a = xEval + yEval;