| /petsc/src/ksp/pc/impls/pbjacobi/ |
| H A D | pbjacobi.c | 11 PetscScalar *yy, x0, x1, x2, x3, x4, x5, x6; in PCApply_PBJacobi() local
34 x2 = xx[3 * i + 2]; in PCApply_PBJacobi()
36 yy[3 * i] = diag[0] * x0 + diag[3] * x1 + diag[6] * x2; in PCApply_PBJacobi()
37 yy[3 * i + 1] = diag[1] * x0 + diag[4] * x1 + diag[7] * x2; in PCApply_PBJacobi()
38 yy[3 * i + 2] = diag[2] * x0 + diag[5] * x1 + diag[8] * x2; in PCApply_PBJacobi()
46 x2 = xx[4 * i + 2]; in PCApply_PBJacobi()
49 yy[4 * i] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_PBJacobi()
50 yy[4 * i + 1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3; in PCApply_PBJacobi()
51 yy[4 * i + 2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3; in PCApply_PBJacobi()
52 yy[4 * i + 3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3; in PCApply_PBJacobi()
[all …]
|
| /petsc/src/ksp/pc/impls/vpbjacobi/ |
| H A D | vpbjacobi.c | 11 PetscScalar *yy, x0, x1, x2, x3, x4, x5, x6; in PCApply_VPBJacobi() local
34 x2 = xx[ncnt + 2]; in PCApply_VPBJacobi()
35 yy[ncnt] = diag[0] * x0 + diag[3] * x1 + diag[6] * x2; in PCApply_VPBJacobi()
36 yy[ncnt + 1] = diag[1] * x0 + diag[4] * x1 + diag[7] * x2; in PCApply_VPBJacobi()
37 yy[ncnt + 2] = diag[2] * x0 + diag[5] * x1 + diag[8] * x2; in PCApply_VPBJacobi()
42 x2 = xx[ncnt + 2]; in PCApply_VPBJacobi()
44 yy[ncnt] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_VPBJacobi()
45 yy[ncnt + 1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3; in PCApply_VPBJacobi()
46 yy[ncnt + 2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3; in PCApply_VPBJacobi()
47 yy[ncnt + 3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3; in PCApply_VPBJacobi()
[all …]
|
| /petsc/src/vec/vec/tutorials/ |
| H A D | ex19.c | 8 Vec x1, x2, *x3ts, *x4ts; in main() local
39 PetscCall(VecCreate(comm, &x2)); in main()
40 PetscCall(PetscObjectSetName((PetscObject)x2, "x2")); in main()
41 PetscCall(VecSetSizes(x2, PETSC_DECIDE, n)); in main()
42 PetscCall(VecSetBlockSize(x2, 2)); in main()
43 PetscCall(VecSetFromOptions(x2)); in main()
44 PetscCall(VecSetRandom(x2, rand)); in main()
45 PetscCall(VecView(x2, viewer)); in main()
60 PetscCall(VecDuplicateVecs(x2, n_timesteps, &x4ts)); in main()
89 PetscCall(VecDuplicate(x2, &x2r)); in main()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 286 PetscScalar x1, x2; in MatMult_SeqBAIJ_2() local
318 x2 = xb[1]; 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()
337 PetscScalar *z = NULL, sum1, sum2, sum3, x1, x2, x3, *zarray; in MatMult_SeqBAIJ_3() local
375 x2 = xb[1]; in MatMult_SeqBAIJ_3()
378 sum1 += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMult_SeqBAIJ_3()
379 sum2 += v[1] * x1 + v[4] * x2 + v[7] * x3; in MatMult_SeqBAIJ_3()
380 sum3 += v[2] * x1 + v[5] * x2 + v[8] * x3; in MatMult_SeqBAIJ_3()
398 PetscScalar *z = NULL, sum1, sum2, sum3, sum4, x1, x2, x3, x4, *zarray; in MatMult_SeqBAIJ_4() local
[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
36 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
42 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
43 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
44 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
67 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
71 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
72 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvnat3.c | 14 PetscScalar *x, s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local
38 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
42 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
61 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
63 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
64 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
65 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
89 PetscScalar s1, s2, s3, x1, x2, x3; in MatSolve_SeqBAIJ_3_NaturalOrdering() local
[all …]
|
| H A D | 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
109 x2 = t[1 + idx]; in MatSolve_SeqBAIJ_7_inplace()
115 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_inplace()
116 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_inplace()
117 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_inplace()
118 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_inplace()
119 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_inplace()
120 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_inplace()
121 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_inplace()
149 x2 = t[1 + idx]; in MatSolve_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvnat2.c | 17 PetscScalar *x, s1, s2, x1, x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
39 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
40 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
41 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
58 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
59 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
60 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace()
80 PetscScalar *x, s1, s2, x1, x2; in MatSolve_SeqBAIJ_2_NaturalOrdering() local
102 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
103 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering()
[all …]
|
| H A D | baijsolvnat4.c | 31 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
54 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
57 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
58 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
59 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
60 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
82 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
85 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
86 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
87 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
[all …]
|
| H A D | baijsolvnat6.c | 10 PetscScalar *x, s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() local
38 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
43 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
44 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
45 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
46 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
47 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
48 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
73 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
78 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace()
[all …]
|
| H A D | 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
40 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
46 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
47 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
48 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
49 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
50 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
51 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
52 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace()
79 x2 = x[1 + idx]; in MatSolve_SeqBAIJ_7_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
21 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
25 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
26 s2 = v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
27 s3 = v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
28 s4 = v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
29 s5 = v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering_inplace()
83 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering() local
95 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
99 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatSolveTranspose_SeqBAIJ_5_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat4.c | 9 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace() local
21 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
24 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
25 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
26 s3 = v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
27 s4 = v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering_inplace()
77 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4, *x; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering() local
89 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
92 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
93 s2 = v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatSolveTranspose_SeqBAIJ_4_NaturalOrdering()
[all …]
|
| H A D | baijsolvtrannat3.c | 9 PetscScalar s1, s2, s3, x1, x2, x3, *x; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace() local
21 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
23 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
24 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
25 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering_inplace()
71 PetscScalar s1, s2, s3, x1, x2, x3, *x; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering() local
83 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
85 s1 = v[0] * x1 + v[1] * x2 + v[2] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
86 s2 = v[3] * x1 + v[4] * x2 + v[5] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
87 s3 = v[6] * x1 + v[7] * x2 + v[8] * x3; in MatSolveTranspose_SeqBAIJ_3_NaturalOrdering()
|
| H A D | baijsolvtrannat2.c | 9 PetscScalar s1, s2, x1, x2, *x; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace() local
21 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
22 s1 = v[0] * x1 + v[1] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
23 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering_inplace()
65 PetscScalar s1, s2, x1, x2, *x; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering() local
77 x2 = x[1 + idx]; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
78 s1 = v[0] * x1 + v[1] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
79 s2 = v[2] * x1 + v[3] * x2; in MatSolveTranspose_SeqBAIJ_2_NaturalOrdering()
|
| /petsc/src/snes/interface/ |
| H A D | snesj.c | 50 Vec j1a, j2a, x2; in SNESComputeJacobianDefault() local
78 PetscCall(DMGetGlobalVector(snes->dm, &x2)); in SNESComputeJacobianDefault()
84 x2 = snes->vwork[2]; in SNESComputeJacobianDefault()
106 x2 = perturbed iterate, j2a = F(x2) in SNESComputeJacobianDefault()
109 PetscCall(VecCopy(x1, x2)); in SNESComputeJacobianDefault()
118 if (x2->ops->setvalues) PetscCall(VecSetValues(x2, 1, &i, &dx, ADD_VALUES)); in SNESComputeJacobianDefault()
120 PetscCall(VecGetArray(x2, &y)); in SNESComputeJacobianDefault()
122 PetscCall(VecRestoreArray(x2, &y)); in SNESComputeJacobianDefault()
127 PetscCall(VecAssemblyBegin(x2)); in SNESComputeJacobianDefault()
128 PetscCall(VecAssemblyEnd(x2)); in SNESComputeJacobianDefault()
[all …]
|
| /petsc/src/vec/vec/impls/seq/ |
| H A D | dvec2.c | 87 PetscScalar sum0 = 0.0, sum1 = 0.0, sum2 = 0.0, sum3, x0, x1, x2, x3; in VecMDot_Seq() local
105 x2 = x[2]; in VecMDot_Seq()
106 sum0 += x2 * PetscConj(yy0[2]); in VecMDot_Seq()
107 sum1 += x2 * PetscConj(yy1[2]); in VecMDot_Seq()
108 sum2 += x2 * PetscConj(yy2[2]); /* fall through */ in VecMDot_Seq()
130 x2 = x[2]; in VecMDot_Seq()
134 …sum0 += x0 * PetscConj(yy0[0]) + x1 * PetscConj(yy0[1]) + x2 * PetscConj(yy0[2]) + x3 * PetscConj(… in VecMDot_Seq()
136 …sum1 += x0 * PetscConj(yy1[0]) + x1 * PetscConj(yy1[1]) + x2 * PetscConj(yy1[2]) + x3 * PetscConj(… in VecMDot_Seq()
138 …sum2 += x0 * PetscConj(yy2[0]) + x1 * PetscConj(yy2[1]) + x2 * PetscConj(yy2[2]) + x3 * PetscConj(… in VecMDot_Seq()
154 x2 = x[2]; in VecMDot_Seq()
[all …]
|
| /petsc/src/ksp/pc/impls/mat/tests/ |
| H A D | ex1.c | 41 Vec b, x, x2; in TestPCMatVersusMat() local
42 Mat B, X, X2; in TestPCMatVersusMat() local
52 PetscCall(VecDuplicate(x, &x2)); in TestPCMatVersusMat()
56 PetscCall(MatDuplicate(A, MAT_DO_NOT_COPY_VALUES, &X2)); in TestPCMatVersusMat()
60 PetscCall(PCApply(pc, b, x2)); in TestPCMatVersusMat()
61 PetscCall(TestVecEquality(x, x2)); in TestPCMatVersusMat()
64 PetscCall(PCApplyTranspose(pc, b, x2)); in TestPCMatVersusMat()
65 PetscCall(TestVecEquality(x, x2)); in TestPCMatVersusMat()
68 PetscCall(PCMatApply(pc, B, X2)); in TestPCMatVersusMat()
69 PetscCall(TestMatEquality(X, X2)); in TestPCMatVersusMat()
[all …]
|
| /petsc/src/vec/vec/tests/ |
| H A D | ex60.c | 9 Vec x, x1, x2, x3; in main() local
28 PetscCall(VecCreateSeqCUDAWithArray(PETSC_COMM_SELF, bs, n, NULL, &x2)); in main()
34 PetscCall(VecCreateSeqKokkosWithArray(PETSC_COMM_SELF, bs, n, NULL, &x2)); in main()
40 PetscCall(VecCreateSeqHIPWithArray(PETSC_COMM_SELF, bs, n, NULL, &x2)); in main()
45 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, bs, n, NULL, &x2)); in main()
51 PetscCall(VecPlaceArray(x2, px + n)); in main()
54 PetscCall(VecSet(x2, 0.0)); in main()
57 PetscCall(VecResetArray(x2)); in main()
61 PetscCall(VecDestroy(&x2)); in main()
|
| H A D | ex39.c | 9 X1,X2,Y1,Y2,F1,F2 and does vector addition on the subvectors F1 = X1 + Y1, F2 = X2 + Y2 so
16 Vec X, Y, F, X1, Y1, X2, Y2, F1, F2; in main() local
39 /* Create subvectors X1,X2,Y1,Y2,F1,F2 */ in main()
43 PetscCall(VecDuplicate(X1, &X2)); in main()
55 PetscCall(VecPlaceArray(X2, x + N / 2)); in main()
63 PetscCall(VecWAXPY(F2, 1.0, X2, Y2)); in main()
67 PetscCall(VecResetArray(X2)); in main()
87 PetscCall(VecDestroy(&X2)); in main()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 290 PetscScalar *z, x1, x2, zero = 0.0; in MatMult_SeqSBAIJ_2() local
309 x2 = xb[1]; 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()
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()
344 PetscScalar *z, x1, x2, x3, zero = 0.0; in MatMult_SeqSBAIJ_3() local
363 x2 = xb[1]; in MatMult_SeqSBAIJ_3()
369 z[3 * i] += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMult_SeqSBAIJ_3()
370 z[3 * i + 1] += v[3] * x1 + v[4] * x2 + v[7] * x3; in MatMult_SeqSBAIJ_3()
[all …]
|
| H A D | sbaijfact2.c | 222 PetscScalar *x, x0, x1, x2, x3, x4, x5, x6, *t, *tp; in MatSolve_SeqSBAIJ_7_inplace() local
251 x2 = tp[2]; in MatSolve_SeqSBAIJ_7_inplace()
259 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3 + v[4] * x4 + v[5] * x5 + v[6] * x6; in MatSolve_SeqSBAIJ_7_inplace()
260 … tp[1] += v[7] * x0 + v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4 + v[12] * x5 + v[13] * x6; in MatSolve_SeqSBAIJ_7_inplace()
261 … tp[2] += v[14] * x0 + v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5 + v[20] * x6; in MatSolve_SeqSBAIJ_7_inplace()
262 … tp[3] += v[21] * x0 + v[22] * x1 + v[23] * x2 + v[24] * x3 + v[25] * x4 + v[26] * x5 + v[27] * x6; in MatSolve_SeqSBAIJ_7_inplace()
263 … tp[4] += v[28] * x0 + v[29] * x1 + v[30] * x2 + v[31] * x3 + v[32] * x4 + v[33] * x5 + v[34] * x6; in MatSolve_SeqSBAIJ_7_inplace()
264 … tp[5] += v[35] * x0 + v[36] * x1 + v[37] * x2 + v[38] * x3 + v[39] * x4 + v[40] * x5 + v[41] * x6; in MatSolve_SeqSBAIJ_7_inplace()
265 … tp[6] += v[42] * x0 + v[43] * x1 + v[44] * x2 + v[45] * x3 + v[46] * x4 + v[47] * x5 + v[48] * x6; in MatSolve_SeqSBAIJ_7_inplace()
274 tp[0] = d[0] * x0 + d[7] * x1 + d[14] * x2 + d[21] * x3 + d[28] * x4 + d[35] * x5 + d[42] * x6; in MatSolve_SeqSBAIJ_7_inplace()
[all …]
|
| /petsc/src/ts/tutorials/hybrid/output/ |
| H A D | ex1fd_1.out | 5 d x1(tf) /d x1(t0) = -0.237478 d x2(tf) / d x1(t0) = 0.031943
8 d x1(tf) /d x2(t0) = -0.319570 d x2(tf) / d x2(t0) = -0.237558
11 Final gradients: d x1(tf) /d p = -0.101525 d x2(tf) / d p = 0.070411
|
| /petsc/src/dm/tests/ |
| H A D | ex2k.kokkos.cxx | 55 static PetscErrorCode Update2(DM da, const PetscScalar *__restrict__ x2, PetscScalar *__restrict__ … in Update2() argument
64 #define X2(k, j, i) x2[(k - gzs) * gym * gxm + (j - gys) * gxm + (i - gxs)] in Update2() macro
70 …i++) Y2(k, j, i) = 6 * X2(k, j, i) - X2(k - 1, j, i) - X2(k, j - 1, i) - X2(k, j, i - 1) - X2(k + … in Update2()
76 #undef X2 in Update2()
92 const PetscScalar *x2; in main() local
130 PetscCall(VecGetArrayRead(x, &x2)); in main()
132 PetscCall(Update2(da, x2, y2, nwarm, nloop, &avgTime)); in main()
133 PetscCall(VecRestoreArrayRead(x, &x2)); in main()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex5.c | 9 Vec x, x2, f, u, b; in main() local
40 PetscCall(VecDuplicate(x, &x2)); in main()
55 PetscCall(PCApply(P, b, x2)); in main()
56 PetscCall(VecAXPY(x2, -1.0, x)); in main()
57 PetscCall(VecNorm(x2, NORM_2, &err)); in main()
66 PetscCall(VecDestroy(&x2)); in main()
|