278   Mat            J, Jpre, quadJ = NULL, quadJp = NULL;  in TSAdjointStepBEuler_Private()  local
297     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSAdjointStepBEuler_Private()
381     if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStepBEuler_Private()
392       if (quadJp) {  in TSAdjointStepBEuler_Private()
393         PetscCall(MatDenseGetColumn(quadJp, nadj, &xarr));  in TSAdjointStepBEuler_Private()
397         PetscCall(MatDenseRestoreColumn(quadJp, &xarr));  in TSAdjointStepBEuler_Private()
423   Mat          J, Jpre, quadJ = NULL, quadJp = NULL;  in TSAdjointStep_Theta()  local
439     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSAdjointStep_Theta()
570       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSAdjointStep_Theta()
574         if (quadJp) {  in TSAdjointStep_Theta()
575           PetscCall(MatDenseGetColumn(quadJp, nadj, &xarr));  in TSAdjointStep_Theta()
579           PetscCall(MatDenseRestoreColumn(quadJp, &xarr));  in TSAdjointStep_Theta()
605       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, adjoint_ptime, th->X0, quadJp));  in TSAdjointStep_Theta()
609         if (quadJp) {  in TSAdjointStep_Theta()
610           PetscCall(MatDenseGetColumn(quadJp, nadj, &xarr));  in TSAdjointStep_Theta()
614           PetscCall(MatDenseRestoreColumn(quadJp, &xarr));  in TSAdjointStep_Theta()
656       if (quadts) PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, th->X, quadJp));  in TSAdjointStep_Theta()
660         if (quadJp) {  in TSAdjointStep_Theta()
661           PetscCall(MatDenseGetColumn(quadJp, nadj, &xarr));  in TSAdjointStep_Theta()
665           PetscCall(MatDenseRestoreColumn(quadJp, &xarr));  in TSAdjointStep_Theta()
746   Mat          J, Jpre, quadJ = NULL, quadJp = NULL;  in TSForwardStep_Theta()  local
759     PetscCall(TSGetRHSJacobianP(quadts, &quadJp, NULL, NULL));  in TSForwardStep_Theta()
811       PetscCall(TSComputeRHSJacobianP(quadts, th->ptime0, th->X0, quadJp));  in TSForwardStep_Theta()
813       PetscCall(MatAXPY(th->MatIntegralSensipTemp, 1, quadJp, SAME_NONZERO_PATTERN));  in TSForwardStep_Theta()
849       PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, th->X, quadJp));  in TSForwardStep_Theta()
851       PetscCall(MatAXPY(th->MatIntegralSensipTemp, 1, quadJp, SAME_NONZERO_PATTERN));  in TSForwardStep_Theta()
856       PetscCall(TSComputeRHSJacobianP(quadts, th->stage_time, ts->vec_sol, quadJp));  in TSForwardStep_Theta()
858       PetscCall(MatAXPY(th->MatIntegralSensipTemp, 1, quadJp, SAME_NONZERO_PATTERN));  in TSForwardStep_Theta()