38 RK::safety=0.9, RK::pGrow=-0.2, RK::pShrink=-0.25, RK::errCon=1.89e-4;
41 RK::a2 = 0.2, RK::a3 = 0.3, RK::a4 = 0.6, RK::a5 = 1.0, RK::a6 = 0.875,
42 RK::b21 = 0.2, RK::b31 = 3.0/40.0, RK::b32 = 9.0/40.0,
43 RK::b41 = 0.3, RK::b42 = -0.9, RK::b43 = 1.2,
44 RK::b51 = -11.0/54.0, RK::b52 = 2.5, RK::b53 = -70.0/27.0,
46 RK::b61 = 1631.0/55296.0, RK::b62 = 175.0/512.0, RK::b63 = 575.0/13824.0,
47 RK::b64 = 44275.0/110592.0, RK::b65 = 253.0/4096.0,
48 RK::c1 = 37.0/378.0, RK::c3 = 250.0/621.0,
49 RK::c4 = 125.0/594.0, RK::c6 = 512.0/1771.0,
50 RK::dc1 = RK::c1 - 2825.0/27648.0, RK::dc3 = RK::c3 - 18575.0/48384.0,
51 RK::dc4 = RK::c4 - 13525.0/55296.0, RK::dc5 = -277.00/14336.0,
52 RK::dc6 = RK::c6 - 0.25;
87 yTemp_[i] = y[i] + b21*h*dydx[i];
94 yTemp_[i] = y[i] + h*(b31*dydx[i] + b32*ak2_[i]);
101 yTemp_[i] = y[i] + h*(b41*dydx[i] + b42*ak2_[i] + b43*ak3_[i]);
109 + h*(b51*dydx[i] + b52*ak2_[i] + b53*ak3_[i] + b54*ak4_[i]);
119 b61*dydx[i] + b62*ak2_[i] + b63*ak3_[i]
120 + b64*ak4_[i] + b65*ak5_[i]
129 + h*(c1*dydx[i] + c3*ak3_[i] + c4*ak4_[i] + c6*ak6_[i]);
137 dc1*dydx[i] + dc3*ak3_[i] + dc4*ak4_[i]
138 + dc5*ak5_[i] + dc6*ak6_[i]
162 solve(ode, x, y, dydx, h, yTemp2_, yErr_);
165 for (
register label i=0; i<n_; i++)
167 maxErr =
max(maxErr,
mag(yErr_[i]/yScale[i]));
177 scalar hTemp = safety*h*
pow(maxErr, pShrink);
178 h = (h >= 0.0 ?
max(hTemp, 0.1*h) :
min(hTemp, 0.1*h));
184 <<
"stepsize underflow"
196 hNext = safety*h*
pow(maxErr, pGrow);