30 void Foam::SIBS::polyExtrapolate
45 for (
register label j=0; j<n; j++)
47 dy[j] = yz[j] = yest[j];
52 for (
register label j=0; j<n; j++)
61 for (
register label k1=0; k1<iest; k1++)
63 scalar delta = 1.0/(x[iest - k1 - 1] - xest);
64 scalar f1 = xest*delta;
65 scalar f2 = x[iest - k1 - 1]*delta;
67 for (
register label j=0; j<n; j++)
78 for (
register label j=0; j<n; j++)