50 std::vector<double>& rCurrentYValues,
51 std::vector<double>& rNextYValues)
80 const double mu1_tilde = 0.256134735558604;
81 const double mu2 = 1.952097590002976;
83 const double mu2_tilde =0.500000000000000;
88 std::vector<double>& w0 = rCurrentYValues;
89 std::vector<double> w1(num_equations);
90 std::vector<double>& w2 = rNextYValues;
91 std::vector<double> F0(num_equations);
92 std::vector<double>& F1 = rNextYValues;
97 for (
unsigned i=0; i<num_equations; i++)
99 w1[i] = w0[i] + mu1_tilde * timeStep * F0[i];
104 for (
unsigned i=0; i<num_equations; i++)
106 w2[i] = (1-mu2) * w0[i] + mu2 * w1[i] + mu2_tilde * timeStep * F1[i];