649   PetscScalar     *wf = mprk->work_fast, *wsb = mprk->work_slowbuffer;  in TSStep_MPRK()  local
660     for (j = 0; j < i; j++) wsb[j] = h * Asb[i * s + j];  in TSStep_MPRK()
663     PetscCall(VecMAXPY(Yslowbuffer, i, wsb, mprk->YdotRHS_slowbuffer));  in TSStep_MPRK()
670       PetscCall(VecMAXPY(Yslow, i, wsb, mprk->YdotRHS_slow));  in TSStep_MPRK()
723   PetscScalar *wf = mprk->work_fast, *ws = mprk->work_slow, *wsb = mprk->work_slowbuffer;  in TSEvaluateStep_MPRKSPLIT()  local
745 …if (tab->rmb[j]) wsb[computedstages - tab->sbase + (tab->rmb[j] - 1) % tab->sbase] += h * tab->bsb…  in TSEvaluateStep_MPRKSPLIT()
746       else wsb[computedstages++] = h * tab->bsb[j];  in TSEvaluateStep_MPRKSPLIT()
749     PetscCall(VecMAXPY(Xslowbuffer, computedstages, wsb, mprk->YdotRHS_slowbuffer));  in TSEvaluateStep_MPRKSPLIT()
753     for (j = 0; j < s; j++) wsb[j] = h * tab->bsb[j];  in TSEvaluateStep_MPRKSPLIT()
755     PetscCall(VecMAXPY(Xslowbuffer, s, wsb, mprk->YdotRHS_slowbuffer));  in TSEvaluateStep_MPRKSPLIT()
794   PetscScalar     *wf = mprk->work_fast, *ws = mprk->work_slow, *wsb = mprk->work_slowbuffer;  in TSStep_MPRKSPLIT()  local
804     for (j = 0; j < i; j++) wsb[j] = h * Asb[i * s + j];  in TSStep_MPRKSPLIT()
816 …       if (tab->rmb[j]) wm[computedstages - tab->sbase + (tab->rmb[j] - 1) % tab->sbase] += wsb[j];  in TSStep_MPRKSPLIT()
817           else wm[computedstages++] = wsb[j];  in TSStep_MPRKSPLIT()
825       PetscCall(VecMAXPY(Yslowbuffer, i, wsb, YdotRHS_slowbuffer));  in TSStep_MPRKSPLIT()
838           if (tab->rsb[j]) ws[tab->rsb[j] - 1] += wsb[j];  in TSStep_MPRKSPLIT()
839           else ws[computedstages++] = wsb[j];  in TSStep_MPRKSPLIT()