39   PetscReal *GammaInv;          /* Inverse of Gamma, used for transformed variables */  member
763     PetscCall(PetscFree5(t->At, t->bt, t->GammaInv, t->GammaZeroDiag, t->GammaExplicitCorr));  in TSRosWRegisterDestroy()
835   PetscScalar    *GammaInv;  in TSRosWRegister()  local
851 …PetscCall(PetscMalloc5(s * s, &t->At, s, &t->bt, s * s, &t->GammaInv, s, &t->GammaZeroDiag, s * s,…  in TSRosWRegister()
868 …PetscCall(PetscMalloc1(s * s, &GammaInv)); /* Need to use Scalar for inverse, then convert back to…  in TSRosWRegister()
869   for (i = 0; i < s * s; i++) GammaInv[i] = Gamma[i];  in TSRosWRegister()
872       GammaInv[i * s + i] = 1.0;  in TSRosWRegister()
881     GammaInv[0] = 1. / GammaInv[0];  in TSRosWRegister()
884     PetscCall(PetscKernel_A_gets_inverse_A_2(GammaInv, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
887     PetscCall(PetscKernel_A_gets_inverse_A_3(GammaInv, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
890     PetscCall(PetscKernel_A_gets_inverse_A_4(GammaInv, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
895     PetscCall(PetscKernel_A_gets_inverse_A_5(GammaInv, ipvt5, work5, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
899     PetscCall(PetscKernel_A_gets_inverse_A_6(GammaInv, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
902     PetscCall(PetscKernel_A_gets_inverse_A_7(GammaInv, 0, PETSC_FALSE, NULL));  in TSRosWRegister()
907   for (i = 0; i < s * s; i++) t->GammaInv[i] = PetscRealPart(GammaInv[i]);  in TSRosWRegister()
908   PetscCall(PetscFree(GammaInv));  in TSRosWRegister()
912 …    t->GammaExplicitCorr[i * s + k] = (t->GammaExplicitCorr[i * s + k]) * (t->GammaInv[k * s + k]);  in TSRosWRegister()
913 …+) t->GammaExplicitCorr[i * s + k] += (t->GammaExplicitCorr[i * s + j]) * (t->GammaInv[j * s + k]);  in TSRosWRegister()
920       for (k = 0; k < s; k++) t->At[i * s + j] += t->A[i * s + k] * t->GammaInv[k * s + j];  in TSRosWRegister()
923     for (j = 0; j < s; j++) t->bt[i] += t->b[j] * t->GammaInv[j * s + i];  in TSRosWRegister()
926       for (j = 0; j < s; j++) t->bembedt[i] += t->bembed[j] * t->GammaInv[j * s + i];  in TSRosWRegister()
1130   const PetscReal *At = tab->At, *Gamma = tab->Gamma, *ASum = tab->ASum, *GammaInv = tab->GammaInv;  in TSStep_RosW()  local
1164       for (j = 0; j < i; j++) w[j] = 1. / h * GammaInv[i * s + j];  in TSStep_RosW()
1248   const PetscReal *GammaInv = ros->tableau->GammaInv;  in TSInterpolate_RosW()  local
1280     for (i = j; i < s; i++) w[j] += bt[i] * GammaInv[i * s + j];  in TSInterpolate_RosW()