| /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
43 x3 = xx[ncnt + 3]; in PCApply_VPBJacobi()
44 yy[ncnt] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_VPBJacobi()
45 yy[ncnt + 1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3; in PCApply_VPBJacobi()
46 yy[ncnt + 2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3; in PCApply_VPBJacobi()
47 yy[ncnt + 3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3; in PCApply_VPBJacobi()
53 x3 = xx[ncnt + 3]; 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()
[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
47 x3 = xx[4 * i + 3]; in PCApply_PBJacobi()
49 yy[4 * i] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_PBJacobi()
50 yy[4 * i + 1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3; in PCApply_PBJacobi()
51 yy[4 * i + 2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3; in PCApply_PBJacobi()
52 yy[4 * i + 3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3; in PCApply_PBJacobi()
61 x3 = xx[5 * i + 3]; 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()
[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
37 x3 = x[2 + 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()
68 x3 = x[2 + 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
39 x3 = x[2 + 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()
74 x3 = x[2 + 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 | baijsolvnat3.c | 14 PetscScalar *x, s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
39 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
42 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
62 x3 = x[2 + idx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
63 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
64 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
65 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
89 PetscScalar s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering() local
[all …]
|
| H A D | baijsolvnat4.c | 31 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
55 x3 = x[2 + 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()
83 x3 = x[2 + 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 | baij2.c | 337 PetscScalar *z = NULL, sum1, sum2, sum3, x1, x2, x3, *zarray; in MatMult_SeqBAIJ_3() local
376 x3 = xb[2]; in MatMult_SeqBAIJ_3()
378 sum1 += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMult_SeqBAIJ_3()
379 sum2 += v[1] * x1 + v[4] * x2 + v[7] * x3; in MatMult_SeqBAIJ_3()
380 sum3 += v[2] * x1 + v[5] * x2 + v[8] * x3; in MatMult_SeqBAIJ_3()
398 PetscScalar *z = NULL, sum1, sum2, sum3, sum4, x1, x2, x3, x4, *zarray; in MatMult_SeqBAIJ_4() local
435 x3 = xb[2]; 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()
[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
41 x3 = x[2 + 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()
80 x3 = x[2 + 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
110 x3 = t[2 + 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()
150 x3 = t[2 + 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
22 x3 = x[2 + 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
90 x3 = x[2 + 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
22 x3 = x[2 + 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
96 x3 = x[2 + 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 | baijsolvtrannat3.c | 9 PetscScalar s1, s2, s3, x1, x2, x3, *x; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
22 x3 = x[2 + idx]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
23 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
24 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
25 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
71 PetscScalar s1, s2, s3, x1, x2, x3, *x; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering() local
84 x3 = x[2 + idx]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
85 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
86 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
87 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
|
| 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
22 x3 = x[2 + 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
102 x3 = x[2 + 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
22 x3 = x[2 + 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
63 x3 = x[2 + 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
44 x3 = t[2 + 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
154 x3 = t[2 + 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
46 x3 = t[2 + 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
43 x3 = t[2 + 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
145 x3 = t[2 + 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
45 x3 = t[2 + 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
163 x3 = t[2 + idx]; in MatSolveTranspose_SeqBAIJ_6()
[all …]
|
| H A D | baijsolvtran3.c | 12 PetscScalar s1, s2, s3, x1, x2, x3, *x, *t; in MatSolveTranspose_SeqBAIJ_3_inplace() local
42 x3 = t[2 + idx]; in MatSolveTranspose_SeqBAIJ_3_inplace()
43 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3_inplace()
44 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_inplace()
45 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3_inplace()
107 PetscScalar s1, s2, s3, x1, x2, x3, *x, *t; in MatSolveTranspose_SeqBAIJ_3() local
136 x3 = t[2 + idx]; in MatSolveTranspose_SeqBAIJ_3()
137 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3()
138 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3()
139 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3()
|
| /petsc/src/vec/vec/tests/ |
| H A D | ex60.c | 9 Vec x, x1, x2, x3; in main() local
29 PetscCall(VecCreateSeqCUDAWithArray(PETSC_COMM_SELF, bs, n, NULL, &x3)); in main()
35 PetscCall(VecCreateSeqKokkosWithArray(PETSC_COMM_SELF, bs, n, NULL, &x3)); in main()
41 PetscCall(VecCreateSeqHIPWithArray(PETSC_COMM_SELF, bs, n, NULL, &x3)); in main()
46 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, bs, n, NULL, &x3)); in main()
52 PetscCall(VecPlaceArray(x3, px + 2 * n)); in main()
55 PetscCall(VecSet(x3, 2.0)); in main()
58 PetscCall(VecResetArray(x3)); in main()
62 PetscCall(VecDestroy(&x3)); in main()
|
| /petsc/src/sys/tests/ |
| H A D | ex35.c | 11 PetscReal x3[] = {39, 9, 19, -39, 29, 309, 209, 390, 12, 11}; in main() local
36 … PetscCall(PetscPrintf(PETSC_COMM_SELF, " %2" PetscInt_FMT " %g\n", index3[i], (double)x3[i])); in main()
38 PetscCall(PetscSortRealWithArrayInt(10, x3, index3)); in main()
39 … PetscCall(PetscPrintf(PETSC_COMM_SELF, " %2" PetscInt_FMT " %g\n", index3[i], (double)x3[i])); in main()
43 PetscCall(PetscFindReal(val, 10, x3, PETSC_SMALL, &loc)); in main()
46 PetscCall(PetscFindReal(val, 10, x3, PETSC_SMALL, &loc)); in main()
49 PetscCall(PetscFindReal(val, 10, x3, 0.21, &loc)); in main()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 344 PetscScalar *z, x1, x2, x3, zero = 0.0; in MatMult_SeqSBAIJ_3() local
364 x3 = xb[2]; in MatMult_SeqSBAIJ_3()
369 z[3 * i] += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMult_SeqSBAIJ_3()
370 z[3 * i + 1] += v[3] * x1 + v[4] * x2 + v[7] * x3; in MatMult_SeqSBAIJ_3()
371 z[3 * i + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMult_SeqSBAIJ_3()
380 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3; in MatMult_SeqSBAIJ_3()
381 z[cval + 1] += v[3] * x1 + v[4] * x2 + v[5] * x3; in MatMult_SeqSBAIJ_3()
382 z[cval + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMult_SeqSBAIJ_3()
402 PetscScalar *z, x1, x2, x3, x4, zero = 0.0; in MatMult_SeqSBAIJ_4() local
422 x3 = xb[2]; 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
252 x3 = tp[3]; 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/sys/classes/draw/impls/x/ |
| H A D | xtone.c | 10 …iangle_X(PetscDraw_X *win, int x1, int y_1, int t1, int x2, int y2, int t2, int x3, int y3, int t3) in PetscDrawInterpolatedTriangle_X() argument
21 if (x1 < 0 && x2 < 0 && x3 < 0) PetscFunctionReturn(PETSC_SUCCESS); in PetscDrawInterpolatedTriangle_X()
23 if (x1 > win->w && x2 > win->w && x3 > win->w) PetscFunctionReturn(PETSC_SUCCESS); in PetscDrawInterpolatedTriangle_X()
46 SWAP(x1, x3); in PetscDrawInterpolatedTriangle_X()
51 SWAP(x2, x3); in PetscDrawInterpolatedTriangle_X()
63 x3_x1 = x3 - x1; in PetscDrawInterpolatedTriangle_X()
105 x3_x1 = x3 - x1; in PetscDrawInterpolatedTriangle_X()
108 x3_x2 = x3 - x2; in PetscDrawInterpolatedTriangle_X()
|