65   IS           is_slow, is_slowbuffer, is_medium, is_mediumbuffer, is_fast;  member
689       if (mprk->is_medium) {  in TSStep_MPRK()
690 …for (j = 0; j < i; j++) PetscCall(VecGetSubVector(YdotRHS[j], mprk->is_medium, &YdotRHS_medium[j])…  in TSStep_MPRK()
691         PetscCall(VecGetSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRK()
693         PetscCall(VecRestoreSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRK()
694 …for (j = 0; j < i; j++) PetscCall(VecRestoreSubVector(YdotRHS[j], mprk->is_medium, &YdotRHS_medium…  in TSStep_MPRK()
742   if (tab->np == 3 && mprk->is_medium) {  in TSEvaluateStep_MPRKSPLIT()
762     if (mprk->is_medium) {  in TSEvaluateStep_MPRKSPLIT()
768       PetscCall(VecGetSubVector(X, mprk->is_medium, &Xmedium));  in TSEvaluateStep_MPRKSPLIT()
770       PetscCall(VecRestoreSubVector(X, mprk->is_medium, &Xmedium));  in TSEvaluateStep_MPRKSPLIT()
807     if (tab->np == 3 && mprk->is_medium) {  in TSStep_MPRKSPLIT()
867       if (mprk->is_medium) {  in TSStep_MPRKSPLIT()
869           PetscCall(VecGetSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRKSPLIT()
870           PetscCall(VecISCopy(Y[tab->rmb[i] - 1], mprk->is_medium, SCATTER_REVERSE, Ymedium));  in TSStep_MPRKSPLIT()
871           PetscCall(VecRestoreSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRKSPLIT()
878           PetscCall(VecGetSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRKSPLIT()
880           PetscCall(VecRestoreSubVector(Y[i], mprk->is_medium, &Ymedium));  in TSStep_MPRKSPLIT()
890 …if (!(tab->np == 3 && mprk->is_medium)) PetscCall(TSComputeRHSFunction(mprk->subts_slowbuffer, t +…  in TSStep_MPRKSPLIT()
924       if (mprk->is_medium) PetscCall(PetscFree(mprk->YdotRHS_medium));  in TSMPRKTableauReset()
974     if (mprk->is_medium) PetscCall(PetscMalloc1(tab->s, &mprk->work_medium));  in TSMPRKTableauSetUp()
989       if (mprk->is_medium) {  in TSMPRKTableauSetUp()
990         PetscCall(VecGetSubVector(ts->vec_sol, mprk->is_medium, &YdotRHS_medium));  in TSMPRKTableauSetUp()
992         PetscCall(VecRestoreSubVector(ts->vec_sol, mprk->is_medium, &YdotRHS_medium));  in TSMPRKTableauSetUp()
1006       if (mprk->is_medium) PetscCall(PetscMalloc1(tab->s, &mprk->YdotRHS_medium));  in TSMPRKTableauSetUp()
1026     PetscCall(TSRHSSplitGetIS(ts, "medium", &mprk->is_medium));  in TSSetUp_MPRK()
1027 …PetscCheck(mprk->is_medium, PetscObjectComm((PetscObject)ts), PETSC_ERR_USER, "Must set up RHSSpli…  in TSSetUp_MPRK()
1030       mprk->is_mediumbuffer = mprk->is_medium;  in TSSetUp_MPRK()
1031       mprk->is_medium       = NULL;  in TSSetUp_MPRK()