30 void Foam::SIBS::SIMPR
43 scalar
h = deltaX/nSteps;
46 for (
register label i=0; i<
n_; i++)
48 for (
register label j=0; j<
n_; j++)
50 a[i][j] = -h*dfdy[i][j];
58 for (
register label i=0; i<
n_; i++)
60 yEnd[i] = h*(dydx[i] + h*dfdx[i]);
68 for (
register label i=0; i<
n_; i++)
70 ytemp[i] = y[i] + del[i];
73 scalar x = xStart +
h;
75 ode.derivatives(x, ytemp, yEnd);
77 for (
register label nn=2; nn<=nSteps; nn++)
79 for (
register label i=0; i<
n_; i++)
81 yEnd[i] = h*yEnd[i] - del[i];
86 for (
register label i=0; i<
n_; i++)
88 ytemp[i] += (del[i] += 2.0*yEnd[i]);
93 ode.derivatives(x, ytemp, yEnd);
95 for (
register label i=0; i<
n_; i++)
97 yEnd[i] = h*yEnd[i] - del[i];
102 for (
register label i=0; i<
n_; i++)