| /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
24 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * 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()
40 x[idx] = s1; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
52 s1 = x[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
25 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * 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()
43 x[idx] = s1; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
56 s1 = x[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
23 s1 = v[0] * x1 + v[1] * x2 + v[2] * 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()
37 x[idx] = s1; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
48 s1 = x[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
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()
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()
46 x[idx] = s1; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace()
60 s1 = x[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
22 s1 = v[0] * x1 + v[1] * 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()
34 x[idx] = s1; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
44 s1 = x[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
78 s1 = v[0] * x1 + v[1] * 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
34 s1 = b[idx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
44 x[idx] = s1; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
53 s1 = x[idt]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
59 s1 -= v[0] * x1 + v[2] * 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
95 s1 = b[idx]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat1.c | 9 PetscScalar *x, *tmp, s1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering() local
26 s1 = tmp[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
27 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
28 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
29 tmp[i] = s1; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
37 s1 = tmp[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
38 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering()
57 PetscScalar s1, *x; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace() local
67 s1 = (*v++) * x[i]; in MatSolveTranspose_SeqBAIJ_1_NaturalOrdering_inplace()
70 while (nz--) x[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_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
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()
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()
49 x[idx] = s1; 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
47 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * 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()
65 t[idx] = s1; in MatSolveTranspose_SeqBAIJ_5_inplace()
78 s1 = t[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
32 s1 = b[idx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
45 x[idx] = s1; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
55 s1 = x[idt]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
63 s1 -= v[0] * x1 + v[3] * x2 + v[6] * 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
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()
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()
73 t[idx] = s1; 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
45 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * 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()
61 t[idx] = s1; in MatSolveTranspose_SeqBAIJ_4_inplace()
73 s1 = t[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
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()
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()
69 t[idx] = s1; in MatSolveTranspose_SeqBAIJ_6_inplace()
83 s1 = t[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
43 s1 = v[0] * x1 + v[1] * x2 + v[2] * 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()
57 t[idx] = s1; in MatSolveTranspose_SeqBAIJ_3_inplace()
68 s1 = t[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 | baijsolvtran1.c | 11 PetscScalar *x, *tmp, s1; in MatSolveTranspose_SeqBAIJ_1() local
33 s1 = tmp[i]; in MatSolveTranspose_SeqBAIJ_1()
34 s1 *= v[nz]; /* multiply by inverse of diagonal entry */ in MatSolveTranspose_SeqBAIJ_1()
35 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
36 tmp[i] = s1; in MatSolveTranspose_SeqBAIJ_1()
44 s1 = tmp[i]; in MatSolveTranspose_SeqBAIJ_1()
45 for (j = 0; j < nz; j++) tmp[vi[j]] -= s1 * v[j]; in MatSolveTranspose_SeqBAIJ_1()
68 PetscScalar s1, *x, *t; in MatSolveTranspose_SeqBAIJ_1_inplace() local
88 s1 = (*v++) * t[i]; in MatSolveTranspose_SeqBAIJ_1_inplace()
91 while (nz--) t[*vi++] -= (*v++) * s1; in MatSolveTranspose_SeqBAIJ_1_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
99 s1 = b[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()
125 t[idx] = s1; in MatSolve_SeqBAIJ_7_inplace()
139 s1 = t[idt]; in MatSolve_SeqBAIJ_7_inplace()
155 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()
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 | baijsolvnat5.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, x1, x2, x3, x4, x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() local
28 s1 = b[idx]; 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()
47 x[idx] = s1; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
59 s1 = x[idt]; 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()
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
47 s1 = b[idx]; 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()
63 x[idx] = s1; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
75 s1 = x[idt]; 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()
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
41 s1 = v[0] * x1 + v[1] * 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()
53 t[idx] = s1; in MatSolveTranspose_SeqBAIJ_2_inplace()
63 s1 = t[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
127 s1 = v[0] * x1 + v[1] * x2; in MatSolveTranspose_SeqBAIJ_2()
[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
29 s1 = b[idx]; 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()
51 x[idx] = s1; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
64 s1 = x[idt]; 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()
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
30 s1 = b[idx]; 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()
55 x[idx] = s1; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
69 s1 = x[idt]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
85 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()
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/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/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
38 PetscCall(MatCreateVecs(A, &s1, &Ax)); in MatMultEqual_Private()
41 PetscCall(MatCreateVecs(A, &Ax, &s1)); in MatMultEqual_Private()
45 PetscCall(VecDuplicate(s1, &Ay)); in MatMultEqual_Private()
61 PetscCall(VecCopy(Ay, s1)); in MatMultEqual_Private()
63 Aadd = s1; in MatMultEqual_Private()
69 PetscCall(MatMultTransposeAdd(A, Ax, Aadd, s1)); in MatMultEqual_Private()
72 PetscCall(MatMultTranspose(A, Ax, s1)); in MatMultEqual_Private()
77 PetscCall(MatMultHermitianTransposeAdd(A, Ax, Aadd, s1)); in MatMultEqual_Private()
80 PetscCall(MatMultHermitianTranspose(A, Ax, s1)); in MatMultEqual_Private()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex48.c | 8 Vec xx, s1, s2, yy; in main() local
30 PetscCall(VecDuplicate(xx, &s1)); in main()
110 PetscCall(MatMult(A, xx, s1)); in main()
112 PetscCall(VecNorm(s1, NORM_2, &s1norm)); in main()
168 PetscCall(MatBackwardSolve(Fact, xx, s1)); in main()
170 PetscCall(VecScale(s1, -1.0)); in main()
171 PetscCall(MatMultAdd(A, s1, yy, yy)); in main()
184 PetscCall(MatSolveAdd(A, xx, yy, s1)); in main()
185 PetscCall(VecNorm(s1, NORM_2, &s1norm)); in main()
194 PetscCall(VecSetRandom(s1, rdm)); in main()
[all …]
|
| H A D | ex53.c | 11 Vec xx, yy, s1, s2; in main() local
42 PetscCall(VecDuplicate(xx, &s1)); in main()
64 PetscCall(MatMult(A, xx, s1)); in main()
66 PetscCall(VecAXPY(s2, -1.0, s1)); in main()
75 PetscCall(MatMultAdd(A, xx, yy, s1)); in main()
77 PetscCall(VecAXPY(s2, -1.0, s1)); in main()
85 PetscCall(MatMultTranspose(A, xx, s1)); in main()
87 PetscCall(VecNorm(s1, NORM_2, &s1norm)); in main()
96 PetscCall(MatMultTransposeAdd(A, xx, yy, s1)); in main()
98 PetscCall(VecNorm(s1, NORM_2, &s1norm)); in main()
[all …]
|