| /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
27 x1 = xx[ncnt + 1]; in PCApply_VPBJacobi()
28 yy[ncnt] = diag[0] * x0 + diag[2] * x1; in PCApply_VPBJacobi()
29 yy[ncnt + 1] = diag[1] * x0 + diag[3] * x1; in PCApply_VPBJacobi()
33 x1 = xx[ncnt + 1]; 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()
41 x1 = xx[ncnt + 1]; in PCApply_VPBJacobi()
44 yy[ncnt] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_VPBJacobi()
[all …]
|
| /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
24 x1 = xx[2 * i + 1]; in PCApply_PBJacobi()
25 yy[2 * i] = diag[0] * x0 + diag[2] * x1; in PCApply_PBJacobi()
26 yy[2 * i + 1] = diag[1] * x0 + diag[3] * x1; in PCApply_PBJacobi()
33 x1 = xx[3 * i + 1]; 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()
45 x1 = xx[4 * i + 1]; in PCApply_PBJacobi()
49 yy[4 * i] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in PCApply_PBJacobi()
[all …]
|
| /petsc/src/vec/vec/tests/ |
| 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()
40 PetscCall(VecCreate(PETSC_COMM_SELF, &X1)); in main()
41 PetscCall(VecSetSizes(X1, N / 2, N / 2)); in main()
42 PetscCall(VecSetFromOptions(X1)); in main()
43 PetscCall(VecDuplicate(X1, &X2)); in main()
44 PetscCall(VecDuplicate(X1, &Y1)); in main()
45 PetscCall(VecDuplicate(X1, &Y2)); in main()
46 PetscCall(VecDuplicate(X1, &F1)); in main()
[all …]
|
| H A D | ex32.c | 15 PetscInt *X, *X1, *XR, *XSO, *W, *Y, *Z, *XP, *X1P; in main() local
38 PetscCall(PetscCalloc6(n, &X, n, &X1, n, &XR, n, &XSO, n, &Y, n, &Z)); in main()
75 PetscCall(PetscArraycpy(X1, order ? XSO : XR, n)); in main()
79 PetscCall(PetscIntSortSemiOrdered(n, X1)); in main()
89 …X1[i], PETSC_COMM_SELF, PETSC_ERR_PLIB, "PetscIntSortSemiOrdered() rep %" PetscInt_FMT " X1[%" Pet… in main()
91 …for (i = 0; i < n - 1; i++) PetscCheck(X1[i] <= X1[i + 1], PETSC_COMM_SELF, PETSC_ERR_PLIB, "Petsc… in main()
93 PetscCall(PetscArrayzero(X1, n)); in main()
110 PetscCall(PetscArraycpy(X1, order ? XSO : XR, n)); in main()
118 PetscCall(PetscIntSortSemiOrderedWithArray(n, X1, X1P)); in main()
138 …X1[%" PetscInt_FMT "]:%" PetscInt_FMT " does not match PetscSortIntWithArray() X[%" PetscInt_FMT "… in main()
[all …]
|
| H A D | ex60.c | 9 Vec x, x1, x2, x3; in main() local
27 PetscCall(VecCreateSeqCUDAWithArray(PETSC_COMM_SELF, bs, n, NULL, &x1)); in main()
33 PetscCall(VecCreateSeqKokkosWithArray(PETSC_COMM_SELF, bs, n, NULL, &x1)); in main()
39 PetscCall(VecCreateSeqHIPWithArray(PETSC_COMM_SELF, bs, n, NULL, &x1)); in main()
44 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, bs, n, NULL, &x1)); in main()
50 PetscCall(VecPlaceArray(x1, px)); in main()
53 PetscCall(VecSet(x1, 1.0)); in main()
56 PetscCall(VecResetArray(x1)); in main()
60 PetscCall(VecDestroy(&x1)); in main()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij2.c | 286 PetscScalar x1, x2; in MatMult_SeqBAIJ_2() local
317 x1 = xb[0]; 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
374 x1 = xb[0]; 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
35 x1 = x[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()
66 x1 = x[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
37 x1 = x[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()
60 x1 = x[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
108 x1 = t[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()
148 x1 = t[idx]; in MatSolve_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvnat4.c | 31 PetscScalar s1, s2, s3, s4, x1, x2, x3, x4; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local
53 x1 = x[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()
81 x1 = x[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 | baijsolvnat2.c | 17 PetscScalar *x, s1, s2, x1, x2; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local
38 x1 = x[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()
57 x1 = x[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
101 x1 = x[jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering()
103 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_NaturalOrdering()
[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
37 x1 = x[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()
72 x1 = x[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
39 x1 = x[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()
78 x1 = x[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
20 x1 = x[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
94 x1 = x[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 …]
|
| /petsc/src/vec/vec/tutorials/ |
| H A D | ex19.c | 8 Vec x1, x2, *x3ts, *x4ts; in main() local
31 PetscCall(VecCreate(comm, &x1)); in main()
32 PetscCall(PetscObjectSetName((PetscObject)x1, "x1")); in main()
33 PetscCall(VecSetSizes(x1, PETSC_DECIDE, n)); in main()
34 PetscCall(VecSetFromOptions(x1)); in main()
35 PetscCall(VecSetRandom(x1, rand)); in main()
36 PetscCall(VecView(x1, viewer)); in main()
51 PetscCall(VecDuplicateVecs(x1, n_timesteps, &x3ts)); in main()
78 PetscCall(VecDuplicate(x1, &x1r)); in main()
79 PetscCall(PetscObjectSetName((PetscObject)x1r, "x1")); in main()
[all …]
|
| /petsc/src/snes/interface/ |
| H A D | snesj.c | 12 . x1 - compute Jacobian at this point
48 PetscErrorCode SNESComputeJacobianDefault(SNES snes, Vec x1, Mat J, Mat B, PetscCtx ctx) in SNESComputeJacobianDefault() argument
70 PetscCall(PetscObjectGetComm((PetscObject)x1, &comm)); in SNESComputeJacobianDefault()
81 PetscCall(VecDuplicateVecs(x1, snes->nvwork, &snes->vwork)); in SNESComputeJacobianDefault()
88 PetscCall(VecGetSize(x1, &N)); in SNESComputeJacobianDefault()
89 PetscCall(VecGetOwnershipRange(x1, &start, &end)); in SNESComputeJacobianDefault()
93 PetscCall(SNESComputeMFFunction(snes, x1, j1a)); in SNESComputeJacobianDefault()
95 …PetscCall(SNESComputeFunction(snes, x1, j1a)); /* does not handle use of SNESSetFunctionDomainErro… in SNESComputeJacobianDefault()
103 if (use_wp) PetscCall(VecNorm(x1, NORM_2, &unorm)); in SNESComputeJacobianDefault()
105 x1 = current iterate, j1a = F(x1) 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
110 x1 = x[1]; in VecMDot_Seq()
111 sum0 += x1 * PetscConj(yy0[1]); in VecMDot_Seq()
112 sum1 += x1 * PetscConj(yy1[1]); in VecMDot_Seq()
113 sum2 += x1 * PetscConj(yy2[1]); /* fall through */ in VecMDot_Seq()
129 x1 = x[1]; 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()
158 x1 = x[1]; in VecMDot_Seq()
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 290 PetscScalar *z, x1, x2, zero = 0.0; in MatMult_SeqSBAIJ_2() local
308 x1 = xb[0]; 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
362 x1 = xb[0]; 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
250 x1 = tp[1]; 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/sys/classes/draw/impls/x/ |
| H A D | xtone.c | 10 PetscErrorCode PetscDrawInterpolatedTriangle_X(PetscDraw_X *win, int x1, int y_1, int t1, int x2, i… in PetscDrawInterpolatedTriangle_X() argument
21 if (x1 < 0 && x2 < 0 && x3 < 0) PetscFunctionReturn(PETSC_SUCCESS); in PetscDrawInterpolatedTriangle_X()
23 if (x1 > win->w && x2 > win->w && x3 > win->w) PetscFunctionReturn(PETSC_SUCCESS); in PetscDrawInterpolatedTriangle_X()
41 SWAP(x1, x2); in PetscDrawInterpolatedTriangle_X()
46 SWAP(x1, x3); in PetscDrawInterpolatedTriangle_X()
61 x2_x1 = x2 - x1; in PetscDrawInterpolatedTriangle_X()
63 x3_x1 = x3 - x1; in PetscDrawInterpolatedTriangle_X()
69 lx = (int)(lfrac * ((PetscReal)x2_x1) + (PetscReal)x1); in PetscDrawInterpolatedTriangle_X()
73 rx = (int)(rfrac * ((PetscReal)x3_x1) + (PetscReal)x1); in PetscDrawInterpolatedTriangle_X()
102 x1 = rx; in PetscDrawInterpolatedTriangle_X()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex97.c | 51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument
81 PetscCall(VecCopy(X1, rtmp[0])); in CheckMatrices()
82 PetscCall(VecCopy(X1, rtmp[1])); in CheckMatrices()
87 PetscCall(MatMultTransposeAdd(A, Y, X1, rtmp[0])); in CheckMatrices()
88 PetscCall(MatMultTransposeAdd(B, Y, X1, rtmp[1])); in CheckMatrices()
100 Vec left, right, X, Y, X1, Y1; in main() local
128 PetscCall(VecDuplicate(right, &X1)); in main()
138 PetscCall(VecSetRandom(X1, NULL)); in main()
145 PetscCall(VecSet(X1, 3.0)); in main()
148 PetscCall(CheckMatrices(Asub, Bsub, left, right, X, Y, X1, Y1)); in main()
[all …]
|
| H A D | ex99.c | 51 static PetscErrorCode CheckMatrices(Mat A, Mat B, Vec left, Vec right, Vec X, Vec Y, Vec X1, Vec Y1) in CheckMatrices() argument
83 PetscCall(VecCopy(X1, rtmp[0])); in CheckMatrices()
84 PetscCall(VecCopy(X1, rtmp[1])); in CheckMatrices()
89 PetscCall(MatMultTransposeAdd(A, Y, X1, rtmp[0])); in CheckMatrices()
90 PetscCall(MatMultTransposeAdd(B, Y, X1, rtmp[1])); in CheckMatrices()
102 Vec left, right, X, Y, X1, Y1; in main() local
129 PetscCall(VecDuplicate(right, &X1)); in main()
139 PetscCall(VecSetRandom(X1, NULL)); in main()
146 PetscCall(VecSet(X1, 3.0)); in main()
149 PetscCall(CheckMatrices(Asub, Bsub, left, right, X, Y, X1, Y1)); in main()
[all …]
|
| /petsc/src/vec/is/sf/tests/ |
| H A D | ex11.c | 12 Vec x, x1, y; in main() local
46 PetscCall(VecDuplicate(x, &x1)); in main()
47 PetscCall(VecCopy(x, x1)); in main()
48 PetscCall(VecEqual(x, x1, &flg)); in main()
49 PetscCheck(flg, PetscObjectComm((PetscObject)x), PETSC_ERR_ARG_WRONG, "x1 != x"); in main()
51 PetscCall(VecScale(x1, 2.0)); in main()
52 PetscCall(VecSet(x1, 10.0)); in main()
53 /* PetscCall(VecView(x1,PETSC_VIEWER_STDOUT_WORLD)); */ in main()
112 PetscCall(VecDestroy(&x1)); in main()
|
| /petsc/src/vec/vec/tutorials/output/ |
| H A D | ex31_1.testout | 10 MPI_Comm_rank 1x1 8 double
11 MPI_Comm_size 1x1 8 double
12 ans 1x1 8 matlab.settings.SettingsGroup
13 x 5x1 40 double
|