Anonymous user
Thiele's interpolation formula: Difference between revisions
Updated D entry
(Updated D entry) |
(Updated D entry) |
||
Line 357:
const real b, e, s;
auto range() const pure /*nothrow*/ @safe /*@nogc*/ {
return iota(b, e + s, s);
}
}
real eval0(alias RY, alias X, alias Y)(in real x) pure nothrow @safe @nogc {
real a = 0.0L;
foreach_reverse (immutable i; 2 .. X.length - 3)
Line 372:
immutable real[] Y, X, rhoY, rhoX;
this(real[] y, real[] x) immutable pure nothrow /*
in {
assert(x.length > 2, "at leat 3 values");
assert(x.length == y.length, "input arrays not of same size");
} body {
this.Y = y.idup;
this.X = x.idup;
rhoY = rhoN(Y, X);
Line 383:
}
this(in real function(real) pure nothrow @safe @nogc f,
Domain d = Domain(0.0L, 1.55L, 0.05L))
immutable pure /*nothrow @safe*/ {
auto xrng = d.range.array;
this(xrng.map!f.array, xrng);
}
auto rhoN(immutable
pure nothrow @safe {
immutable int N = x.length;
auto p = new real[][](N, N);
| |||