be consistend and use one sided difference if outside

This commit is contained in:
Ahmet Inan 2020-01-17 19:02:15 +01:00
commit e4b687d1f6

View file

@ -102,7 +102,7 @@ struct CubicHermiteSpline
: k < n ?
eval(X, Y, (x - X[k-1]) / (X[k] - X[k-1]), k-1, n)
:
eval(X, Y, (x - X[n-2]) / (X[n-1] - X[n-2]), n-2, n);
eval(X, Y, (x - X[n-2]) / (X[n-1] - X[n-2]), n-1, n);
}
static constexpr OTYPE left(const OTYPE *Y)
{