| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvtrannat4.c | 9 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace() local
25 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
34 x[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
35 x[oidx + 1] -= v[4] * s1 + v[5] * s2 + v[6] * s3 + v[7] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
36 x[oidx + 2] -= v[8] * s1 + v[9] * s2 + v[10] * s3 + v[11] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
37 x[oidx + 3] -= v[12] * s1 + v[13] * s2 + v[14] * s3 + v[15] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
41 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
53 s2 = x[1 + idt]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
58 x[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
59 x[idx + 1] -= v[4] * s1 + v[5] * s2 + v[6] * s3 + v[7] * s4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_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
26 s2 = v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
36 x[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
37 x[oidx + 1] -= v[5] * s1 + v[6] * s2 + v[7] * s3 + v[8] * s4 + v[9] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
38 x[oidx + 2] -= v[10] * s1 + v[11] * s2 + v[12] * s3 + v[13] * s4 + v[14] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
39 x[oidx + 3] -= v[15] * s1 + v[16] * s2 + v[17] * s3 + v[18] * s4 + v[19] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
40 x[oidx + 4] -= v[20] * s1 + v[21] * s2 + v[22] * s3 + v[23] * s4 + v[24] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
44 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
57 s2 = x[1 + idt]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
63 x[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvtrannat3.c | 9 PetscScalar s1, s2, s3, x1, x2, x3, *x; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
24 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
32 x[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
33 x[oidx + 1] -= v[3] * s1 + v[4] * s2 + v[5] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
34 x[oidx + 2] -= v[6] * s1 + v[7] * s2 + v[8] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
38 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
49 s2 = x[1 + idt]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
53 x[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
54 x[idx + 1] -= v[3] * s1 + v[4] * s2 + v[5] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
55 x[idx + 2] -= v[6] * s1 + v[7] * s2 + v[8] * s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
[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
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()
38 x[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5 + v[5] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
39 x[oidx + 1] -= v[6] * s1 + v[7] * s2 + v[8] * s3 + v[9] * s4 + v[10] * s5 + v[11] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
40 x[oidx + 2] -= v[12] * s1 + v[13] * s2 + v[14] * s3 + v[15] * s4 + v[16] * s5 + v[17] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
41 x[oidx + 3] -= v[18] * s1 + v[19] * s2 + v[20] * s3 + v[21] * s4 + v[22] * s5 + v[23] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
42 x[oidx + 4] -= v[24] * s1 + v[25] * s2 + v[26] * s3 + v[27] * s4 + v[28] * s5 + v[29] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
43 x[oidx + 5] -= v[30] * s1 + v[31] * s2 + v[32] * s3 + v[33] * s4 + v[34] * s5 + v[35] * s6; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
47 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
61 s2 = x[1 + idt]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvtrannat2.c | 9 PetscScalar s1, s2, x1, x2, *x; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace() local
23 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
30 x[oidx] -= v[0] * s1 + v[1] * s2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
31 x[oidx + 1] -= v[2] * s1 + v[3] * s2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
35 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
45 s2 = x[1 + idt]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
48 x[idx] -= v[0] * s1 + v[1] * s2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
49 x[idx + 1] -= v[2] * s1 + v[3] * s2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
65 PetscScalar s1, s2, x1, x2, *x; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering() local
79 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
[all …]
|
| H A D | baijsolvnat2.c | 17 PetscScalar *x, s1, s2, x1, x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
35 s2 = b[1 + idx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
45 x[1 + idx] = s2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
54 s2 = x[1 + idt]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
60 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
64 x[idt] = v[0] * s1 + v[2] * s2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
65 x[1 + idt] = v[1] * s1 + v[3] * s2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
80 PetscScalar *x, s1, s2, x1, x2; in MatSolve_SeqBAIJ_2_NaturalOrdering() local
96 s2 = b[1 + idx]; in MatSolve_SeqBAIJ_2_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
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()
40 x[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5 + v[5] * s6 + v[6] * s7; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
41 …x[oidx + 1] -= v[7] * s1 + v[8] * s2 + v[9] * s3 + v[10] * s4 + v[11] * s5 + v[12] * s6 + v[13] * … in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
42 …x[oidx + 2] -= v[14] * s1 + v[15] * s2 + v[16] * s3 + v[17] * s4 + v[18] * s5 + v[19] * s6 + v[20]… in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
43 …x[oidx + 3] -= v[21] * s1 + v[22] * s2 + v[23] * s3 + v[24] * s4 + v[25] * s5 + v[26] * s6 + v[27]… in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
44 …x[oidx + 4] -= v[28] * s1 + v[29] * s2 + v[30] * s3 + v[31] * s4 + v[32] * s5 + v[33] * s6 + v[34]… in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
45 …x[oidx + 5] -= v[35] * s1 + v[36] * s2 + v[37] * s3 + v[38] * s4 + v[39] * s5 + v[40] * s6 + v[41]… in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
46 …x[oidx + 6] -= v[42] * s1 + v[43] * s2 + v[44] * s3 + v[45] * s4 + v[46] * s5 + v[47] * s6 + v[48]… in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
50 x[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
[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
48 s2 = v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatSolveTranspose_SeqBAIJ_5_inplace()
58 t[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
59 t[oidx + 1] -= v[5] * s1 + v[6] * s2 + v[7] * s3 + v[8] * s4 + v[9] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
60 t[oidx + 2] -= v[10] * s1 + v[11] * s2 + v[12] * s3 + v[13] * s4 + v[14] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
61 t[oidx + 3] -= v[15] * s1 + v[16] * s2 + v[17] * s3 + v[18] * s4 + v[19] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
62 t[oidx + 4] -= v[20] * s1 + v[21] * s2 + v[22] * s3 + v[23] * s4 + v[24] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
66 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_5_inplace()
79 s2 = t[1 + idt]; in MatSolveTranspose_SeqBAIJ_5_inplace()
85 t[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5; in MatSolveTranspose_SeqBAIJ_5_inplace()
[all …]
|
| H A D | baijsolvnat3.c | 14 PetscScalar *x, s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
33 s2 = b[1 + idx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
46 x[1 + idx] = s2; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
56 s2 = x[1 + idt]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
64 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
69 x[idt] = v[0] * s1 + v[3] * s2 + v[6] * s3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
70 x[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
71 x[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
89 PetscScalar s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering() local
[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
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()
64 t[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5 + v[5] * s6 + v[6] * s7; in MatSolveTranspose_SeqBAIJ_7_inplace()
65 …t[oidx + 1] -= v[7] * s1 + v[8] * s2 + v[9] * s3 + v[10] * s4 + v[11] * s5 + v[12] * s6 + v[13] * … in MatSolveTranspose_SeqBAIJ_7_inplace()
66 …t[oidx + 2] -= v[14] * s1 + v[15] * s2 + v[16] * s3 + v[17] * s4 + v[18] * s5 + v[19] * s6 + v[20]… in MatSolveTranspose_SeqBAIJ_7_inplace()
67 …t[oidx + 3] -= v[21] * s1 + v[22] * s2 + v[23] * s3 + v[24] * s4 + v[25] * s5 + v[26] * s6 + v[27]… in MatSolveTranspose_SeqBAIJ_7_inplace()
68 …t[oidx + 4] -= v[28] * s1 + v[29] * s2 + v[30] * s3 + v[31] * s4 + v[32] * s5 + v[33] * s6 + v[34]… in MatSolveTranspose_SeqBAIJ_7_inplace()
69 …t[oidx + 5] -= v[35] * s1 + v[36] * s2 + v[37] * s3 + v[38] * s4 + v[39] * s5 + v[40] * s6 + v[41]… in MatSolveTranspose_SeqBAIJ_7_inplace()
70 …t[oidx + 6] -= v[42] * s1 + v[43] * s2 + v[44] * s3 + v[45] * s4 + v[46] * s5 + v[47] * s6 + v[48]… in MatSolveTranspose_SeqBAIJ_7_inplace()
74 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtran4.c | 12 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x, *t; in MatSolveTranspose_SeqBAIJ_4_inplace() local
46 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_inplace()
55 t[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
56 t[oidx + 1] -= v[4] * s1 + v[5] * s2 + v[6] * s3 + v[7] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
57 t[oidx + 2] -= v[8] * s1 + v[9] * s2 + v[10] * s3 + v[11] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
58 t[oidx + 3] -= v[12] * s1 + v[13] * s2 + v[14] * s3 + v[15] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
62 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_4_inplace()
74 s2 = t[1 + idt]; in MatSolveTranspose_SeqBAIJ_4_inplace()
79 t[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
80 t[idx + 1] -= v[4] * s1 + v[5] * s2 + v[6] * s3 + v[7] * s4; in MatSolveTranspose_SeqBAIJ_4_inplace()
[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
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()
61 t[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3 + v[3] * s4 + v[4] * s5 + v[5] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
62 t[oidx + 1] -= v[6] * s1 + v[7] * s2 + v[8] * s3 + v[9] * s4 + v[10] * s5 + v[11] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
63 t[oidx + 2] -= v[12] * s1 + v[13] * s2 + v[14] * s3 + v[15] * s4 + v[16] * s5 + v[17] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
64 t[oidx + 3] -= v[18] * s1 + v[19] * s2 + v[20] * s3 + v[21] * s4 + v[22] * s5 + v[23] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
65 t[oidx + 4] -= v[24] * s1 + v[25] * s2 + v[26] * s3 + v[27] * s4 + v[28] * s5 + v[29] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
66 t[oidx + 5] -= v[30] * s1 + v[31] * s2 + v[32] * s3 + v[33] * s4 + v[34] * s5 + v[35] * s6; in MatSolveTranspose_SeqBAIJ_6_inplace()
70 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_6_inplace()
84 s2 = t[1 + idt]; in MatSolveTranspose_SeqBAIJ_6_inplace()
[all …]
|
| H A D | baijsolvtran3.c | 12 PetscScalar s1, s2, s3, x1, x2, x3, *x, *t; in MatSolveTranspose_SeqBAIJ_3_inplace() local
44 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_inplace()
52 t[oidx] -= v[0] * s1 + v[1] * s2 + v[2] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
53 t[oidx + 1] -= v[3] * s1 + v[4] * s2 + v[5] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
54 t[oidx + 2] -= v[6] * s1 + v[7] * s2 + v[8] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
58 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_3_inplace()
69 s2 = t[1 + idt]; in MatSolveTranspose_SeqBAIJ_3_inplace()
73 t[idx] -= v[0] * s1 + v[1] * s2 + v[2] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
74 t[idx + 1] -= v[3] * s1 + v[4] * s2 + v[5] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
75 t[idx + 2] -= v[6] * s1 + v[7] * s2 + v[8] * s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
[all …]
|
| 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
29 s2 = b[1 + idx]; 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()
48 x[1 + idx] = s2; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
60 s2 = x[1 + idt]; 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()
79 x[idt] = v[0] * s1 + v[5] * s2 + v[10] * s3 + v[15] * s4 + v[20] * s5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
80 x[1 + idt] = v[1] * s1 + v[6] * s2 + v[11] * s3 + v[16] * s4 + v[21] * s5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
81 x[2 + idt] = v[2] * s1 + v[7] * s2 + v[12] * s3 + v[17] * s4 + v[22] * s5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
82 x[3 + idt] = v[3] * s1 + v[8] * s2 + v[13] * s3 + v[18] * s4 + v[23] * s5; in MatSolve_SeqBAIJ_5_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
48 s2 = b[1 + idx]; 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()
64 x[1 + idx] = s2; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
76 s2 = x[1 + idt]; 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()
92 x[idt] = v[0] * s1 + v[4] * s2 + v[8] * s3 + v[12] * s4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
93 x[1 + idt] = v[1] * s1 + v[5] * s2 + v[9] * s3 + v[13] * s4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
94 x[2 + idt] = v[2] * s1 + v[6] * s2 + v[10] * s3 + v[14] * s4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
95 x[3 + idt] = v[3] * s1 + v[7] * s2 + v[11] * s3 + v[15] * s4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvtran2.c | 12 PetscScalar s1, s2, x1, x2, *x, *t; in MatSolveTranspose_SeqBAIJ_2_inplace() local
42 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2_inplace()
49 t[oidx] -= v[0] * s1 + v[1] * s2; in MatSolveTranspose_SeqBAIJ_2_inplace()
50 t[oidx + 1] -= v[2] * s1 + v[3] * s2; in MatSolveTranspose_SeqBAIJ_2_inplace()
54 t[1 + idx] = s2; in MatSolveTranspose_SeqBAIJ_2_inplace()
64 s2 = t[1 + idt]; in MatSolveTranspose_SeqBAIJ_2_inplace()
67 t[idx] -= v[0] * s1 + v[1] * s2; in MatSolveTranspose_SeqBAIJ_2_inplace()
68 t[idx + 1] -= v[2] * s1 + v[3] * s2; in MatSolveTranspose_SeqBAIJ_2_inplace()
99 PetscScalar s1, s2, x1, x2, *x, *t; in MatSolveTranspose_SeqBAIJ_2() local
128 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2()
[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
100 s2 = b[1 + idx]; 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()
126 t[1 + idx] = s2; in MatSolve_SeqBAIJ_7_inplace()
140 s2 = t[1 + idt]; in MatSolve_SeqBAIJ_7_inplace()
156 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()
166 …x[idc] = t[idt] = v[0] * s1 + v[7] * s2 + v[14] * s3 + v[21] * s4 + v[28] * s5 + v[35] * s6 + v[42… in MatSolve_SeqBAIJ_7_inplace()
167 …x[1 + idc] = t[1 + idt] = v[1] * s1 + v[8] * s2 + v[15] * s3 + v[22] * s4 + v[29] * s5 + v[36] * s… in MatSolve_SeqBAIJ_7_inplace()
168 …x[2 + idc] = t[2 + idt] = v[2] * s1 + v[9] * s2 + v[16] * s3 + v[23] * s4 + v[30] * s5 + v[37] * s… in MatSolve_SeqBAIJ_7_inplace()
169 …x[3 + idc] = t[3 + idt] = v[3] * s1 + v[10] * s2 + v[17] * s3 + v[24] * s4 + v[31] * s5 + v[38] * … in MatSolve_SeqBAIJ_7_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
30 s2 = b[1 + idx]; 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()
52 x[1 + idx] = s2; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
65 s2 = x[1 + idt]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
79 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()
87 x[idt] = v[0] * s1 + v[6] * s2 + v[12] * s3 + v[18] * s4 + v[24] * s5 + v[30] * s6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
88 x[1 + idt] = v[1] * s1 + v[7] * s2 + v[13] * s3 + v[19] * s4 + v[25] * s5 + v[31] * s6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
89 x[2 + idt] = v[2] * s1 + v[8] * s2 + v[14] * s3 + v[20] * s4 + v[26] * s5 + v[32] * s6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
90 x[3 + idt] = v[3] * s1 + v[9] * s2 + v[15] * s3 + v[21] * s4 + v[27] * s5 + v[33] * s6; 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
31 s2 = b[1 + idx]; 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()
56 x[1 + idx] = s2; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
70 s2 = x[1 + idt]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
86 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()
95 …x[idt] = v[0] * s1 + v[7] * s2 + v[14] * s3 + v[21] * s4 + v[28] * s5 + v[35] * s6 + v[42] * s… in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
96 …x[1 + idt] = v[1] * s1 + v[8] * s2 + v[15] * s3 + v[22] * s4 + v[29] * s5 + v[36] * s6 + v[43] * s… in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
97 …x[2 + idt] = v[2] * s1 + v[9] * s2 + v[16] * s3 + v[23] * s4 + v[30] * s5 + v[37] * s6 + v[44] * s… in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
98 …x[3 + idt] = v[3] * s1 + v[10] * s2 + v[17] * s3 + v[24] * s4 + v[31] * s5 + v[38] * s6 + v[45] * … in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 11 Vec Ax = NULL, Bx = NULL, s1 = NULL, s2 = NULL, Ay = NULL, By = NULL; in MatMultEqual_Private() local
39 PetscCall(MatCreateVecs(B, &s2, &Bx)); in MatMultEqual_Private()
42 PetscCall(MatCreateVecs(B, &Bx, &s2)); in MatMultEqual_Private()
46 PetscCall(VecDuplicate(s2, &By)); in MatMultEqual_Private()
62 PetscCall(VecCopy(By, s2)); in MatMultEqual_Private()
64 Badd = s2; in MatMultEqual_Private()
70 PetscCall(MatMultTransposeAdd(B, Bx, Badd, s2)); in MatMultEqual_Private()
73 PetscCall(MatMultTranspose(B, Bx, s2)); in MatMultEqual_Private()
78 PetscCall(MatMultHermitianTransposeAdd(B, Bx, Badd, s2)); in MatMultEqual_Private()
81 PetscCall(MatMultHermitianTranspose(B, Bx, s2)); in MatMultEqual_Private()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex53.c | 11 Vec xx, yy, s1, s2; in main() local
43 PetscCall(VecDuplicate(xx, &s2)); in main()
65 PetscCall(MatMult(B, xx, s2)); in main()
66 PetscCall(VecAXPY(s2, -1.0, s1)); in main()
67 PetscCall(VecNorm(s2, NORM_2, &rnorm)); in main()
76 PetscCall(MatMultAdd(B, xx, yy, s2)); in main()
77 PetscCall(VecAXPY(s2, -1.0, s1)); in main()
78 PetscCall(VecNorm(s2, NORM_2, &rnorm)); in main()
86 PetscCall(MatMultTranspose(B, xx, s2)); in main()
88 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
[all …]
|
| H A D | ex48.c | 8 Vec xx, s1, s2, yy; in main() local
31 PetscCall(VecDuplicate(xx, &s2)); in main()
109 PetscCall(VecSet(s2, 0.0)); in main()
111 PetscCall(MatMultAdd(A, xx, s2, s2)); in main()
113 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
183 PetscCall(MatSolveAdd(B, xx, yy, s2)); in main()
186 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
195 PetscCall(VecCopy(s2, s1)); in main()
196 PetscCall(MatSolveAdd(B, xx, s2, s2)); in main()
199 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
[all …]
|
| H A D | ex91.c | 13 Vec xx, s1, s2; in main() local
87 PetscCall(VecDuplicate(xx, &s2)); in main()
91 PetscCall(MatMult(sA, xx, s2)); in main()
93 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
99 PetscCall(VecDestroy(&s2)); in main()
144 PetscCall(VecDuplicate(xx, &s2)); in main()
148 PetscCall(MatMult(submatsA[i], xx, s2)); in main()
150 PetscCall(VecNorm(s2, NORM_2, &s2norm)); in main()
156 PetscCall(VecDestroy(&s2)); in main()
168 PetscCall(VecDuplicate(xx, &s2)); in main()
[all …]
|
| /petsc/src/binding/petsc4py/src/ |
| H A D | cython.h | 19 const char *desc, *s1, *s2; in Cython_ImportFunction() local
23 s1 = desc; s2 = signature; in Cython_ImportFunction()
24 while (*s1 != '\0' && *s1 == *s2) { s1++; s2++; } in Cython_ImportFunction()
25 if (*s1 != *s2) { in Cython_ImportFunction()
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex24.c | 76 PetscSection s1 = NULL, s2 = NULL, tpws = NULL; in main() local
123 PetscCall(PetscSectionCreate(comm, &s2)); in main()
124 PetscCall(PetscPartitionerDMPlexPartition(part2, dm2, tpws, s2, &is2)); in main()
140 PetscCall(PetscSectionViewFromOptions(s2, NULL, "-seq_s2_view")); in main()
142 PetscCall(PetscSectionCompare(s1, s2, &flg)); in main()
155 PetscCall(PetscSectionDestroy(&s2)); in main()
203 PetscCall(PetscSectionCreate(comm, &s2)); in main()
204 PetscCall(PetscPartitionerDMPlexPartition(part2, dmdist2, tpws, s2, &is2)); in main()
220 PetscCall(PetscSectionViewFromOptions(s2, NULL, "-dist_s2_view")); in main()
222 PetscCall(PetscSectionCompare(s1, s2, &flg)); in main()
[all …]
|