222     PetscCall(PetscFree(rk->YdotRHS_fast));  in TSReset_RK_MultirateSplit()
299     PetscCall(VecMAXPY(Xfast, s, w, rk->YdotRHS_fast));  in TSEvaluateStep_RK_MultirateSplit()
312   Vec             *YdotRHS = rk->YdotRHS, *YdotRHS_fast = rk->YdotRHS_fast;  in TSStepRefine_RK_MultirateSplit()  local
323     for (i = 0; i < s; i++) PetscCall(VecGetSubVector(YdotRHS[i], rk->is_fast, &YdotRHS_fast[i]));  in TSStepRefine_RK_MultirateSplit()
339       PetscCall(VecMAXPY(Yfast, i, w, YdotRHS_fast));  in TSStepRefine_RK_MultirateSplit()
343 …eRHSFunction(subts_fast, t + k * h * rk->dtratio + h / rk->dtratio * c[i], Y[i], YdotRHS_fast[i]));  in TSStepRefine_RK_MultirateSplit()
349 … for (i = 0; i < s; i++) PetscCall(VecRestoreSubVector(YdotRHS[i], rk->is_fast, &YdotRHS_fast[i]));  in TSStepRefine_RK_MultirateSplit()
369   Vec             *YdotRHS_fast = rk->YdotRHS_fast, *YdotRHS_slow = rk->YdotRHS_slow;  in TSStep_RK_MultirateSplit()  local
381     PetscCall(VecGetSubVector(YdotRHS[i], rk->is_fast, &YdotRHS_fast[i]));  in TSStep_RK_MultirateSplit()
392     PetscCall(VecMAXPY(Yfast, i, w, YdotRHS_fast));  in TSStep_RK_MultirateSplit()
398     PetscCall(TSComputeRHSFunction(rk->subts_fast, t + h * c[i], Y[i], YdotRHS_fast[i]));  in TSStep_RK_MultirateSplit()
406     PetscCall(VecRestoreSubVector(YdotRHS[i], rk->is_fast, &YdotRHS_fast[i]));  in TSStep_RK_MultirateSplit()
435     PetscCall(PetscMalloc1(rk->tableau->s, &currentlevelrk->YdotRHS_fast));  in TSSetUp_RK_MultirateSplit()