278   Mat            J, Jpre, quadJ = NULL, quadJp = NULL;  in TSAdjointStepBEuler_Private()  local
296     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSAdjointStepBEuler_Private()
304   if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStepBEuler_Private()
309     if (quadJ) {  in TSAdjointStepBEuler_Private()
310       PetscCall(MatDenseGetColumn(quadJ, nadj, &xarr));  in TSAdjointStepBEuler_Private()
314       PetscCall(MatDenseRestoreColumn(quadJ, &xarr));  in TSAdjointStepBEuler_Private()
423   Mat          J, Jpre, quadJ = NULL, quadJp = NULL;  in TSAdjointStep_Theta()  local
438     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSAdjointStep_Theta()
455       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSAdjointStep_Theta()
457       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSAdjointStep_Theta()
464     if (quadJ) {  in TSAdjointStep_Theta()
465       PetscCall(MatDenseGetColumn(quadJ, nadj, &xarr));  in TSAdjointStep_Theta()
469       PetscCall(MatDenseRestoreColumn(quadJ, &xarr));  in TSAdjointStep_Theta()
530     if (quadts) PetscCall(TSComputeRHSJacobian(quadts, adjoint_ptime, th->X0, quadJ, NULL));  in TSAdjointStep_Theta()
533       if (quadJ) {  in TSAdjointStep_Theta()
534         PetscCall(MatDenseGetColumn(quadJ, nadj, &xarr));  in TSAdjointStep_Theta()
538         PetscCall(MatDenseRestoreColumn(quadJ, &xarr));  in TSAdjointStep_Theta()
640     if (quadts) PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSAdjointStep_Theta()
644       if (quadJ) {  in TSAdjointStep_Theta()
645         PetscCall(MatDenseGetColumn(quadJ, nadj, &xarr));  in TSAdjointStep_Theta()
649         PetscCall(MatDenseRestoreColumn(quadJ, &xarr));  in TSAdjointStep_Theta()
746   Mat          J, Jpre, quadJ = NULL, quadJp = NULL;  in TSForwardStep_Theta()  local
758     PetscCall(TSGetRHSJacobian(quadts, &quadJ, NULL, NULL, NULL));  in TSForwardStep_Theta()
810       PetscCall(TSComputeRHSJacobian(quadts, th->ptime0, th->X0, quadJ, NULL));  in TSForwardStep_Theta()
812 …PetscCall(MatTransposeMatMult(ts->mat_sensip, quadJ, MAT_REUSE_MATRIX, PETSC_DETERMINE, &th->MatIn…  in TSForwardStep_Theta()
848       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, th->X, quadJ, NULL));  in TSForwardStep_Theta()
850 …PetscCall(MatTransposeMatMult(ts->mat_sensip, quadJ, MAT_REUSE_MATRIX, PETSC_DETERMINE, &th->MatIn…  in TSForwardStep_Theta()
855       PetscCall(TSComputeRHSJacobian(quadts, th->stage_time, ts->vec_sol, quadJ, NULL));  in TSForwardStep_Theta()
857 …PetscCall(MatTransposeMatMult(ts->mat_sensip, quadJ, MAT_REUSE_MATRIX, PETSC_DETERMINE, &th->MatIn…  in TSForwardStep_Theta()