Home
last modified time | relevance | path

Searched refs:v (Results 1 – 25 of 968) sorted by relevance

12345678910>>...39

/petsc/src/mat/impls/baij/seq/
H A Dbaijsolvnat15.c11 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2() local
40 v = aa + bs2 * ai[i]; in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
77 …v[0] * x1 + v[15] * x2 + v[30] * x3 + v[45] * x4 + v[60] * x5 + v[75] * x6 + v[90] * x7 + v[105] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
78 …v[1] * x1 + v[16] * x2 + v[31] * x3 + v[46] * x4 + v[61] * x5 + v[76] * x6 + v[91] * x7 + v[106] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
79 …v[2] * x1 + v[17] * x2 + v[32] * x3 + v[47] * x4 + v[62] * x5 + v[77] * x6 + v[92] * x7 + v[107] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
80 …v[3] * x1 + v[18] * x2 + v[33] * x3 + v[48] * x4 + v[63] * x5 + v[78] * x6 + v[93] * x7 + v[108] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
81 …v[4] * x1 + v[19] * x2 + v[34] * x3 + v[49] * x4 + v[64] * x5 + v[79] * x6 + v[94] * x7 + v[109] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
82 …v[5] * x1 + v[20] * x2 + v[35] * x3 + v[50] * x4 + v[65] * x5 + v[80] * x6 + v[95] * x7 + v[110] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
83 …v[6] * x1 + v[21] * x2 + v[36] * x3 + v[51] * x4 + v[66] * x5 + v[81] * x6 + v[96] * x7 + v[111] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
84 …v[7] * x1 + v[22] * x2 + v[37] * x3 + v[52] * x4 + v[67] * x5 + v[82] * x6 + v[97] * x7 + v[112] *… in MatSolve_SeqBAIJ_15_NaturalOrdering_ver2()
[all …]
H A Dbaijsolvtrannat7.c8 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace() local
18 v = aa + 49 * diag[i]; 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()
34 v += 49; in MatSolveTranspose_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
H A Dbaijsolvtran7.c11 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_7_inplace() local
42 v = aa + 49 * diag[i]; 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()
58 v += 49; in MatSolveTranspose_SeqBAIJ_7_inplace()
[all …]
H A Dbaijsolvnat7.c9 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() local
26 v = aa + 49 * ai[i]; 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()
53 v += 49; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
[all …]
H A Dbaijsolvtrannat6.c8 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_6_NaturalOrdering_inplace() local
18 v = aa + 36 * diag[i]; 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()
32 v += 36; 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()
[all …]
H A Dbaijsolvnat6.c9 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() local
25 v = aa + 36 * ai[i]; 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()
49 v += 36; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
60 v = aa + 36 * diag[i] + 36; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
[all …]
H A Dbaijsolvtrannat5.c8 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace() local
18 v = aa + 25 * diag[i]; 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()
30 v += 25; 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()
[all …]
H A Dbaijsolvtran6.c11 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_6_inplace() local
41 v = aa + 36 * diag[i]; 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()
55 v += 36; 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()
[all …]
H A Dbaijsolvnat5.c9 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() local
24 v = aa + 25 * ai[i]; 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()
45 v += 25; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
55 v = aa + 25 * diag[i] + 25; 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()
[all …]
H A Dbaijsolv.c12 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_N_inplace() local
30 v = aa + bs2 * ai[i]; in MatSolve_SeqBAIJ_N_inplace()
36 PetscKernel_v_gets_v_minus_A_times_w(bs, s, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
37 v += bs2; in MatSolve_SeqBAIJ_N_inplace()
43 v = aa + bs2 * (a->diag[i] + 1); in MatSolve_SeqBAIJ_N_inplace()
48 PetscKernel_v_gets_v_minus_A_times_w(bs, ls, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
49 v += bs2; in MatSolve_SeqBAIJ_N_inplace()
70 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_7_inplace() local
95 v = aa + 49 * ai[i]; 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()
[all …]
H A Dbaijsolvtran5.c11 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_5_inplace() local
40 v = aa + 25 * diag[i]; 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()
52 v += 25; 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()
[all …]
H A Dbaijsolvtrannat4.c8 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace() local
18 v = aa + 16 * diag[i]; 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()
28 v += 16; 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()
[all …]
H A Dbaij2.c241 const MatScalar *v; in MatMult_SeqBAIJ_1() local
262 v = a->a + ii[0]; in MatMult_SeqBAIJ_1()
266 PetscPrefetchBlock(v + 1 * n, 1 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqBAIJ_1()
268 PetscSparseDensePlusDot(sum, x, v, idx, n); in MatMult_SeqBAIJ_1()
287 const MatScalar *v; in MatMult_SeqBAIJ_2() local
296 v = a->a; in MatMult_SeqBAIJ_2()
314 PetscPrefetchBlock(v + 4 * n, 4 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqBAIJ_2()
319 sum1 += v[0] * x1 + v[2] * x2; in MatMult_SeqBAIJ_2()
320 sum2 += v[1] * x1 + v[3] * x2; in MatMult_SeqBAIJ_2()
321 v += 4; in MatMult_SeqBAIJ_2()
[all …]
H A Dbaijsolvnat4.c32 const MatScalar *v; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
43 v = aa + 16 * ai[i]; 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()
61 v += 16; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
72 v = aa + ai16 + 16; 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()
[all …]
H A Dbaijsolvtran4.c11 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_4_inplace() local
39 v = aa + 16 * diag[i]; 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()
49 v += 16; 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()
[all …]
H A Dbaijsolvtrannat3.c8 const MatScalar *aa = a->a, *v; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
18 v = aa + 9 * diag[i]; 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()
26 v += 9; 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()
35 v += 9; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
[all …]
H A Dbaijsolvnat14.c12 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_14_NaturalOrdering() local
23 v = aa + bs2 * ai[i]; in MatSolve_SeqBAIJ_14_NaturalOrdering()
45 x[idt] -= v[0] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
46 x[1 + idt] -= v[1] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
47 x[2 + idt] -= v[2] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
48 x[3 + idt] -= v[3] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
49 x[4 + idt] -= v[4] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
50 x[5 + idt] -= v[5] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
51 x[6 + idt] -= v[6] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
52 x[7 + idt] -= v[7] * xv; in MatSolve_SeqBAIJ_14_NaturalOrdering()
[all …]
H A Dbaijsolvnat3.c13 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
28 v = aa + 9 * ai[i]; 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()
43 v += 9; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
51 v = aa + 9 * diag[i] + 9; 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()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij2.c292 const MatScalar *v; in MatMult_SeqSBAIJ_2() local
303 v = a->a; in MatMult_SeqSBAIJ_2()
314 z[2 * i] += v[0] * x1 + v[2] * x2; in MatMult_SeqSBAIJ_2()
315 z[2 * i + 1] += v[2] * x1 + v[3] * x2; in MatMult_SeqSBAIJ_2()
316 v += 4; in MatMult_SeqSBAIJ_2()
320 PetscPrefetchBlock(v + 4 * n, 4 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_2()
324 z[cval] += v[0] * x1 + v[1] * x2; in MatMult_SeqSBAIJ_2()
325 z[cval + 1] += v[2] * x1 + v[3] * x2; in MatMult_SeqSBAIJ_2()
327 z[2 * i] += v[0] * x[cval] + v[2] * x[cval + 1]; in MatMult_SeqSBAIJ_2()
328 z[2 * i + 1] += v[1] * x[cval] + v[3] * x[cval + 1]; in MatMult_SeqSBAIJ_2()
[all …]
/petsc/include/petsc/private/kernels/
H A Dblockmatmult.h5 #define PetscKernel_v_gets_A_times_w_1_exp(v, A, w, exp) \ argument
7 v[0] exp A[0] * w[0]; \
10 #define PetscKernel_v_gets_A_times_w_2_exp(v, A, w, exp) \ argument
12 v[0] exp A[0] * w[0] + A[2] * w[1]; \
13 v[1] exp A[1] * w[0] + A[3] * w[1]; \
16 #define PetscKernel_v_gets_A_times_w_3_exp(v, A, w, exp) \ argument
18 v[0] exp A[0] * w[0] + A[3] * w[1] + A[6] * w[2]; \
19 v[1] exp A[1] * w[0] + A[4] * w[1] + A[7] * w[2]; \
20 v[2] exp A[2] * w[0] + A[5] * w[1] + A[8] * w[2]; \
23 #define PetscKernel_v_gets_A_times_w_4_exp(v, A, w, exp) \ argument
[all …]
/petsc/src/sys/tests/linter/output/
H A DtestValidLogicalCollective.patch7 PetscViewer v; /* dummy variable to satisfy the PetscObject for the following */
11 - PetscValidLogicalCollectiveInt(v, d, 2);
12 + PetscValidLogicalCollectiveBool(v, d, 4);
13 PetscValidLogicalCollectiveEnum(v, e, 3);
14 PetscValidLogicalCollectiveMPIInt(v, f, 4);
17 PetscValidLogicalCollectiveInt(v, d, 2);
18 - PetscValidLogicalCollectiveEnum(v, e, 3);
19 + PetscValidLogicalCollectiveScalar(v, e, 5);
20 PetscValidLogicalCollectiveMPIInt(v, f, 4);
21 PetscValidLogicalCollectiveScalar(v, a, 5);
[all …]
/petsc/src/mat/ftn-mod/
H A Dpetscmatmod.F9060 subroutine MatDenseGetArrayF901d(v, array, ierr) argument
63 Mat v
64 call MatDenseGetArray(v, array, ierr)
70 subroutine MatDenseRestoreArrayF901d(v, array, ierr) argument
73 Mat v
74 call MatDenseRestoreArray(v, array, ierr)
80 subroutine MatDenseGetArrayReadF901d(v, array, ierr) argument
83 Mat v
84 call MatDenseGetArrayRead(v, array, ierr)
90 subroutine MatDenseRestoreArrayReadF901d(v, array, ierr) argument
[all …]
/petsc/src/vec/vec/tutorials/
H A Dex1.c16 PetscReal norm, v, v1, v2, maxval; in main() local
78 v = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
79 if (v > -PETSC_SMALL && v < PETSC_SMALL) v = 0.0; in main()
80 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "VecScale %g\n", (double)v)); in main()
84 v = norm - 2.0 * PetscSqrtReal((PetscReal)n); in main()
85 if (v > -PETSC_SMALL && v < PETSC_SMALL) v = 0.0; in main()
86 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "VecCopy %g\n", (double)v)); in main()
90 v = norm - 8.0 * PetscSqrtReal((PetscReal)n); in main()
91 if (v > -PETSC_SMALL && v < PETSC_SMALL) v = 0.0; in main()
92 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "VecAXPY %g\n", (double)v)); in main()
[all …]
H A Dex20f90.F9018 PetscReal norm, v, v1, v2, tol
99 v = abs(norm - 2.0*sqrt(nfloat))
100 if (v > -tol .and. v < tol) v = 0.0
101 if (rank == 0) write (6, 130) v
106 v = abs(norm - 2.0*sqrt(nfloat))
107 if (v > -tol .and. v < tol) v = 0.0
108 if (rank == 0) write (6, 140) v
113 v = abs(norm - 8.0*sqrt(nfloat))
114 if (v > -tol .and. v < tol) v = 0.0
115 if (rank == 0) write (6, 150) v
[all …]
H A Dex1f90.F9017 PetscReal norm, v, v1, v2
97 v = abs(norm - 2.0*sqrt(nfloat))
98 if (v > -1.d-10 .and. v < 1.d-10) v = 0.0
99 if (rank == 0) write (6, 130) v
104 v = abs(norm - 2.0*sqrt(nfloat))
105 if (v > -1.d-10 .and. v < 1.d-10) v = 0.0
106 if (rank == 0) write (6, 140) v
111 v = abs(norm - 8.0*sqrt(nfloat))
112 if (v > -1.d-10 .and. v < 1.d-10) v = 0.0
113 if (rank == 0) write (6, 150) v
[all …]

12345678910>>...39