| /petsc/src/dm/field/impls/ds/ |
| H A D | dmfieldds.c | 151 PetscInt nq, nc, dim, meshDim, numCells, feDim, i, K = H ? 2 : (D ? 1 : (B ? 0 : -1)); in DMFieldEvaluateFE_DS() local
163 PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &nq, &qpoints, NULL)); in DMFieldEvaluateFE_DS()
177 PetscCall(PetscFECreateTabulation(fe, 1, nq, qpoints, K, &T)); in DMFieldEvaluateFE_DS()
189 PetscScalar *cB = &((PetscScalar *)B)[nc * nq * i]; in DMFieldEvaluateFE_DS()
191 DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, (PetscScalar)); in DMFieldEvaluateFE_DS()
193 PetscReal *cB = &((PetscReal *)B)[nc * nq * i]; in DMFieldEvaluateFE_DS()
195 DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, PetscRealPart); in DMFieldEvaluateFE_DS()
200 PetscScalar *cD = &((PetscScalar *)D)[nc * nq * dim * i]; in DMFieldEvaluateFE_DS()
202 DMFieldDSdot(cD, T->T[1], elem, nq, feDim, (nc * dim), (PetscScalar)); in DMFieldEvaluateFE_DS()
204 PetscReal *cD = &((PetscReal *)D)[nc * nq * dim * i]; in DMFieldEvaluateFE_DS()
[all …]
|
| /petsc/src/dm/field/impls/da/ |
| H A D | dmfieldda.c | 216 PetscInt nq, nc; in DMFieldEvaluateFE_DA() local
245 PetscCall(PetscQuadratureGetData(points, NULL, NULL, &nq, &q, NULL)); in DMFieldEvaluateFE_DA()
247 PetscCall(DMGetWorkArray(dm, nq * dim, MPIU_SCALAR, &qs)); in DMFieldEvaluateFE_DA()
248 for (i = 0; i < nq * dim; i++) qs[i] = q[i]; in DMFieldEvaluateFE_DA()
267 cB = PetscSafePointerPlusOffset((PetscScalar *)B, nc * nq * c); in DMFieldEvaluateFE_DA()
268 cD = PetscSafePointerPlusOffset((PetscScalar *)D, nc * nq * dim * c); in DMFieldEvaluateFE_DA()
269 cH = PetscSafePointerPlusOffset((PetscScalar *)H, nc * nq * dim * dim * c); in DMFieldEvaluateFE_DA()
271 cB = PetscSafePointerPlusOffset((PetscReal *)B, nc * nq * c); in DMFieldEvaluateFE_DA()
272 cD = PetscSafePointerPlusOffset((PetscReal *)D, nc * nq * dim * c); in DMFieldEvaluateFE_DA()
273 cH = PetscSafePointerPlusOffset((PetscReal *)H, nc * nq * dim * dim * c); in DMFieldEvaluateFE_DA()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex10.c | 506 static PetscErrorCode RDGetQuadrature(RD rd, PetscReal hx, PetscInt *nq, PetscReal weight[], PetscR… in RDGetQuadrature() argument
521 *nq = 1; in RDGetQuadrature()
536 *nq = 2; in RDGetQuadrature()
553 *nq = 3; in RDGetQuadrature()
573 *nq = 4; in RDGetQuadrature()
588 *nq = 2; in RDGetQuadrature()
605 *nq = 3; in RDGetQuadrature()
614 for (q = 0; q < *nq; q++) { in RDGetQuadrature()
634 PetscInt i, j, q, nq; in RDIFunction_FE() local
646 PetscCall(RDGetQuadrature(rd, hx, &nq, weight, interp, deriv)); in RDIFunction_FE()
[all …]
|
| /petsc/src/dm/impls/forest/p4est/ |
| H A D | pforest.h | 3260 p4est_quadrant_t nq; in DMPforestLabelsInitialize() local
3264 PetscCallP4est(p4est_quadrant_face_neighbor, (q, l, &nq)); in DMPforestLabelsInitialize()
3265 PetscCallP4estReturn(isInside, p4est_quadrant_is_inside_root, (&nq)); in DMPforestLabelsInitialize()
3286 p4est_quadrant_t nq; in DMPforestLabelsInitialize() local
3290 PetscCallP4est(p8est_quadrant_edge_neighbor, (q, l, &nq)); in DMPforestLabelsInitialize()
3291 PetscCallP4estReturn(isInside, p4est_quadrant_is_inside_root, (&nq)); in DMPforestLabelsInitialize()
3303 PetscCallP4estReturn(isOutsideFace, p4est_quadrant_is_outside_face, (&nq)); in DMPforestLabelsInitialize()
3307 if (nq.x < 0) { in DMPforestLabelsInitialize()
3309 } else if (nq.x >= P4EST_ROOT_LEN) { in DMPforestLabelsInitialize()
3311 } else if (nq.y < 0) { in DMPforestLabelsInitialize()
[all …]
|
| /petsc/src/dm/field/tutorials/ |
| H A D | ex1.c | 64 PetscInt dim, i, nc, nq; in TestEvaluateFE() local
88 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &nq, NULL, NULL)); in TestEvaluateFE()
89 N = n * nq * nc; in TestEvaluateFE()
|