13   PetscReal det = DMatrix_Determinant_2x2_Internal(inmat);  in DMatrix_Invert_2x2_Internal()  local
15     outmat[0] = inmat[3] / det;  in DMatrix_Invert_2x2_Internal()
16     outmat[1] = -inmat[1] / det;  in DMatrix_Invert_2x2_Internal()
17     outmat[2] = -inmat[2] / det;  in DMatrix_Invert_2x2_Internal()
18     outmat[3] = inmat[0] / det;  in DMatrix_Invert_2x2_Internal()
20   if (determinant) *determinant = det;  in DMatrix_Invert_2x2_Internal()
31   PetscReal det = DMatrix_Determinant_3x3_Internal(inmat);  in DMatrix_Invert_3x3_Internal()  local
33     outmat[0] = (inmat[8] * inmat[4] - inmat[7] * inmat[5]) / det;  in DMatrix_Invert_3x3_Internal()
34     outmat[1] = -(inmat[8] * inmat[1] - inmat[7] * inmat[2]) / det;  in DMatrix_Invert_3x3_Internal()
35     outmat[2] = (inmat[5] * inmat[1] - inmat[4] * inmat[2]) / det;  in DMatrix_Invert_3x3_Internal()
36     outmat[3] = -(inmat[8] * inmat[3] - inmat[6] * inmat[5]) / det;  in DMatrix_Invert_3x3_Internal()
37     outmat[4] = (inmat[8] * inmat[0] - inmat[6] * inmat[2]) / det;  in DMatrix_Invert_3x3_Internal()
38     outmat[5] = -(inmat[5] * inmat[0] - inmat[3] * inmat[2]) / det;  in DMatrix_Invert_3x3_Internal()
39     outmat[6] = (inmat[7] * inmat[3] - inmat[6] * inmat[4]) / det;  in DMatrix_Invert_3x3_Internal()
40     outmat[7] = -(inmat[7] * inmat[0] - inmat[6] * inmat[1]) / det;  in DMatrix_Invert_3x3_Internal()
41     outmat[8] = (inmat[4] * inmat[0] - inmat[3] * inmat[1]) / det;  in DMatrix_Invert_3x3_Internal()
43   if (determinant) *determinant = det;  in DMatrix_Invert_3x3_Internal()
54   PetscReal det = DMatrix_Determinant_4x4_Internal(inmat);  in DMatrix_Invert_4x4_Internal()  local
56 … inmat[11] * inmat[14] - inmat[6] * inmat[9] * inmat[15] - inmat[7] * inmat[10] * inmat[13]) / det;  in DMatrix_Invert_4x4_Internal()
57 … inmat[10] * inmat[15] - inmat[2] * inmat[11] * inmat[13] - inmat[3] * inmat[9] * inmat[14]) / det;  in DMatrix_Invert_4x4_Internal()
58 … * inmat[7] * inmat[14] - inmat[2] * inmat[5] * inmat[15] - inmat[3] * inmat[6] * inmat[13]) / det;  in DMatrix_Invert_4x4_Internal()
59 …] * inmat[6] * inmat[11] - inmat[2] * inmat[7] * inmat[9] - inmat[3] * inmat[5] * inmat[10]) / det;  in DMatrix_Invert_4x4_Internal()
60 … inmat[10] * inmat[15] - inmat[6] * inmat[11] * inmat[12] - inmat[7] * inmat[8] * inmat[14]) / det;  in DMatrix_Invert_4x4_Internal()
61 … inmat[11] * inmat[14] - inmat[2] * inmat[8] * inmat[15] - inmat[3] * inmat[10] * inmat[12]) / det;  in DMatrix_Invert_4x4_Internal()
62 … * inmat[6] * inmat[15] - inmat[2] * inmat[7] * inmat[12] - inmat[3] * inmat[4] * inmat[14]) / det;  in DMatrix_Invert_4x4_Internal()
63 …] * inmat[7] * inmat[10] - inmat[2] * inmat[4] * inmat[11] - inmat[3] * inmat[6] * inmat[8]) / det;  in DMatrix_Invert_4x4_Internal()
64 …* inmat[11] * inmat[13] - inmat[5] * inmat[8] * inmat[15] - inmat[7] * inmat[9] * inmat[12]) / det;  in DMatrix_Invert_4x4_Internal()
65 …* inmat[9] * inmat[15] - inmat[1] * inmat[11] * inmat[12] - inmat[3] * inmat[8] * inmat[13]) / det;  in DMatrix_Invert_4x4_Internal()
66 … * inmat[7] * inmat[13] - inmat[1] * inmat[4] * inmat[15] - inmat[3] * inmat[5] * inmat[12]) / det;  in DMatrix_Invert_4x4_Internal()
67 …0] * inmat[5] * inmat[11] - inmat[1] * inmat[7] * inmat[8] - inmat[3] * inmat[4] * inmat[9]) / det;  in DMatrix_Invert_4x4_Internal()
68 …* inmat[9] * inmat[14] - inmat[5] * inmat[10] * inmat[12] - inmat[6] * inmat[8] * inmat[13]) / det;  in DMatrix_Invert_4x4_Internal()
69 …* inmat[10] * inmat[13] - inmat[1] * inmat[8] * inmat[14] - inmat[2] * inmat[9] * inmat[12]) / det;  in DMatrix_Invert_4x4_Internal()
70 … * inmat[5] * inmat[14] - inmat[1] * inmat[6] * inmat[12] - inmat[2] * inmat[4] * inmat[13]) / det;  in DMatrix_Invert_4x4_Internal()
71 …0] * inmat[6] * inmat[9] - inmat[1] * inmat[4] * inmat[10] - inmat[2] * inmat[5] * inmat[8]) / det;  in DMatrix_Invert_4x4_Internal()
73   if (determinant) *determinant = det;  in DMatrix_Invert_4x4_Internal()