MultiRate Runge-Kutta, Order 2, Constantinescu

Solver/TESTCASES/ConvergenceMultirate.lua

@@ -82,6 +82,7 @@ for i=1,nRefSteps do
 end
 
 multirateRK=dgRungeKuttaMultirate43(GC,law)
+--multirateRK=dgRungeKuttaMultirate22(GC,law)
 integrator=dgFunctionIntegrator(functionLua(1, 'L2Norm', {'Solution'}):getName())
 
 oldnorm=1

Solver/TESTCASES/MultirateAdvection.lua

@@ -86,6 +86,7 @@ dt=dt
 print('DT=',dt)
 RK=dgRungeKutta()
 multirateRK=dgRungeKuttaMultirate43(GC,law)
+--multirateRK=dgRungeKuttaMultirate22(GC,law)
 norm1=solution:norm()
 norm2=solution2:norm()
 print('Norm: ',norm1,norm2)