| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvnat5.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, x1, x2, x3, x4, x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() local
39 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
42 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
43 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
44 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
70 x5 = x[4 + idx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
71 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
72 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvnat6.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() local
41 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
43 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
44 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
45 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
46 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
47 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
48 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
76 x5 = x[4 + idx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
78 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvnat7.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() local
43 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
46 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
47 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
48 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
49 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
50 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
51 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
52 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
82 x5 = x[4 + idx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolv.c | 71 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x, *t; in MatSolve_SeqBAIJ_7_inplace() local
112 x5 = t[4 + idx]; in MatSolve_SeqBAIJ_7_inplace()
115 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_inplace()
116 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_inplace()
117 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_inplace()
118 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_inplace()
119 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_inplace()
120 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_inplace()
121 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_inplace()
152 x5 = t[4 + idx]; in MatSolve_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtrannat5.c | 9 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace() local
24 x5 = x[4 + idx]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
25 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
26 s2 = v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
27 s3 = v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
28 s4 = v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
29 s5 = v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
83 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering() local
98 x5 = x[4 + idx]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
99 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat6.c | 9 PetscScalar s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6, *x; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace() local
24 x5 = x[4 + idx]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
26 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
27 s2 = v[6] * x1 + v[7] * x2 + v[8] * x3 + v[9] * x4 + v[10] * x5 + v[11] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
28 s3 = v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4 + v[16] * x5 + v[17] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
29 s4 = v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[22] * x5 + v[23] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
30 s5 = v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[29] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
31 s6 = v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
89 PetscScalar s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6, *x; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering() local
104 x5 = x[4 + idx]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat7.c | 9 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace() local
24 x5 = x[4 + idx]; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
27 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6 + v[6] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
28 s2 = v[7] * x1 + v[8] * x2 + v[9] * x3 + v[10] * x4 + v[11] * x5 + v[12] * x6 + v[13] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
29 s3 = v[14] * x1 + v[15] * x2 + v[16] * x3 + v[17] * x4 + v[18] * x5 + v[19] * x6 + v[20] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
30 s4 = v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[25] * x5 + v[26] * x6 + v[27] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
31 s5 = v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[33] * x6 + v[34] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
32 s6 = v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[41] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
33 s7 = v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48] * x7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
94 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering() local
[all …]
|
| H A D | baijsolvnat15.c | 13 PetscScalar x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2() local
65 x5 = x[4 + idx]; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
77 …s1 -= v[0] * x1 + v[15] * x2 + v[30] * x3 + v[45] * x4 + v[60] * x5 + v[75] * x6 + v[90] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
78 …s2 -= v[1] * x1 + v[16] * x2 + v[31] * x3 + v[46] * x4 + v[61] * x5 + v[76] * x6 + v[91] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
79 …s3 -= v[2] * x1 + v[17] * x2 + v[32] * x3 + v[47] * x4 + v[62] * x5 + v[77] * x6 + v[92] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
80 …s4 -= v[3] * x1 + v[18] * x2 + v[33] * x3 + v[48] * x4 + v[63] * x5 + v[78] * x6 + v[93] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
81 …s5 -= v[4] * x1 + v[19] * x2 + v[34] * x3 + v[49] * x4 + v[64] * x5 + v[79] * x6 + v[94] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
82 …s6 -= v[5] * x1 + v[20] * x2 + v[35] * x3 + v[50] * x4 + v[65] * x5 + v[80] * x6 + v[95] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
83 …s7 -= v[6] * x1 + v[21] * x2 + v[36] * x3 + v[51] * x4 + v[66] * x5 + v[81] * x6 + v[96] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
84 …s8 -= v[7] * x1 + v[22] * x2 + v[37] * x3 + v[52] * x4 + v[67] * x5 + v[82] * x6 + v[97] * x7 + v[… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
[all …]
|
| H A D | baijsolvtran5.c | 12 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x, *t; in MatSolveTranspose_SeqBAIJ_5_inplace() local
46 x5 = t[4 + idx]; in MatSolveTranspose_SeqBAIJ_5_inplace()
47 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
48 s2 = v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
49 s3 = v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
50 s4 = v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
51 s5 = v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
123 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x, *t; in MatSolveTranspose_SeqBAIJ_5() local
156 x5 = t[4 + idx]; in MatSolveTranspose_SeqBAIJ_5()
157 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5()
[all …]
|
| H A D | baijsolvtran7.c | 12 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x, *t; in MatSolveTranspose_SeqBAIJ_7_inplace() local
48 x5 = t[4 + idx]; in MatSolveTranspose_SeqBAIJ_7_inplace()
51 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6 + v[6] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
52 s2 = v[7] * x1 + v[8] * x2 + v[9] * x3 + v[10] * x4 + v[11] * x5 + v[12] * x6 + v[13] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
53 s3 = v[14] * x1 + v[15] * x2 + v[16] * x3 + v[17] * x4 + v[18] * x5 + v[19] * x6 + v[20] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
54 s4 = v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[25] * x5 + v[26] * x6 + v[27] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
55 s5 = v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[33] * x6 + v[34] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
56 s6 = v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[41] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
57 s7 = v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48] * x7; in MatSolveTranspose_SeqBAIJ_7_inplace()
138 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x, *t; in MatSolveTranspose_SeqBAIJ_7() local
[all …]
|
| H A D | baijsolvtran6.c | 12 PetscScalar s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6, *x, *t; in MatSolveTranspose_SeqBAIJ_6_inplace() local
47 x5 = t[4 + idx]; in MatSolveTranspose_SeqBAIJ_6_inplace()
49 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
50 s2 = v[6] * x1 + v[7] * x2 + v[8] * x3 + v[9] * x4 + v[10] * x5 + v[11] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
51 s3 = v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4 + v[16] * x5 + v[17] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
52 s4 = v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[22] * x5 + v[23] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
53 s5 = v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[29] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
54 s6 = v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatSolveTranspose_SeqBAIJ_6_inplace()
131 PetscScalar s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6, *x, *t; in MatSolveTranspose_SeqBAIJ_6() local
165 x5 = t[4 + idx]; in MatSolveTranspose_SeqBAIJ_6()
[all …]
|
| H A D | baij2.c | 459 PetscScalar *z = NULL, sum1, sum2, sum3, sum4, sum5, x1, x2, x3, x4, x5, *zarray; in MatMult_SeqBAIJ_5() local
499 x5 = xb[4]; in MatMult_SeqBAIJ_5()
500 sum1 += v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatMult_SeqBAIJ_5()
501 sum2 += v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatMult_SeqBAIJ_5()
502 sum3 += v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatMult_SeqBAIJ_5()
503 sum4 += v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatMult_SeqBAIJ_5()
504 sum5 += v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatMult_SeqBAIJ_5()
526 PetscScalar x1, x2, x3, x4, x5, x6, *zarray; in MatMult_SeqBAIJ_6() local
566 x5 = xb[4]; in MatMult_SeqBAIJ_6()
568 sum1 += v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatMult_SeqBAIJ_6()
[all …]
|
| H A D | baijfact11.c | 21 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local
98 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
115 pc[4] = m5 = p1 * x5 + p5 * x6 + p9 * x7 + p13 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
116 pc[5] = m6 = p2 * x5 + p6 * x6 + p10 * x7 + p14 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
117 pc[6] = m7 = p3 * x5 + p7 * x6 + p11 * x7 + p15 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
118 pc[7] = m8 = p4 * x5 + p8 * x6 + p12 * x7 + p16 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
137 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
155 x[4] -= m1 * x5 + m5 * x6 + m9 * x7 + m13 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
156 x[5] -= m2 * x5 + m6 * x6 + m10 * x7 + m14 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
157 x[6] -= m3 * x5 + m7 * x6 + m11 * x7 + m15 * x8; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
[all …]
|
| H A D | baijfact13.c | 21 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
80 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
89 pc[3] = m4 = p1 * x4 + p4 * x5 + p7 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
90 pc[4] = m5 = p2 * x4 + p5 * x5 + p8 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
91 pc[5] = m6 = p3 * x4 + p6 * x5 + p9 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
103 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
113 x[3] -= m1 * x4 + m4 * x5 + m7 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
114 x[4] -= m2 * x4 + m5 * x5 + m8 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
115 x[5] -= m3 * x4 + m6 * x5 + m9 * x6; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
281 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() local
[all …]
|
| H A D | baijfact7.c | 18 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16; in MatILUFactorNumeric_SeqBAIJ_6_inplace() local
141 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
173 pc[0] = m1 = p1 * x1 + p7 * x2 + p13 * x3 + p19 * x4 + p25 * x5 + p31 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
174 pc[1] = m2 = p2 * x1 + p8 * x2 + p14 * x3 + p20 * x4 + p26 * x5 + p32 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
175 pc[2] = m3 = p3 * x1 + p9 * x2 + p15 * x3 + p21 * x4 + p27 * x5 + p33 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
176 pc[3] = m4 = p4 * x1 + p10 * x2 + p16 * x3 + p22 * x4 + p28 * x5 + p34 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
177 pc[4] = m5 = p5 * x1 + p11 * x2 + p17 * x3 + p23 * x4 + p29 * x5 + p35 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
178 pc[5] = m6 = p6 * x1 + p12 * x2 + p18 * x3 + p24 * x4 + p30 * x5 + p36 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
222 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
255 x[0] -= m1 * x1 + m7 * x2 + m13 * x3 + m19 * x4 + m25 * x5 + m31 * x6; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
[all …]
|
| H A D | baijfact9.c | 22 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local
129 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
150 pc[0] = m1 = p1 * x1 + p6 * x2 + p11 * x3 + p16 * x4 + p21 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
151 pc[1] = m2 = p2 * x1 + p7 * x2 + p12 * x3 + p17 * x4 + p22 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
152 pc[2] = m3 = p3 * x1 + p8 * x2 + p13 * x3 + p18 * x4 + p23 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
153 pc[3] = m4 = p4 * x1 + p9 * x2 + p14 * x3 + p19 * x4 + p24 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
154 pc[4] = m5 = p5 * x1 + p10 * x2 + p15 * x3 + p20 * x4 + p25 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
187 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
209 x[0] -= m1 * x1 + m6 * x2 + m11 * x3 + m16 * x4 + m21 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
210 x[1] -= m2 * x1 + m7 * x2 + m12 * x3 + m17 * x4 + m22 * x5; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
[all …]
|
| H A D | baijfact5.c | 18 MatScalar p5, p6, p7, p8, p9, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16; in MatILUFactorNumeric_SeqBAIJ_7_inplace() local
171 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
216 pc[0] = m1 = p1 * x1 + p8 * x2 + p15 * x3 + p22 * x4 + p29 * x5 + p36 * x6 + p43 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
217 pc[1] = m2 = p2 * x1 + p9 * x2 + p16 * x3 + p23 * x4 + p30 * x5 + p37 * x6 + p44 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
218 pc[2] = m3 = p3 * x1 + p10 * x2 + p17 * x3 + p24 * x4 + p31 * x5 + p38 * x6 + p45 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
219 pc[3] = m4 = p4 * x1 + p11 * x2 + p18 * x3 + p25 * x4 + p32 * x5 + p39 * x6 + p46 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
220 pc[4] = m5 = p5 * x1 + p12 * x2 + p19 * x3 + p26 * x4 + p33 * x5 + p40 * x6 + p47 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
221 pc[5] = m6 = p6 * x1 + p13 * x2 + p20 * x3 + p27 * x4 + p34 * x5 + p41 * x6 + p48 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
222 pc[6] = m7 = p7 * x1 + p14 * x2 + p21 * x3 + p28 * x4 + p35 * x5 + p42 * x6 + p49 * x7; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
279 x5 = pv[4]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
[all …]
|
| /petsc/src/ksp/pc/impls/pbjacobi/ |
| H A D | pbjacobi.c | 11 PetscScalar *yy, x0, x1, x2, x3, x4, x5, x6; in PCApply_PBJacobi() local
79 x5 = xx[6 * i + 5]; in PCApply_PBJacobi()
81 … = diag[0] * x0 + diag[6] * x1 + diag[12] * x2 + diag[18] * x3 + diag[24] * x4 + diag[30] * x5; in PCApply_PBJacobi()
82 …+ 1] = diag[1] * x0 + diag[7] * x1 + diag[13] * x2 + diag[19] * x3 + diag[25] * x4 + diag[31] * x5; in PCApply_PBJacobi()
83 …+ 2] = diag[2] * x0 + diag[8] * x1 + diag[14] * x2 + diag[20] * x3 + diag[26] * x4 + diag[32] * x5; in PCApply_PBJacobi()
84 …+ 3] = diag[3] * x0 + diag[9] * x1 + diag[15] * x2 + diag[21] * x3 + diag[27] * x4 + diag[33] * x5; in PCApply_PBJacobi()
85 … 4] = diag[4] * x0 + diag[10] * x1 + diag[16] * x2 + diag[22] * x3 + diag[28] * x4 + diag[34] * x5; in PCApply_PBJacobi()
86 … 5] = diag[5] * x0 + diag[11] * x1 + diag[17] * x2 + diag[23] * x3 + diag[29] * x4 + diag[35] * x5; in PCApply_PBJacobi()
97 x5 = xx[7 * i + 5]; in PCApply_PBJacobi()
100 … x0 + diag[7] * x1 + diag[14] * x2 + diag[21] * x3 + diag[28] * x4 + diag[35] * x5 + diag[42] * x6; in PCApply_PBJacobi()
[all …]
|
| /petsc/src/ksp/pc/impls/vpbjacobi/ |
| H A D | vpbjacobi.c | 11 PetscScalar *yy, x0, x1, x2, x3, x4, x5, x6; in PCApply_VPBJacobi() local
67 x5 = xx[ncnt + 5]; in PCApply_VPBJacobi()
68 … = diag[0] * x0 + diag[6] * x1 + diag[12] * x2 + diag[18] * x3 + diag[24] * x4 + diag[30] * x5; in PCApply_VPBJacobi()
69 …+ 1] = diag[1] * x0 + diag[7] * x1 + diag[13] * x2 + diag[19] * x3 + diag[25] * x4 + diag[31] * x5; in PCApply_VPBJacobi()
70 …+ 2] = diag[2] * x0 + diag[8] * x1 + diag[14] * x2 + diag[20] * x3 + diag[26] * x4 + diag[32] * x5; in PCApply_VPBJacobi()
71 …+ 3] = diag[3] * x0 + diag[9] * x1 + diag[15] * x2 + diag[21] * x3 + diag[27] * x4 + diag[33] * x5; in PCApply_VPBJacobi()
72 … 4] = diag[4] * x0 + diag[10] * x1 + diag[16] * x2 + diag[22] * x3 + diag[28] * x4 + diag[34] * x5; in PCApply_VPBJacobi()
73 … 5] = diag[5] * x0 + diag[11] * x1 + diag[17] * x2 + diag[23] * x3 + diag[29] * x4 + diag[35] * x5; in PCApply_VPBJacobi()
81 x5 = xx[ncnt + 5]; in PCApply_VPBJacobi()
83 … x0 + diag[7] * x1 + diag[14] * x2 + diag[21] * x3 + diag[28] * x4 + diag[35] * x5 + diag[42] * x6; in PCApply_VPBJacobi()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 464 PetscScalar *z, x1, x2, x3, x4, x5, zero = 0.0; in MatMult_SeqSBAIJ_5() local
486 x5 = xb[4]; in MatMult_SeqSBAIJ_5()
491 z[5 * i] += v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatMult_SeqSBAIJ_5()
492 z[5 * i + 1] += v[5] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatMult_SeqSBAIJ_5()
493 z[5 * i + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatMult_SeqSBAIJ_5()
494 z[5 * i + 3] += v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[23] * x5; in MatMult_SeqSBAIJ_5()
495 z[5 * i + 4] += v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatMult_SeqSBAIJ_5()
504 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatMult_SeqSBAIJ_5()
505 z[cval + 1] += v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatMult_SeqSBAIJ_5()
506 z[cval + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatMult_SeqSBAIJ_5()
[all …]
|
| H A D | sbaijfact2.c | 222 PetscScalar *x, x0, x1, x2, x3, x4, x5, x6, *t, *tp; in MatSolve_SeqSBAIJ_7_inplace() local
254 x5 = tp[5]; in MatSolve_SeqSBAIJ_7_inplace()
259 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3 + v[4] * x4 + v[5] * x5 + v[6] * x6; in MatSolve_SeqSBAIJ_7_inplace()
260 … tp[1] += v[7] * x0 + v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4 + v[12] * x5 + v[13] * x6; in MatSolve_SeqSBAIJ_7_inplace()
261 … tp[2] += v[14] * x0 + v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5 + v[20] * x6; in MatSolve_SeqSBAIJ_7_inplace()
262 … tp[3] += v[21] * x0 + v[22] * x1 + v[23] * x2 + v[24] * x3 + v[25] * x4 + v[26] * x5 + v[27] * x6; in MatSolve_SeqSBAIJ_7_inplace()
263 … tp[4] += v[28] * x0 + v[29] * x1 + v[30] * x2 + v[31] * x3 + v[32] * x4 + v[33] * x5 + v[34] * x6; in MatSolve_SeqSBAIJ_7_inplace()
264 … tp[5] += v[35] * x0 + v[36] * x1 + v[37] * x2 + v[38] * x3 + v[39] * x4 + v[40] * x5 + v[41] * x6; in MatSolve_SeqSBAIJ_7_inplace()
265 … tp[6] += v[42] * x0 + v[43] * x1 + v[44] * x2 + v[45] * x3 + v[46] * x4 + v[47] * x5 + v[48] * x6; in MatSolve_SeqSBAIJ_7_inplace()
274 tp[0] = d[0] * x0 + d[7] * x1 + d[14] * x2 + d[21] * x3 + d[28] * x4 + d[35] * x5 + d[42] * x6; in MatSolve_SeqSBAIJ_7_inplace()
[all …]
|
| /petsc/src/vec/pf/impls/string/ |
| H A D | cstringbase.template | 8 PetscScalar x,y,z,f = NOTLIKELY,x1 = 0,x2 = 0,x3 = 0,x4 = 0,x5 = 0;
19 (void)x5;
37 x5 = in[_NIN_*i+4];
|
| /petsc/src/snes/utils/dm/ |
| H A D | dminterpolatesnes.c | 641 const PetscScalar x5 = vertices[15]; in HexMap_Private() local
665 const PetscScalar f_02 = x5 - x1 - x4 + x0; in HexMap_Private()
668 const PetscScalar f_012 = x6 - x0 + x1 - x2 + x3 + x4 - x5 - x7; in HexMap_Private()
710 const PetscScalar x5 = vertices[15]; in HexJacobian_Private() local
725 const PetscScalar f_xz = x5 - x1 - x4 + x0; in HexJacobian_Private()
728 const PetscScalar f_xyz = x6 - x0 + x1 - x2 + x3 + x4 - x5 - x7; in HexJacobian_Private()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | inode.c | 2498 PetscScalar *x, tmp4, tmp5, x1, x2, x3, x4, x5; in MatSOR_SeqAIJ_Inode() local
3607 x5 = x[row + 4]; in MatSOR_SeqAIJ_Inode()
3608 …ag[cnt] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 15] + x5 * bdiag[cnt + 20]; in MatSOR_SeqAIJ_Inode()
3609 …nt + 1] + x2 * bdiag[cnt + 6] + x3 * bdiag[cnt + 11] + x4 * bdiag[cnt + 16] + x5 * bdiag[cnt + 21]; in MatSOR_SeqAIJ_Inode()
3610 …nt + 2] + x2 * bdiag[cnt + 7] + x3 * bdiag[cnt + 12] + x4 * bdiag[cnt + 17] + x5 * bdiag[cnt + 22]; in MatSOR_SeqAIJ_Inode()
3611 …nt + 3] + x2 * bdiag[cnt + 8] + x3 * bdiag[cnt + 13] + x4 * bdiag[cnt + 18] + x5 * bdiag[cnt + 23]; in MatSOR_SeqAIJ_Inode()
3612 …nt + 4] + x2 * bdiag[cnt + 9] + x3 * bdiag[cnt + 14] + x4 * bdiag[cnt + 19] + x5 * bdiag[cnt + 24]; in MatSOR_SeqAIJ_Inode()
3811 PetscScalar *x, tmp1, tmp2, tmp3, tmp4, tmp5, x1, x2, x3, x4, x5; in MatMultDiagonalBlock_SeqAIJ_Inode() local
3870 x5 = b[row + 4]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3871 …ag[cnt] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 15] + x5 * bdiag[cnt + 20]; in MatMultDiagonalBlock_SeqAIJ_Inode()
[all …]
|