| /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
26 s3 = v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * 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()
42 x[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
54 s3 = x[2 + 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
27 s3 = v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * 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()
45 x[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
58 s3 = x[2 + 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
25 s3 = v[6] * x1 + v[7] * x2 + v[8] * 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()
39 x[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
50 s3 = x[2 + 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
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()
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()
48 x[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
62 s3 = x[2 + idt]; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
[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
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()
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()
51 x[2 + idx] = s3; 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
49 s3 = v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * 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()
67 t[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_5_inplace()
80 s3 = t[2 + 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
34 s3 = b[2 + idx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
42 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
47 x[2 + idx] = s3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
57 s3 = x[2 + idt]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
65 s3 -= v[2] * x1 + v[5] * x2 + v[8] * 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
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()
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()
75 t[2 + idx] = s3; 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
47 s3 = v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * 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()
63 t[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_4_inplace()
75 s3 = t[2 + 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
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()
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()
71 t[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_6_inplace()
85 s3 = t[2 + 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
45 s3 = v[6] * x1 + v[7] * x2 + v[8] * 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()
59 t[2 + idx] = s3; in MatSolveTranspose_SeqBAIJ_3_inplace()
70 s3 = t[2 + 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
30 s3 = b[2 + idx]; 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()
49 x[2 + idx] = s3; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
61 s3 = x[2 + idt]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
73 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * 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
49 s3 = b[2 + idx]; 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()
65 x[2 + idx] = s3; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
77 s3 = x[2 + idt]; 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()
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 | baijsolvnat6.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() local
31 s3 = b[2 + idx]; 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()
53 x[2 + idx] = s3; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
66 s3 = x[2 + idt]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
80 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()
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 | 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
101 s3 = b[2 + idx]; 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()
127 t[2 + idx] = s3; in MatSolve_SeqBAIJ_7_inplace()
141 s3 = t[2 + idt]; in MatSolve_SeqBAIJ_7_inplace()
157 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()
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 | 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
32 s3 = b[2 + idx]; 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()
57 x[2 + idx] = s3; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
71 s3 = x[2 + idt]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
87 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()
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 …]
|
| H A D | baijsolvnat15.c | 12 PetscScalar s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2() local
46 s3 = b[2 + idt]; 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()
97 x[2 + idt] = s3; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
119 s3 = x[2 + idt]; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
153 …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()
170 …x[idt] = v[0] * s1 + v[15] * s2 + v[30] * s3 + v[45] * s4 + v[60] * s5 + v[75] * s6 + v[90] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
171 …x[1 + idt] = v[1] * s1 + v[16] * s2 + v[31] * s3 + v[46] * s4 + v[61] * s5 + v[76] * s6 + v[91] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
172 …x[2 + idt] = v[2] * s1 + v[17] * s2 + v[32] * s3 + v[47] * s4 + v[62] * s5 + v[77] * s6 + v[92] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
173 …x[3 + idt] = v[3] * s1 + v[18] * s2 + v[33] * s3 + v[48] * s4 + v[63] * s5 + v[78] * s6 + v[93] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
[all …]
|
| /petsc/config/PETSc/options/ |
| H A D | petscclone.py | 28 …(o3, e3, s3) = self.executeShellCommand([self.sourceControl.git, 'log', '-1', '--pretty=format:%ci…
31 if s2 or s3 or s4:
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 118 Vec Ax, Bx, Cx, s1, s2, s3; in MatMatMultEqual_Private() local
168 PetscCall(MatCreateVecs(C, &Cx, &s3)); in MatMatMultEqual_Private()
185 PetscCall(MatMult(C, Cx, s3)); in MatMatMultEqual_Private()
189 PetscCall(VecNorm(s3, NORM_INFINITY, &r1)); in MatMatMultEqual_Private()
191 PetscCall(VecAXPY(s2, none, s3)); in MatMatMultEqual_Private()
207 PetscCall(VecDestroy(&s3)); in MatMatMultEqual_Private()
|
| /petsc/src/mat/impls/sell/mpi/ |
| H A D | mpisell.c | 801 PetscInt s1, s2, s3; in MatDiagonalScale_MPISELL() local
804 PetscCall(MatGetLocalSize(mat, &s2, &s3)); in MatDiagonalScale_MPISELL()
807 …PetscCheck(s1 == s3, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "right vector non-conforming local size"); in MatDiagonalScale_MPISELL()
|
| /petsc/src/mat/impls/dense/mpi/ |
| H A D | mpidense.c | 978 PetscInt i, j, s2a, s3a, s2, s3, m = mdn->A->rmap->n, n = mdn->A->cmap->n, lda; in MatDiagonalScale_MPIDense() local
983 PetscCall(MatGetLocalSize(A, &s2, &s3)); in MatDiagonalScale_MPIDense()
1003 …k(s3a == s3, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Right scaling vec non-conforming local size, %" … in MatDiagonalScale_MPIDense()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 514 PetscReal c1, s1, c2, s2, c3, s3; in DMPlexBasisTransformSetUp_Rotation_Internal() local
535 s3 = PetscSinReal(rc->gamma); in DMPlexBasisTransformSetUp_Rotation_Internal()
536 rc->R[0] = c1 * c3 - c2 * s1 * s3; in DMPlexBasisTransformSetUp_Rotation_Internal()
537 rc->R[1] = c3 * s1 + c1 * c2 * s3; in DMPlexBasisTransformSetUp_Rotation_Internal()
538 rc->R[2] = s2 * s3; in DMPlexBasisTransformSetUp_Rotation_Internal()
539 rc->R[3] = -c1 * s3 - c2 * c3 * s1; in DMPlexBasisTransformSetUp_Rotation_Internal()
540 rc->R[4] = c1 * c2 * c3 - s1 * s3; in DMPlexBasisTransformSetUp_Rotation_Internal()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 1571 PetscInt s1, s2, s3; in MatDiagonalScale_MPIBAIJ() local
1574 PetscCall(MatGetLocalSize(mat, &s2, &s3)); in MatDiagonalScale_MPIBAIJ()
1577 …PetscCheck(s1 == s3, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "right vector non-conforming local size"); in MatDiagonalScale_MPIBAIJ()
|
| /petsc/doc/manual/ |
| H A D | other.md | 312 diagonal information `s1,s2,s3,...` indicates which (block) diagonals
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiaij.c | 1982 PetscInt s1, s2, s3; in MatDiagonalScale_MPIAIJ() local
1985 PetscCall(MatGetLocalSize(mat, &s2, &s3)); in MatDiagonalScale_MPIAIJ()
1988 …PetscCheck(s1 == s3, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "right vector non-conforming local size"); in MatDiagonalScale_MPIAIJ()
|