51   Vec         *YdotRHS_slow;         /* Function evaluations by slow tableau for slow components  */  member
644 …otRHS = mprk->YdotRHS, *YdotRHS_fast = mprk->YdotRHS_fast, *YdotRHS_slow = mprk->YdotRHS_slow, *Yd…  in TSStep_MPRK()  local
668 …   for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_slow, &YdotRHS_slow[j]));  in TSStep_MPRK()
670       PetscCall(VecMAXPY(Yslow, i, wsb, mprk->YdotRHS_slow));  in TSStep_MPRK()
672 …or (j = 0; j < i; j++) PetscCall(VecRestoreSubVector(YdotRHS[j], mprk->is_slow, &YdotRHS_slow[j]));  in TSStep_MPRK()
738     PetscCall(VecMAXPY(Xslow, computedstages, ws, mprk->YdotRHS_slow));  in TSEvaluateStep_MPRKSPLIT()
790 …         *Y = mprk->Y, *YdotRHS_fast = mprk->YdotRHS_fast, *YdotRHS_slow = mprk->YdotRHS_slow, *Yd…  in TSStep_MPRKSPLIT()  local
842         PetscCall(VecMAXPY(Yslow, computedstages, ws, YdotRHS_slow));  in TSStep_MPRKSPLIT()
845 …PetscCall(TSComputeRHSFunction(mprk->subts_slow, t + h * csb[i], Y[i], YdotRHS_slow[computedstages…  in TSStep_MPRKSPLIT()
915     PetscCall(VecDestroyVecs(tab->s, &mprk->YdotRHS_slow));  in TSMPRKTableauReset()
921     if (mprk->is_slow) PetscCall(PetscFree(mprk->YdotRHS_slow));  in TSMPRKTableauReset()
967   Vec         YdotRHS_slow, YdotRHS_slowbuffer, YdotRHS_medium, YdotRHS_mediumbuffer, YdotRHS_fast;  in TSMPRKTableauSetUp()  local
981       PetscCall(VecGetSubVector(ts->vec_sol, mprk->is_slow, &YdotRHS_slow));  in TSMPRKTableauSetUp()
982       PetscCall(VecDuplicateVecs(YdotRHS_slow, tab->s, &mprk->YdotRHS_slow));  in TSMPRKTableauSetUp()
983       PetscCall(VecRestoreSubVector(ts->vec_sol, mprk->is_slow, &YdotRHS_slow));  in TSMPRKTableauSetUp()
1003     if (mprk->is_slow) PetscCall(PetscMalloc1(tab->s, &mprk->YdotRHS_slow));  in TSMPRKTableauSetUp()