43 constexpr
G4double inv_STEPFAC1 = 1.0 / STEPFAC1;
44 constexpr
G4double inv_STEPFAC4 = 1.0 / STEPFAC4;
49 : fnvar(nvar), m_eps_rel(eps_rel), m_midpoint(equation,nvar),
50 m_last_step_rejected(
false), m_first(
true), m_dt_last(0.0), m_max_dt(max_dt)
92 return step_result::fail;
142 / static_cast<G4double>(
m_cost[k]);
174 / static_cast<G4double>(
m_cost[k]);
230 return reject ? step_result::fail : step_result::success;
244 for(
G4int j = k - 1 ; j > 0; --j)
263 const G4double expo = 1.0 / (2 * k + 1);
264 const G4double facmin = std::pow(STEPFAC3, expo);
274 fac = STEPFAC2 * std::pow(error * inv_STEPFAC1 , -expo);
297 else if( (
work[k] < KFAC2 *
work[k-1])
331 return error * e2 * e2 > d *
d;
337 return error * e * e > d *
d;