39   Vec         *YdotI; /* Work vectors to hold the residual evaluation */  member
273   Vec         *YdotI = irk->YdotI;  in TSEvaluateStep_IRK()  local
293   PetscCall(VecMAXPY(U, irk->nstages, w, YdotI));  in TSEvaluateStep_IRK()
335       PetscCall(VecZeroEntries(irk->YdotI[i]));  in TSStep_IRK()
336 …for (j = 0; j < nstages; j++) PetscCall(VecAXPY(irk->YdotI[i], A_inv[i + j * nstages] / ts->time_s…  in TSStep_IRK()
337       PetscCall(VecAXPY(irk->YdotI[i], -A_inv_rowsum[i] / ts->time_step, irk->U));  in TSStep_IRK()
394   PetscCall(VecMAXPY(U, nstages, bt, irk->YdotI));  in TSInterpolate_IRK()
420   PetscCall(VecDestroyVecs(irk->nstages, &irk->YdotI));  in TSReset_IRK()
464   Vec                U, *YdotI = irk->YdotI, Ydot = irk->Ydot, *Y = irk->Y;  in SNESTSFormFunction_IRK()  local
478 …PetscCall(TSComputeIFunction(ts, ts->ptime + ts->time_step * c[i], Y[i], Ydot, YdotI[i], PETSC_FAL…  in SNESTSFormFunction_IRK()
480   PetscCall(VecStrideScatterAll(YdotI, FC, INSERT_VALUES));  in SNESTSFormFunction_IRK()
580   if (!irk->YdotI) PetscCall(VecDuplicateVecs(ts->vec_sol, irk->nstages, &irk->YdotI));  in TSSetUp_IRK()