| /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
54 x4 = xx[ncnt + 4]; in PCApply_VPBJacobi()
55 yy[ncnt] = diag[0] * x0 + diag[5] * x1 + diag[10] * x2 + diag[15] * x3 + diag[20] * x4; in PCApply_VPBJacobi()
56 yy[ncnt + 1] = diag[1] * x0 + diag[6] * x1 + diag[11] * x2 + diag[16] * x3 + diag[21] * x4; in PCApply_VPBJacobi()
57 yy[ncnt + 2] = diag[2] * x0 + diag[7] * x1 + diag[12] * x2 + diag[17] * x3 + diag[22] * x4; in PCApply_VPBJacobi()
58 yy[ncnt + 3] = diag[3] * x0 + diag[8] * x1 + diag[13] * x2 + diag[18] * x3 + diag[23] * x4; in PCApply_VPBJacobi()
59 yy[ncnt + 4] = diag[4] * x0 + diag[9] * x1 + diag[14] * x2 + diag[19] * x3 + diag[24] * x4; in PCApply_VPBJacobi()
66 x4 = xx[ncnt + 4]; in PCApply_VPBJacobi()
68 …yy[ncnt] = diag[0] * x0 + diag[6] * x1 + diag[12] * x2 + diag[18] * x3 + diag[24] * x4 + diag[… in PCApply_VPBJacobi()
69 …yy[ncnt + 1] = diag[1] * x0 + diag[7] * x1 + diag[13] * x2 + diag[19] * x3 + diag[25] * x4 + diag[… in PCApply_VPBJacobi()
[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
62 x4 = xx[5 * i + 4]; in PCApply_PBJacobi()
64 yy[5 * i] = diag[0] * x0 + diag[5] * x1 + diag[10] * x2 + diag[15] * x3 + diag[20] * x4; in PCApply_PBJacobi()
65 yy[5 * i + 1] = diag[1] * x0 + diag[6] * x1 + diag[11] * x2 + diag[16] * x3 + diag[21] * x4; in PCApply_PBJacobi()
66 yy[5 * i + 2] = diag[2] * x0 + diag[7] * x1 + diag[12] * x2 + diag[17] * x3 + diag[22] * x4; in PCApply_PBJacobi()
67 yy[5 * i + 3] = diag[3] * x0 + diag[8] * x1 + diag[13] * x2 + diag[18] * x3 + diag[23] * x4; in PCApply_PBJacobi()
68 yy[5 * i + 4] = diag[4] * x0 + diag[9] * x1 + diag[14] * x2 + diag[19] * x3 + diag[24] * x4; in PCApply_PBJacobi()
78 x4 = xx[6 * i + 4]; in PCApply_PBJacobi()
81 …yy[6 * i] = diag[0] * x0 + diag[6] * x1 + diag[12] * x2 + diag[18] * x3 + diag[24] * x4 + diag… in PCApply_PBJacobi()
82 …yy[6 * i + 1] = diag[1] * x0 + diag[7] * x1 + diag[13] * x2 + diag[19] * x3 + diag[25] * x4 + diag… in PCApply_PBJacobi()
[all …]
|
| /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
38 x4 = x[3 + 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()
69 x4 = x[3 + 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
40 x4 = x[3 + 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()
75 x4 = x[3 + 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 | baijsolvnat4.c | 31 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
56 x4 = x[3 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
57 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
58 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
59 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
60 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
84 x4 = x[3 + idx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
85 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
86 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
87 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_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
42 x4 = x[3 + 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()
81 x4 = x[3 + idx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
|
| H A D | baij2.c | 398 PetscScalar *z = NULL, sum1, sum2, sum3, sum4, x1, x2, x3, x4, *zarray; in MatMult_SeqBAIJ_4() local
436 x4 = xb[3]; in MatMult_SeqBAIJ_4()
437 sum1 += v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatMult_SeqBAIJ_4()
438 sum2 += v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatMult_SeqBAIJ_4()
439 sum3 += v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatMult_SeqBAIJ_4()
440 sum4 += v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatMult_SeqBAIJ_4()
459 PetscScalar *z = NULL, sum1, sum2, sum3, sum4, sum5, x1, x2, x3, x4, x5, *zarray; in MatMult_SeqBAIJ_5() local
498 x4 = xb[3]; 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()
[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
111 x4 = t[3 + 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()
151 x4 = t[3 + idx]; in MatSolve_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtrannat4.c | 9 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace() local
23 x4 = x[3 + idx]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
24 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
25 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
26 s3 = v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
27 s4 = v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
77 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering() local
91 x4 = x[3 + idx]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
92 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
93 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
[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
23 x4 = x[3 + 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
97 x4 = x[3 + 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
23 x4 = x[3 + 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
103 x4 = x[3 + 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
23 x4 = x[3 + 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
64 x4 = x[3 + 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
45 x4 = t[3 + 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
155 x4 = t[3 + 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
47 x4 = t[3 + 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 | baijsolvtran4.c | 12 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x, *t; in MatSolveTranspose_SeqBAIJ_4_inplace() local
44 x4 = t[3 + idx]; in MatSolveTranspose_SeqBAIJ_4_inplace()
45 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4_inplace()
46 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_inplace()
47 s3 = v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatSolveTranspose_SeqBAIJ_4_inplace()
48 s4 = v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatSolveTranspose_SeqBAIJ_4_inplace()
115 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x, *t; in MatSolveTranspose_SeqBAIJ_4() local
146 x4 = t[3 + idx]; in MatSolveTranspose_SeqBAIJ_4()
147 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4()
148 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4()
[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
46 x4 = t[3 + 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
164 x4 = t[3 + idx]; in MatSolveTranspose_SeqBAIJ_6()
[all …]
|
| H A D | baijfact11.c | 20 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local
97 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
110 pc[0] = m1 = p1 * x1 + p5 * x2 + p9 * x3 + p13 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
111 pc[1] = m2 = p2 * x1 + p6 * x2 + p10 * x3 + p14 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
112 pc[2] = m3 = p3 * x1 + p7 * x2 + p11 * x3 + p15 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
113 pc[3] = m4 = p4 * x1 + p8 * x2 + p12 * x3 + p16 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
136 x4 = pv[3]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
150 x[0] -= m1 * x1 + m5 * x2 + m9 * x3 + m13 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
151 x[1] -= m2 * x1 + m6 * x2 + m10 * x3 + m14 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
152 x[2] -= m3 * x1 + m7 * x2 + m11 * x3 + m15 * x4; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
[all …]
|
| H A D | baijfact13.c | 20 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
79 x4 = pv[3]; 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()
102 x4 = pv[3]; 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()
280 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_3_NaturalOrdering_inplace() local
[all …]
|
| H A D | baijfact7.c | 17 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_6_inplace() local
140 x4 = pv[3]; 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()
221 x4 = pv[3]; 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 | 21 MatScalar p1, p2, p3, p4, m1, m2, m3, m4, m5, m6, m7, m8, m9, x1, x2, x3, x4; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local
128 x4 = pv[3]; 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()
186 x4 = pv[3]; 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 …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 402 PetscScalar *z, x1, x2, x3, x4, zero = 0.0; in MatMult_SeqSBAIJ_4() local
423 x4 = xb[3]; in MatMult_SeqSBAIJ_4()
428 z[4 * i] += v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatMult_SeqSBAIJ_4()
429 z[4 * i + 1] += v[4] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatMult_SeqSBAIJ_4()
430 z[4 * i + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[14] * x4; in MatMult_SeqSBAIJ_4()
431 z[4 * i + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMult_SeqSBAIJ_4()
440 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatMult_SeqSBAIJ_4()
441 z[cval + 1] += v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatMult_SeqSBAIJ_4()
442 z[cval + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatMult_SeqSBAIJ_4()
443 z[cval + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMult_SeqSBAIJ_4()
[all …]
|
| H A D | sbaijfact2.c | 222 PetscScalar *x, x0, x1, x2, x3, x4, x5, x6, *t, *tp; in MatSolve_SeqSBAIJ_7_inplace() local
253 x4 = tp[4]; 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;
18 (void)x4;
34 x4 = in[_NIN_*i+3];
|
| /petsc/src/dm/tutorials/ |
| H A D | ex11f90.F90 | 11 PetscScalar, pointer :: x3(:, :, :), x4(:, :, :, :)
146 PetscCallA(DMDAVecGetArray(ada, g, x4, ierr))
151 x4(0, i, j, k) = i + j + k
152 x4(1, i, j, k) = -(i + j + k)
153 x4(2, i, j, k) = i + j + k
158 PetscCallA(DMDAVecRestoreArray(ada, g, x4, ierr))
|