Home
last modified time | relevance | path

Searched refs:d (Results 1 – 25 of 517) sorted by relevance

12345678910>>...21

/petsc/share/petsc/datafiles/meshes/
H A Dgmsh-3d-binary.msh248 �����?P��<���?��������d���^������?^������?��������e���A������?������?��������f���������?������?…
62 ���g������������Z������2���d��� ���g������������[������2��� ���d���j������������\������…
65 …d������2������b���m������������e������2���������h������������f������2������d���h����…
H A Dwedge_cylinder.cas26 4 2b 1 2d 2c 1 8
27 4 1 2 2e 2d 2 9
32 4 5 6 32 31 6 d
33 4 1 2d 37 7 8 9
44 4 5 31 3b b c d
45 4 5 6 c b d 37
47 4 6 32 3c c d e
51 4 35 7 d 3f f 3a
52 4 b c 3c 3b d 14
53 4 7 37 41 d f 10
[all …]
H A Dsquare.cas111 6 3d 0 50
114 4 3e 1d 0
115 3e d 1f 0
116 d 3b 1b 0
126 20 21 d 0
141 51 3d 4d 0
148 4 3f 1e 1d
149 3e 3f 1d 22
151 d c 1f 5
152 d e 5 1b
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact12.c13 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
59 d = ba + i * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
112 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
113 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
114 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
115 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
116 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
117 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
118 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
120 …uik[7] = -(d[0] * u7 + d[7] * u8 + d[14] * u9 + d[21] * u10 + d[28] * u11 + d[35] * u12 + d[42] *… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
H A Dsbaijfact11.c12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
96 d = ba + i * 49; in MatCholeskyFactorNumeric_SeqSBAIJ_7()
149 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
150 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
151 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
152 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
153 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
154 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
155 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
157 …uik[7] = -(d[0] * u7 + d[7] * u8 + d[14] * u9 + d[21] * u10 + d[28] * u11 + d[35] * u12 + d[42] *… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
[all …]
H A Dsbaijfact8.c13 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
55 d = ba + i * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
58 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
59 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
60 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
61 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
62 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
64 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
65 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
66 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
H A Dsbaijfact9.c12 MatScalar *u, *d, *w, *wp, u0, u1, u2, u3, u4, u5, u6, u7, u8, u9, u10, u11, u12; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
95 d = ba + i * 36; in MatCholeskyFactorNumeric_SeqSBAIJ_6()
135 uik[0] = -(d[0] * u0 + d[6] * u1 + d[12] * u2 + d[18] * u3 + d[24] * u4 + d[30] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
136 uik[1] = -(d[1] * u0 + d[7] * u1 + d[13] * u2 + d[19] * u3 + d[25] * u4 + d[31] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
137 uik[2] = -(d[2] * u0 + d[8] * u1 + d[14] * u2 + d[20] * u3 + d[26] * u4 + d[32] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
138 uik[3] = -(d[3] * u0 + d[9] * u1 + d[15] * u2 + d[21] * u3 + d[27] * u4 + d[33] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
139 uik[4] = -(d[4] * u0 + d[10] * u1 + d[16] * u2 + d[22] * u3 + d[28] * u4 + d[34] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
140 uik[5] = -(d[5] * u0 + d[11] * u1 + d[17] * u2 + d[23] * u3 + d[29] * u4 + d[35] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
142 uik[6] = -(d[0] * u6 + d[6] * u7 + d[12] * u8 + d[18] * u9 + d[24] * u10 + d[30] * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
143 uik[7] = -(d[1] * u6 + d[7] * u7 + d[13] * u8 + d[19] * u9 + d[25] * u10 + d[31] * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
[all …]
H A Dsbaijfact7.c12 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
93 d = ba + i * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
96 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
97 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
98 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
99 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
100 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
102 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
103 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
104 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
[all …]
/petsc/src/dm/interface/
H A Ddmperiodicity.c49 PetscInt dim, d; in DMSetPeriodicity() local
59 for (d = 0; d < dim; ++d) dm->maxCell[d] = maxCell[d]; in DMSetPeriodicity()
66 for (d = 0; d < dim; ++d) dm->Lstart[d] = Lstart[d]; in DMSetPeriodicity()
73 for (d = 0; d < dim; ++d) dm->L[d] = L[d]; in DMSetPeriodicity()
99 PetscInt dim, d; in DMLocalizeCoordinate() local
104 for (d = 0; d < dim; ++d) out[d] = in[d]; in DMLocalizeCoordinate()
107 for (d = 0; d < dim; ++d) { in DMLocalizeCoordinate()
108 …l(PetscRealPart(in[d]) / dm->L[d] - PetscFloorReal(PetscRealPart(in[d]) / dm->L[d])) < PETSC_SMALL… in DMLocalizeCoordinate()
109 out[d] = in[d] - dm->L[d] * (PetscFloorReal(PetscRealPart(in[d]) / dm->L[d]) - 1); in DMLocalizeCoordinate()
111 out[d] = in[d] - dm->L[d] * PetscFloorReal(PetscRealPart(in[d]) / dm->L[d]); in DMLocalizeCoordinate()
[all …]
H A Ddmgeommodel.c19 for (PetscInt d = 0; d < dE; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelBall() local
21 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d] / norm; in DMSnapToGeomModelBall() local
23 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelBall() local
39 for (PetscInt d = 0; d < dE - 1; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelCylinder() local
41 for (PetscInt d = 0; d < dE - 1; ++d) gcoords[d] = mcoords[d] * gmax[0] / norm; in DMSnapToGeomModelCylinder() local
44 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelCylinder() local
183 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModel() local
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Ddfovec.m45 v = [4.0d0,2.0d0,1.0d0,5.0d-1,2.5d-1,1.67d-1,1.25d-1,1.0d-1,8.33d-2,...
46 7.14d-2,6.25d-2];
47 y1 = [1.4d-1,1.8d-1,2.2d-1,2.5d-1,2.9d-1,3.2d-1,3.5d-1,3.9d-1,3.7d-1,...
48 5.8d-1,7.3d-1,9.6d-1,1.34d0,2.1d0,4.39d0];
49 y2 = [1.957d-1,1.947d-1,1.735d-1,1.6d-1,8.44d-2,6.27d-2,4.56d-2,3.42d-2,...
50 3.23d-2,2.35d-2,2.46d-2];
53 y4 = [8.44d-1,9.08d-1,9.32d-1,9.36d-1,9.25d-1,9.08d-1,8.81d-1,8.5d-1,...
54 8.18d-1,7.84d-1,7.51d-1,7.18d-1,6.85d-1,6.58d-1,6.28d-1,6.03d-1,...
55 5.8d-1,5.58d-1,5.38d-1,5.22d-1,5.06d-1,4.9d-1,4.78d-1,4.67d-1,...
56 4.57d-1,4.48d-1,4.38d-1,4.31d-1,4.24d-1,4.2d-1,4.14d-1,4.11d-1,...
[all …]
/petsc/include/petsc/private/
H A Dpetsclegacycupmblas.h9 …#define cublasXaxpy(a, b, c, d, e, f, g) cublasCaxpy((a), (b), (cuComplex *)(… argument
10 …#define cublasXscal(a, b, c, d, e) cublasCscal((a), (b), (cuComplex *)(… argument
11 …#define cublasXdotu(a, b, c, d, e, f, g) cublasCdotu((a), (b), (cuComplex *)(… argument
12 …#define cublasXdot(a, b, c, d, e, f, g) cublasCdotc((a), (b), (cuComplex *)(… argument
13 …#define cublasXswap(a, b, c, d, e, f) cublasCswap((a), (b), (cuComplex *)(… argument
14 …#define cublasXnrm2(a, b, c, d, e) cublasScnrm2((a), (b), (cuComplex *)… argument
15 …#define cublasIXamax(a, b, c, d, e) cublasIcamax((a), (b), (cuComplex *)… argument
16 …#define cublasXasum(a, b, c, d, e) cublasScasum((a), (b), (cuComplex *)… argument
17 …#define cublasXgemv(a, b, c, d, e, f, g, h, i, j, k, l) cublasCgemv((a), (b), (c), (d), (cuC… argument
18 …#define cublasXgemm(a, b, c, d, e, f, g, h, i, j, k, l, m, n) cublasCgemm((a), (b), (c), (d), (e),… argument
[all …]
/petsc/src/dm/impls/stag/
H A Dstag2d.c60 PetscInt xc, yc, mc, nc, nExtraxc, nExtrayc, i, j, d; in DMStagRestrictSimple_2d() local
85 for (d = 0; d < dof[0]; ++d) in DMStagRestrictSimple_2d()
90 arrc[j][i][idownleftc + d] = arrf[jj][ii][idownleftf + d]; in DMStagRestrictSimple_2d()
93 for (d = 0; d < dof[1]; ++d) in DMStagRestrictSimple_2d()
98 …if (factory % 2 == 0) arrc[j][i][ileftc + d] = 0.5 * (arrf[jj - 1][ii][ileftf + d] + arrf[jj][ii][… in DMStagRestrictSimple_2d()
99 else arrc[j][i][ileftc + d] = arrf[jj][ii][ileftf + d]; in DMStagRestrictSimple_2d()
102 for (d = 0; d < dof[1]; ++d) in DMStagRestrictSimple_2d()
107 …if (factorx % 2 == 0) arrc[j][i][idownc + d] = 0.5 * (arrf[jj][ii - 1][idownf + d] + arrf[jj][ii][… in DMStagRestrictSimple_2d()
108 else arrc[j][i][idownc + d] = arrf[jj][ii][idownf + d]; in DMStagRestrictSimple_2d()
111 for (d = 0; d < dof[2]; ++d) in DMStagRestrictSimple_2d()
[all …]
H A Dstag3d.c67 PetscInt xc, yc, zc, mc, nc, pc, nExtraxc, nExtrayc, nExtrazc, i, j, k, d; in DMStagRestrictSimple_3d() local
102 for (d = 0; d < dof[0]; ++d) in DMStagRestrictSimple_3d()
108 arrc[k][j][i][ibackdownleftc + d] = arrf[kk][jj][ii][ibackdownleftf + d]; in DMStagRestrictSimple_3d()
111 for (d = 0; d < dof[1]; ++d) in DMStagRestrictSimple_3d()
117 …2 == 0) arrc[k][j][i][ibackdownc + d] = 0.5 * (arrf[kk][jj][ii - 1][ibackdownf + d] + arrf[kk][jj]… in DMStagRestrictSimple_3d()
118 else arrc[k][j][i][ibackdownc + d] = arrf[kk][jj][ii][ibackdownf + d]; in DMStagRestrictSimple_3d()
121 for (d = 0; d < dof[1]; ++d) in DMStagRestrictSimple_3d()
127 …2 == 0) arrc[k][j][i][ibackleftc + d] = 0.5 * (arrf[kk][jj - 1][ii][ibackleftf + d] + arrf[kk][jj]… in DMStagRestrictSimple_3d()
128 else arrc[k][j][i][ibackleftc + d] = arrf[kk][jj][ii][ibackleftf + d]; in DMStagRestrictSimple_3d()
131 for (d = 0; d < dof[1]; ++d) in DMStagRestrictSimple_3d()
[all …]
H A Dstag.c12 PetscInt f0, f1, f2, f3, dof0, dof1, dof2, dof3, n_entries, k, d, cnt, n_fields, dim; in DMCreateFieldDecomposition_Stag() local
38 for (d = 0; d < dof0; ++d) { in DMCreateFieldDecomposition_Stag()
39 stencil0[dof0 * k + d].i = 0; in DMCreateFieldDecomposition_Stag()
40 stencil0[dof0 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
41 stencil0[dof0 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
45 for (d = 0; d < dof1; ++d) { in DMCreateFieldDecomposition_Stag()
46 stencil1[dof1 * k + d].i = 0; in DMCreateFieldDecomposition_Stag()
47 stencil1[dof1 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
48 stencil1[dof1 * k + d].j = 0; in DMCreateFieldDecomposition_Stag()
53 for (d = 0; d < dof2; ++d) { in DMCreateFieldDecomposition_Stag()
[all …]
/petsc/src/dm/impls/plex/tutorials/
H A Dex20.c59 for (PetscInt d = 0; d < dim; d++) { in main() local
66 PetscCall(DMCreateGlobalVector(dm, &Vecxyz[d])); in main()
67 …PetscCall(PetscObjectSetName((PetscObject)Vecxyz[d], d == 0 ? "x_coordinate" : (d == 1 ? "y_coordi… in main()
81 xArray[sOff] = coordArray[cOff + d]; in main()
86 PetscCall(DMLocalToGlobal(dm, loc_vec, INSERT_VALUES, Vecxyz[d])); in main()
88 PetscCall(VecGetSize(Vecxyz[d], &n_state_global)); in main()
110 for (PetscInt d = 0; d < dim; d++) { in main() local
111 PetscCall(DMCreateGlobalVector(dm, &Vecxyz[d])); in main()
112 PetscCall(VecSetFromOptions(Vecxyz[d])); in main()
113 …PetscCall(PetscObjectSetName((PetscObject)Vecxyz[d], d == 0 ? "x_coordinate" : (d == 1 ? "y_coordi… in main()
[all …]
/petsc/src/dm/impls/plex/tests/
H A Dex8.c111 PetscInt vStart, vEnd, v, d, coordSize; in ChangeCoordinates() local
134 for (d = 0; d < spaceDim; ++d) coords[off + d] = vertexCoords[(v - vStart) * spaceDim + d]; in ChangeCoordinates()
149 PetscInt d, i, j; in CheckFEMGeometry() local
153 for (d = 0; d < spaceDim; ++d) { in CheckFEMGeometry()
154 if (v0[d] != v0Ex[d]) { in CheckFEMGeometry()
179 PetscInt d; in CheckFVMGeometry() local
183 for (d = 0; d < spaceDim; ++d) { in CheckFVMGeometry()
185 …d], centroidEx[d]) < tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid centro… in CheckFVMGeometry()
186 …d], normalEx[d]) < tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid normal[%… in CheckFVMGeometry()
198 PetscInt coneSize, f, dim, cdim, d; in CheckGaussLaw() local
[all …]
H A Dex100.c14 for (PetscInt d = 0; d < dim; d++) x_tot += x[d]; in project_function() local
88 for (PetscInt d = 0; d < cdim; ++d) { in main() local
89 min[d] = PetscMin(min[d], PetscRealPart(cell_coords[n * cdim + d])); in main()
90 max[d] = PetscMax(max[d], PetscRealPart(cell_coords[n * cdim + d])); in main()
94 for (PetscInt d = 0; d < cdim; d++) cell_bounding_box_size[d] = max[d] - min[d]; in main() local
97 for (PetscInt d = 0; d < cdim; d++) { in main() local
98 …e = PetscAbsReal((ref_cell_bounding_box_size[d] - cell_bounding_box_size[d]) / ref_cell_bounding_b… in main()
100 … in dimension %" PetscInt_FMT ". Relative difference = %g\n", cell, cStart, d, (double)difference); in main()
/petsc/src/dm/impls/swarm/
H A Ddata_ex.c87 DMSwarmDataEx d; in DMSwarmDataExCreate() local
90 PetscCall(PetscNew(&d)); in DMSwarmDataExCreate()
91 PetscCallMPI(MPI_Comm_dup(comm, &d->comm)); in DMSwarmDataExCreate()
92 PetscCallMPI(MPI_Comm_rank(d->comm, &d->rank)); in DMSwarmDataExCreate()
94 d->instance = count; in DMSwarmDataExCreate()
96 d->topology_status = DEOBJECT_STATE_UNKNOWN; in DMSwarmDataExCreate()
97 d->message_lengths_status = DEOBJECT_STATE_UNKNOWN; in DMSwarmDataExCreate()
98 d->packer_status = DEOBJECT_STATE_UNKNOWN; in DMSwarmDataExCreate()
99 d->communication_status = DEOBJECT_STATE_UNKNOWN; in DMSwarmDataExCreate()
101 d->n_neighbour_procs = -1; in DMSwarmDataExCreate()
[all …]
/petsc/src/dm/impls/product/
H A Dproduct.c6 PetscInt d; in DMDestroy_Product() local
9 for (d = 0; d < DMPRODUCT_MAX_DIM; ++d) PetscCall(DMDestroy(&product->dm[d])); in DMDestroy_Product()
17 PetscInt d; in DMView_Product() local
20 for (d = 0; d < DMPRODUCT_MAX_DIM; ++d) { in DMView_Product()
21 if (product->dm[d]) { in DMView_Product()
22 … PetscCall(PetscViewerASCIIPrintf(viewer, " DM that defines dimension %" PetscInt_FMT "\n", d)); in DMView_Product()
25 PetscCall(DMView(product->dm[d], viewer)); in DMView_Product()
66 PetscInt d; in DMCreate_Product() local
73 for (d = 0; d < DMPRODUCT_MAX_DIM; ++d) product->dm[d] = NULL; in DMCreate_Product()
74 for (d = 0; d < DMPRODUCT_MAX_DIM; ++d) product->dim[d] = -1; in DMCreate_Product()
/petsc/src/dm/impls/stag/tests/
H A Dex7.c10 …PetscInt startx, starty, startz, nx, ny, nz, i, j, k, d, is, js, ks, dof0, dof1, dof2, dof3,… in main() local
44 for (d = 0; d < dofTotal; ++d) { in main()
45 …][j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expecting %g)… in main()
46 a2[k][j][i][d] = 0.0; in main()
49 …for (is = -stencilWidth; is <= stencilWidth; ++is) a2[k][j][i][d] += a1[k + ks][j + js][i + is][d]; in main()
67 for (d = 0; d < 3; ++d) expected *= (2 * stencilWidth + 1); in main()
90 for (d = 0; d < dof0; ++d) { in main()
91 if (a2[k][j][i][d] != expected) { in main()
92 …"] Unexpected value %g (expecting %g)\n", rank, i, j, k, d, (double)PetscRealPart(a2[k][j][i][d]),… in main()
99 for (d = dof0; d < dof0 + dof1; ++d) { in main()
[all …]
H A Dex10.c10 …PetscInt startx, starty, nx, ny, i, j, d, is, js, dof0, dof1, dof2, dofTotal, stencilWidth, … in main() local
42 for (d = 0; d < dofTotal; ++d) { in main()
43 …(a1[j][i][d] != 1.0) PetscCall(PetscPrintf(PETSC_COMM_SELF, "[%d] Unexpected value %g (expecting %… in main()
44 a2[j][i][d] = 0.0; in main()
46 for (is = -stencilWidth; is <= stencilWidth; ++is) a2[j][i][d] += a1[j + js][i + is][d]; in main()
62 for (d = 0; d < 2; ++d) expected *= (2 * stencilWidth + 1); in main()
82 for (d = 0; d < dof0; ++d) { in main()
83 if (a2[j][i][d] != expected) { in main()
84 …_FMT "] Unexpected value %g (expecting %g)\n", rank, i, j, d, (double)PetscRealPart(a2[j][i][d]), … in main()
91 for (d = dof0; d < dof0 + dof1; ++d) { in main()
[all …]
/petsc/src/dm/dt/tests/
H A Dex8.c9 PetscInt d, n, maxdim = 4; in main() local
16 for (d = 0; d <= maxdim; d++) { in main()
17 for (n = 0; n <= d + 2; n++) { in main()
20 PetscCall(PetscDTBinomialInt(d + n, d, &Nk)); in main()
24 PetscCall(PetscDTIndexToBary(d + 1, n, k, btup)); in main()
25 for (j = 0, sum = 0; j < d + 1; j++) { in main()
26 …ary, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FMT " negative entry", d, n, k); in main()
29 …Bary, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FMT " incorrect sum", d, n, k); in main()
30 PetscCall(PetscDTBaryToIndex(d + 1, n, btup, &kchk)); in main()
31 …ryToIndex, d = %" PetscInt_FMT ", n = %" PetscInt_FMT ", k = %" PetscInt_FMT " mismatch", d, n, k); in main()
[all …]
/petsc/src/snes/tests/
H A Dex8.c152 PetscInt d = user->dir; in constant() local
155 for (d = 0; d < Nc; ++d) u[d] = user->constants[d]; in constant()
157 u[0] = user->constants[d]; in constant()
164 PetscInt d = user->dir; in constantDer() local
167 for (d = 0; d < Nc; ++d) u[d] = 0.0; in constantDer()
169 u[0] = user->constants[d]; in constantDer()
178 PetscInt d = user->dir; in linear() local
181 for (d = 0; d < Nc; ++d) u[d] = coords[d]; in linear()
183 u[0] = coords[d]; in linear()
190 PetscInt d = user->dir; in linearDer() local
[all …]
/petsc/src/vec/is/sf/utils/
H A Dsfcoord.c6 for (PetscInt d = 0; d < dim; d++) { in GetBoundingBox_Internal() local
7 bbox[0 * dim + d] = PETSC_MAX_REAL; in GetBoundingBox_Internal()
8 bbox[1 * dim + d] = PETSC_MIN_REAL; in GetBoundingBox_Internal()
11 for (PetscInt d = 0; d < dim; d++) { in GetBoundingBox_Internal() local
12 bbox[0 * dim + d] = PetscMin(bbox[0 * dim + d], coords[i * dim + d]); in GetBoundingBox_Internal()
13 bbox[1 * dim + d] = PetscMax(bbox[1 * dim + d], coords[i * dim + d]); in GetBoundingBox_Internal()
21 for (PetscInt d = 0; d < dim; d++) { in IntersectBoundingBox_Internal() local
22 …if (a[1 * dim + d] + tol < b[0 * dim + d] || b[1 * dim + d] + tol < a[0 * dim + d]) return PETSC_F… in IntersectBoundingBox_Internal()
29 for (PetscInt d = 0; d < dim; d++) { in InBoundingBox_Internal() local
30 if (x[d] + tol < bbox[0 * dim + d] || bbox[1 * dim + d] + tol < x[d]) return PETSC_FALSE; in InBoundingBox_Internal()
[all …]

12345678910>>...21