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