Skip to content
Snippets Groups Projects
Commit a25e4e17 authored by Samuel Scherrer's avatar Samuel Scherrer Committed by Timo Koch
Browse files

[common] fix segment index issue in spline extrapolation

parent 6d430f90
No related branches found
No related tags found
1 merge request!2235[common] fix segment index issue in spline extrapolation
......@@ -145,12 +145,12 @@ public:
// handle extrapolation
if (extrapolate) {
if (x < xMin()) {
Scalar m = evalDerivative(xMin(), /*segmentIdx=*/0);
Scalar m = evalDerivative_(xMin(), /*segmentIdx=*/0);
Scalar y0 = y_(0);
return y0 + m*(x - xMin());
}
else if (x > xMax()) {
Scalar m = evalDerivative(xMax(), /*segmentIdx=*/numSamples_()-1);
Scalar m = evalDerivative_(xMax(), /*segmentIdx=*/numSamples_()-2);
Scalar y0 = y_(numSamples_()-1);
return y0 + m*(x - xMax());
}
......@@ -179,7 +179,7 @@ public:
if (Dune::FloatCmp::le(x, xMin()))
return evalDerivative_(xMin(), 0);
else if (Dune::FloatCmp::ge(x, xMax()))
return evalDerivative_(xMax(), numSamples_() - 1);
return evalDerivative_(xMax(), numSamples_() - 2);
}
return evalDerivative_(x, segmentIdx_(x));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment