14 PetscErrorCode PetscKernel_A_gets_inverse_A_2(MatScalar *a, PetscReal shift, PetscBool allowzeropiv…  in PetscKernel_A_gets_inverse_A_2()  argument
23   shift = .25 * shift * (1.e-12 + PetscAbsScalar(a[0]) + PetscAbsScalar(a[3]));  in PetscKernel_A_gets_inverse_A_2()
26   a -= 3;  in PetscKernel_A_gets_inverse_A_2()
35   aa   = &a[k4];  in PetscKernel_A_gets_inverse_A_2()
48   if (a[l + k3] == 0.0) {  in PetscKernel_A_gets_inverse_A_2()
54       a[l + k3] = shift;  in PetscKernel_A_gets_inverse_A_2()
60     stmp      = a[l + k3];  in PetscKernel_A_gets_inverse_A_2()
61     a[l + k3] = a[k4];  in PetscKernel_A_gets_inverse_A_2()
62     a[k4]     = stmp;  in PetscKernel_A_gets_inverse_A_2()
66   stmp = -1. / a[k4];  in PetscKernel_A_gets_inverse_A_2()
68   aa   = &a[1 + k4];  in PetscKernel_A_gets_inverse_A_2()
72   ax = &a[k4 + 1];  in PetscKernel_A_gets_inverse_A_2()
75     stmp = a[l + j3];  in PetscKernel_A_gets_inverse_A_2()
77       a[l + j3] = a[k + j3];  in PetscKernel_A_gets_inverse_A_2()
78       a[k + j3] = stmp;  in PetscKernel_A_gets_inverse_A_2()
82     ay   = &a[1 + k + j3];  in PetscKernel_A_gets_inverse_A_2()
87   if (a[6] == 0.0) {  in PetscKernel_A_gets_inverse_A_2()
98     a[k4] = 1.0 / a[k4];  in PetscKernel_A_gets_inverse_A_2()
99     stmp  = -a[k4];  in PetscKernel_A_gets_inverse_A_2()
101     aa    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_2()
108       stmp      = a[k + j3];  in PetscKernel_A_gets_inverse_A_2()
109       a[k + j3] = 0.0;  in PetscKernel_A_gets_inverse_A_2()
110       ay        = &a[j3 + 1];  in PetscKernel_A_gets_inverse_A_2()
119   aa  = a + k3;  in PetscKernel_A_gets_inverse_A_2()
126     ax   = &a[2 * j + 1];  in PetscKernel_A_gets_inverse_A_2()
127     ay   = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_2()
133     ax    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_2()
134     ay    = &a[2 * l + 1];  in PetscKernel_A_gets_inverse_A_2()
146 PetscErrorCode PetscKernel_A_gets_inverse_A_9(MatScalar *a, PetscReal shift, PetscBool allowzeropiv…  in PetscKernel_A_gets_inverse_A_9()  argument
157   a -= 10;  in PetscKernel_A_gets_inverse_A_9()
166     aa   = &a[k4];  in PetscKernel_A_gets_inverse_A_9()
179     if (a[l + k3] == 0.0) {  in PetscKernel_A_gets_inverse_A_9()
185         a[l + k3] = shift;  in PetscKernel_A_gets_inverse_A_9()
191       stmp      = a[l + k3];  in PetscKernel_A_gets_inverse_A_9()
192       a[l + k3] = a[k4];  in PetscKernel_A_gets_inverse_A_9()
193       a[k4]     = stmp;  in PetscKernel_A_gets_inverse_A_9()
197     stmp = -1. / a[k4];  in PetscKernel_A_gets_inverse_A_9()
199     aa   = &a[1 + k4];  in PetscKernel_A_gets_inverse_A_9()
203     ax = &a[k4 + 1];  in PetscKernel_A_gets_inverse_A_9()
206       stmp = a[l + j3];  in PetscKernel_A_gets_inverse_A_9()
208         a[l + j3] = a[k + j3];  in PetscKernel_A_gets_inverse_A_9()
209         a[k + j3] = stmp;  in PetscKernel_A_gets_inverse_A_9()
213       ay   = &a[1 + k + j3];  in PetscKernel_A_gets_inverse_A_9()
218   if (a[90] == 0.0) {  in PetscKernel_A_gets_inverse_A_9()
229     a[k4] = 1.0 / a[k4];  in PetscKernel_A_gets_inverse_A_9()
230     stmp  = -a[k4];  in PetscKernel_A_gets_inverse_A_9()
232     aa    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_9()
239       stmp      = a[k + j3];  in PetscKernel_A_gets_inverse_A_9()
240       a[k + j3] = 0.0;  in PetscKernel_A_gets_inverse_A_9()
241       ay        = &a[j3 + 1];  in PetscKernel_A_gets_inverse_A_9()
251     aa  = a + k3;  in PetscKernel_A_gets_inverse_A_9()
258       ax   = &a[9 * j + 1];  in PetscKernel_A_gets_inverse_A_9()
259       ay   = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_9()
272       ax    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_9()
273       ay    = &a[9 * l + 1];  in PetscKernel_A_gets_inverse_A_9()
317 PetscErrorCode PetscKernel_A_gets_inverse_A_15(MatScalar *a, PetscInt *ipvt, MatScalar *work, Petsc…  in PetscKernel_A_gets_inverse_A_15()  argument
328   a -= 16;  in PetscKernel_A_gets_inverse_A_15()
337     aa   = &a[k4];  in PetscKernel_A_gets_inverse_A_15()
350     if (a[l + k3] == 0.0) {  in PetscKernel_A_gets_inverse_A_15()
356         a[l + k3] = shift;  in PetscKernel_A_gets_inverse_A_15()
362       stmp      = a[l + k3];  in PetscKernel_A_gets_inverse_A_15()
363       a[l + k3] = a[k4];  in PetscKernel_A_gets_inverse_A_15()
364       a[k4]     = stmp;  in PetscKernel_A_gets_inverse_A_15()
368     stmp = -1. / a[k4];  in PetscKernel_A_gets_inverse_A_15()
370     aa   = &a[1 + k4];  in PetscKernel_A_gets_inverse_A_15()
374     ax = &a[k4 + 1];  in PetscKernel_A_gets_inverse_A_15()
377       stmp = a[l + j3];  in PetscKernel_A_gets_inverse_A_15()
379         a[l + j3] = a[k + j3];  in PetscKernel_A_gets_inverse_A_15()
380         a[k + j3] = stmp;  in PetscKernel_A_gets_inverse_A_15()
384       ay   = &a[1 + k + j3];  in PetscKernel_A_gets_inverse_A_15()
389   if (a[240] == 0.0) {  in PetscKernel_A_gets_inverse_A_15()
400     a[k4] = 1.0 / a[k4];  in PetscKernel_A_gets_inverse_A_15()
401     stmp  = -a[k4];  in PetscKernel_A_gets_inverse_A_15()
403     aa    = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_15()
410       stmp      = a[k + j3];  in PetscKernel_A_gets_inverse_A_15()
411       a[k + j3] = 0.0;  in PetscKernel_A_gets_inverse_A_15()
412       ay        = &a[j3 + 1];  in PetscKernel_A_gets_inverse_A_15()
422     aa  = a + k3;  in PetscKernel_A_gets_inverse_A_15()
429       ax   = &a[15 * j + 1];  in PetscKernel_A_gets_inverse_A_15()
430       ay   = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_15()
449       ax     = &a[k3 + 1];  in PetscKernel_A_gets_inverse_A_15()
450       ay     = &a[15 * l + 1];  in PetscKernel_A_gets_inverse_A_15()