14 PetscErrorCode PetscKernel_A_gets_inverse_A_3(MatScalar *a, PetscReal shift, PetscBool allowzeropiv…  in PetscKernel_A_gets_inverse_A_3()  argument
23 …shift = .333 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[4]) + PetscAbsScalar(a[8]…  in PetscKernel_A_gets_inverse_A_3()
26   a -= 4;  in PetscKernel_A_gets_inverse_A_3()
35     aa   = &a[k4];  in PetscKernel_A_gets_inverse_A_3()
48     if (a[l + k3] == 0.0) {  in PetscKernel_A_gets_inverse_A_3()
55         a[l + k3] = shift;  in PetscKernel_A_gets_inverse_A_3()
61       stmp      = a[l + k3];  in PetscKernel_A_gets_inverse_A_3()
62       a[l + k3] = a[k4];  in PetscKernel_A_gets_inverse_A_3()
63       a[k4]     = stmp;  in PetscKernel_A_gets_inverse_A_3()
67     stmp = -1. / a[k4];  in PetscKernel_A_gets_inverse_A_3()
69     aa   = &a[1 + k4];  in PetscKernel_A_gets_inverse_A_3()
73     ax = &a[k4 + 1];  in PetscKernel_A_gets_inverse_A_3()
76       stmp = a[l + j3];  in PetscKernel_A_gets_inverse_A_3()
78         a[l + j3] = a[k + j3];  in PetscKernel_A_gets_inverse_A_3()
79         a[k + j3] = stmp;  in PetscKernel_A_gets_inverse_A_3()
83       ay   = &a[1 + k + j3];  in PetscKernel_A_gets_inverse_A_3()
88   if (a[12] == 0.0) {  in PetscKernel_A_gets_inverse_A_3()
99     a[k4] = 1.0 / a[k4];  in PetscKernel_A_gets_inverse_A_3()
100     stmp  = -a[k4];  in PetscKernel_A_gets_inverse_A_3()
102     aa    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_3()
109       stmp      = a[k + j3];  in PetscKernel_A_gets_inverse_A_3()
110       a[k + j3] = 0.0;  in PetscKernel_A_gets_inverse_A_3()
111       ay        = &a[j3 + 1];  in PetscKernel_A_gets_inverse_A_3()
121     aa  = a + k3;  in PetscKernel_A_gets_inverse_A_3()
128       ax   = &a[3 * j + 1];  in PetscKernel_A_gets_inverse_A_3()
129       ay   = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_3()
136       ax    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_3()
137       ay    = &a[3 * l + 1];  in PetscKernel_A_gets_inverse_A_3()