6 static inline PetscReal DMatrix_Determinant_2x2_Internal(const PetscReal inmat[2 * 2])  in DMatrix_Determinant_2x2_Internal()
8   return inmat[0] * inmat[3] - inmat[1] * inmat[2];  in DMatrix_Determinant_2x2_Internal()
11 static inline PetscErrorCode DMatrix_Invert_2x2_Internal(const PetscReal *inmat, PetscReal *outmat,…  in DMatrix_Invert_2x2_Internal()  argument
13   PetscReal det = DMatrix_Determinant_2x2_Internal(inmat);  in DMatrix_Invert_2x2_Internal()
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()
24 static inline PetscReal DMatrix_Determinant_3x3_Internal(const PetscReal inmat[3 * 3])  in DMatrix_Determinant_3x3_Internal()
26 …inmat[0] * (inmat[8] * inmat[4] - inmat[7] * inmat[5]) - inmat[3] * (inmat[8] * inmat[1] - inmat[7…  in DMatrix_Determinant_3x3_Internal()
29 static inline PetscErrorCode DMatrix_Invert_3x3_Internal(const PetscReal *inmat, PetscReal *outmat,…  in DMatrix_Invert_3x3_Internal()  argument
31   PetscReal det = DMatrix_Determinant_3x3_Internal(inmat);  in DMatrix_Invert_3x3_Internal()
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()
47 static inline PetscReal DMatrix_Determinant_4x4_Internal(PetscReal inmat[4 * 4])  in DMatrix_Determinant_4x4_Internal()
49 …inmat[0 + 0 * 4] * (inmat[1 + 1 * 4] * (inmat[2 + 2 * 4] * inmat[3 + 3 * 4] - inmat[2 + 3 * 4] * i…  in DMatrix_Determinant_4x4_Internal()
52 static inline PETSC_UNUSED PetscErrorCode DMatrix_Invert_4x4_Internal(PetscReal *inmat, PetscReal *…  in DMatrix_Invert_4x4_Internal()  argument
54   PetscReal det = DMatrix_Determinant_4x4_Internal(inmat);  in DMatrix_Invert_4x4_Internal()
56 …inmat[5] * inmat[10] * inmat[15] + inmat[6] * inmat[11] * inmat[13] + inmat[7] * inmat[9] * inmat[…  in DMatrix_Invert_4x4_Internal()
57 …inmat[1] * inmat[11] * inmat[14] + inmat[2] * inmat[9] * inmat[15] + inmat[3] * inmat[10] * inmat[…  in DMatrix_Invert_4x4_Internal()
58 …inmat[1] * inmat[6] * inmat[15] + inmat[2] * inmat[7] * inmat[13] + inmat[3] * inmat[5] * inmat[14…  in DMatrix_Invert_4x4_Internal()
59 …inmat[1] * inmat[7] * inmat[10] + inmat[2] * inmat[5] * inmat[11] + inmat[3] * inmat[6] * inmat[9]…  in DMatrix_Invert_4x4_Internal()
60 …inmat[4] * inmat[11] * inmat[14] + inmat[6] * inmat[8] * inmat[15] + inmat[7] * inmat[10] * inmat[…  in DMatrix_Invert_4x4_Internal()
61 …inmat[0] * inmat[10] * inmat[15] + inmat[2] * inmat[11] * inmat[12] + inmat[3] * inmat[8] * inmat[…  in DMatrix_Invert_4x4_Internal()
62 …inmat[0] * inmat[7] * inmat[14] + inmat[2] * inmat[4] * inmat[15] + inmat[3] * inmat[6] * inmat[12…  in DMatrix_Invert_4x4_Internal()
63 …inmat[0] * inmat[6] * inmat[11] + inmat[2] * inmat[7] * inmat[8] + inmat[3] * inmat[4] * inmat[10]…  in DMatrix_Invert_4x4_Internal()
64 …inmat[4] * inmat[9] * inmat[15] + inmat[5] * inmat[11] * inmat[12] + inmat[7] * inmat[8] * inmat[1…  in DMatrix_Invert_4x4_Internal()
65 …inmat[0] * inmat[11] * inmat[13] + inmat[1] * inmat[8] * inmat[15] + inmat[3] * inmat[9] * inmat[1…  in DMatrix_Invert_4x4_Internal()
66 …inmat[0] * inmat[5] * inmat[15] + inmat[1] * inmat[7] * inmat[12] + inmat[3] * inmat[4] * inmat[13…  in DMatrix_Invert_4x4_Internal()
67 …inmat[0] * inmat[7] * inmat[9] + inmat[1] * inmat[4] * inmat[11] + inmat[3] * inmat[5] * inmat[8] …  in DMatrix_Invert_4x4_Internal()
68 …inmat[4] * inmat[10] * inmat[13] + inmat[5] * inmat[8] * inmat[14] + inmat[6] * inmat[9] * inmat[1…  in DMatrix_Invert_4x4_Internal()
69 …inmat[0] * inmat[9] * inmat[14] + inmat[1] * inmat[10] * inmat[12] + inmat[2] * inmat[8] * inmat[1…  in DMatrix_Invert_4x4_Internal()
70 …inmat[0] * inmat[6] * inmat[13] + inmat[1] * inmat[4] * inmat[14] + inmat[2] * inmat[5] * inmat[12…  in DMatrix_Invert_4x4_Internal()
71 …inmat[0] * inmat[5] * inmat[10] + inmat[1] * inmat[6] * inmat[8] + inmat[2] * inmat[4] * inmat[9] …  in DMatrix_Invert_4x4_Internal()