Home
last modified time | relevance | path

Searched refs:quad (Results 1 – 25 of 60) sorted by relevance

123

/petsc/src/binding/petsc4py/demo/legacy/kspsolve/
H A Dpetsc-cg.tex4 & \quad i \Leftarrow 0 \\
5 & \quad r \Leftarrow b - A x \\
6 & \quad d \Leftarrow r \\
7 & \quad \delta_{0} \Leftarrow r^T r \\
8 & \quad \delta_{ } \Leftarrow \delta_{0} \\
9 & \quad \text{while}\;\; i < i_{max} \text{ and } \\
10 & \quad\quad\qquad \delta_{ } > \delta_{0} \epsilon^2 \text{ do} :\\
11 & \quad\quad\quad q \Leftarrow Ad \\
12 & \quad\quad\quad \alpha \Leftarrow \frac{\delta_{ }}{d^T q} \\
13 & \quad\quad\quad x \Leftarrow x + \alpha d\\
[all …]
/petsc/src/dm/dt/space/impls/point/
H A Dspacepoint.c14 PetscCall(PetscQuadratureView(pt->quad, viewer)); in PetscSpacePointView_Ascii()
16 …space in dimension %" PetscInt_FMT " on %" PetscInt_FMT " points\n", sp->Nv, pt->quad->numPoints)); in PetscSpacePointView_Ascii()
37 if (!pt->quad->points && sp->degree >= 0) { in PetscSpaceSetUp_Point()
38 PetscCall(PetscQuadratureDestroy(&pt->quad)); in PetscSpaceSetUp_Point()
39 …PetscDTStroudConicalQuadrature(sp->Nv, sp->Nc, PetscMax(sp->degree + 1, 1), -1.0, 1.0, &pt->quad)); in PetscSpaceSetUp_Point()
49 PetscCall(PetscQuadratureDestroy(&pt->quad)); in PetscSpaceDestroy_Point()
59 *dim = pt->quad->numPoints; in PetscSpaceGetDimension_Point()
66 PetscInt dim = sp->Nv, pdim = pt->quad->numPoints, d, p, i, c; in PetscSpaceEvaluate_Point()
69 …>quad->numPoints, PETSC_COMM_SELF, PETSC_ERR_SUP, "Cannot evaluate Point space on %" PetscInt_FMT … in PetscSpaceEvaluate_Point()
74 if (PetscAbsReal(points[p * dim + d] - pt->quad->points[p * dim + d]) > 1.0e-10) break; in PetscSpaceEvaluate_Point()
[all …]
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DDT.pyx6 self.obj = <PetscObject*> &self.quad
7 self.quad = NULL
26 CHKERR(PetscQuadratureView(self.quad, vwr))
46 CHKERR(PetscCLEAR(self.obj)); self.quad = newquad
60 CHKERR(PetscQuadratureDuplicate(self.quad, &newquad.quad))
73 CHKERR(PetscQuadratureDestroy(&self.quad))
98 CHKERR(PetscQuadratureGetData(self.quad, &cdim, &cnc, &cnpoints, &cpoints, &cweights))
115 CHKERR(PetscQuadratureGetNumComponents(self.quad, &cnc))
134 CHKERR(PetscQuadratureSetNumComponents(self.quad, cnc))
147 CHKERR(PetscQuadratureGetOrder(self.quad, &corder))
[all …]
H A DFE.pyx222 cdef Quad quad = Quad()
223 CHKERR(PetscFEGetQuadrature(self.fe, &quad.quad))
224 CHKERR(PetscINCREF(quad.obj))
225 return quad
372 cdef Quad quad = Quad()
373 CHKERR(PetscFEGetFaceQuadrature(self.fe, &quad.quad))
374 CHKERR(PetscINCREF(quad.obj))
375 return quad
377 def setQuadrature(self, Quad quad) -> Self: argument
384 quad
[all …]
H A DSpace.pyx493 def setPointPoints(self, Quad quad) -> None: argument
494 """Set the evaluation points for the space to be based on a quad.
503 quad
511 CHKERR(PetscSpacePointSetPoints(self.space, quad.quad))
514 """Return the evaluation points for the space as the points of a quad.
523 cdef Quad quad = Quad()
524 CHKERR(PetscSpacePointGetPoints(self.space, &quad.quad))
525 CHKERR(PetscINCREF(quad.obj))
526 return quad
843 CHKERR(PetscDualSpaceGetFunctional(self.dualspace, ci, &functional.quad))
[all …]
/petsc/src/dm/field/interface/
H A Ddmfield.c378 PetscErrorCode DMFieldCreateDefaultQuadrature(DMField field, IS pointIS, PetscQuadrature *quad) in DMFieldCreateDefaultQuadrature() argument
383 PetscAssertPointer(quad, 3); in DMFieldCreateDefaultQuadrature()
385 *quad = NULL; in DMFieldCreateDefaultQuadrature()
386 PetscTryTypeMethod(field, createDefaultQuadrature, pointIS, quad); in DMFieldCreateDefaultQuadrature()
406 PetscErrorCode DMFieldCreateDefaultFaceQuadrature(DMField field, IS pointIS, PetscQuadrature *quad) in DMFieldCreateDefaultFaceQuadrature() argument
411 PetscAssertPointer(quad, 3); in DMFieldCreateDefaultFaceQuadrature()
413 *quad = NULL; in DMFieldCreateDefaultFaceQuadrature()
414 PetscTryTypeMethod(field, createDefaultFaceQuadrature, pointIS, quad); in DMFieldCreateDefaultFaceQuadrature()
439 PetscErrorCode DMFieldCreateFEGeom(DMField field, IS pointIS, PetscQuadrature quad, PetscFEGeomMode… in DMFieldCreateFEGeom() argument
450 PetscValidHeader(quad, 3); in DMFieldCreateFEGeom()
[all …]
/petsc/src/dm/dt/tests/
H A Dex15.c16 PetscQuadrature quad; in testQuadrature() local
19 PetscCall(PetscDTSimplexQuadrature(dim, degree, type, &quad)); in testQuadrature()
20 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &num_points, &points, &weights)); in testQuadrature()
33 …PetscCall(PetscInfo(quad, "Dimension %" PetscInt_FMT ", degree %" PetscInt_FMT ", method %s, error… in testQuadrature()
38 PetscCall(PetscQuadratureDestroy(&quad)); in testQuadrature()
H A Dex17.c10 PetscQuadrature quad; in main() local
39 PetscCall(PetscDualSpaceGetAllData(joint, &quad, &mat)); in main()
40 PetscCall(PetscQuadratureView(quad, PETSC_VIEWER_STDOUT_SELF)); in main()
H A Dex16.c8 PetscQuadrature quad; in PetscFEVectorTest() local
19 PetscCall(PetscFEGetQuadrature(orig_fe, &quad)); in PetscFEVectorTest()
20 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &num_points, &points, NULL)); in PetscFEVectorTest()
/petsc/src/dm/field/impls/shell/
H A Ddmfieldshell.c31 PetscErrorCode DMFieldShellEvaluateFEDefault(DMField field, IS pointIS, PetscQuadrature quad, Petsc… in DMFieldShellEvaluateFEDefault() argument
43 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldShellEvaluateFEDefault()
45 PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &Nq, NULL, NULL)); in DMFieldShellEvaluateFEDefault()
152 PetscQuadrature quad; in DMFieldShellEvaluateFVDefault() local
160 PetscCall(DMFieldGetFVQuadrature_Internal(coordField, pointIS, &quad)); in DMFieldShellEvaluateFVDefault()
161 …PetscCheck(quad, comm, PETSC_ERR_ARG_WRONGSTATE, "coordinate field must have default quadrature fo… in DMFieldShellEvaluateFVDefault()
162 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, NULL, NULL)); in DMFieldShellEvaluateFVDefault()
164 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldShellEvaluateFVDefault()
170 PetscCall(PetscQuadratureDestroy(&quad)); in DMFieldShellEvaluateFVDefault()
/petsc/src/dm/dt/space/impls/ptrimmed/tests/
H A Dex1.c13 PetscQuadrature quad; in test() local
39 PetscCall(PetscDTStroudConicalQuadrature(dim, 1, maxD + 1, -1., 1., &quad)); in test()
40 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &npoints, &points, NULL)); in test()
51 PetscCall(PetscQuadratureDestroy(&quad)); in test()
/petsc/src/dm/field/tutorials/
H A Dex1.c61 …(DMField field, PetscInt n, PetscInt cStart, PetscInt cEnd, PetscQuadrature quad, PetscRandom rand) in TestEvaluateFE() argument
88 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &nq, NULL, NULL)); in TestEvaluateFE()
91 PetscCall(DMFieldEvaluateFE(field, cellIS, quad, PETSC_SCALAR, B, D, H)); in TestEvaluateFE()
92 PetscCall(DMFieldEvaluateFE(field, cellIS, quad, PETSC_REAL, rB, rD, rH)); in TestEvaluateFE()
95 PetscCall(PetscObjectSetName((PetscObject)quad, "Test quadrature")); in TestEvaluateFE()
96 PetscCall(PetscQuadratureView(quad, viewer)); in TestEvaluateFE()
230 PetscQuadrature quad = NULL; in main() local
275 if (simplex) PetscCall(PetscDTStroudConicalQuadrature(dim, 1, pointsPerEdge, -1.0, 1.0, &quad)); in main()
276 else PetscCall(PetscDTGaussTensorQuadrature(dim, 1, pointsPerEdge, -1.0, 1.0, &quad)); in main()
336 PetscCall(PetscDTGaussTensorQuadrature(dim, 1, pointsPerEdge, -1.0, 1.0, &quad)); in main()
[all …]
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c145 static PetscErrorCode DMFieldEvaluateFE_DS(DMField field, IS pointIS, PetscQuadrature quad, PetscDa… in DMFieldEvaluateFE_DS() argument
163 PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &nq, &qpoints, NULL)); in DMFieldEvaluateFE_DS()
286 PetscQuadrature quad; in DMFieldEvaluate_DS() local
296 PetscCall(PetscQuadratureCreate(PETSC_COMM_SELF, &quad)); in DMFieldEvaluate_DS()
299 PetscCall(PetscQuadratureSetData(quad, dimR, 0, nq, quadPoints, NULL)); in DMFieldEvaluate_DS()
300 PetscCall(DMPlexComputeCellGeometryFEM(field->dm, c, quad, v, J, invJ, detJ)); in DMFieldEvaluate_DS()
301 PetscCall(PetscQuadratureDestroy(&quad)); in DMFieldEvaluate_DS()
455 PetscQuadrature quad = NULL; in DMFieldEvaluateFV_DS() local
484 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, pointIS, &quad)); in DMFieldEvaluateFV_DS()
485 if (!quad) PetscCall(DMFieldCreateDefaultQuadrature(field, pointIS, &quad)); in DMFieldEvaluateFV_DS()
[all …]
/petsc/src/dm/dt/fe/impls/vector/
H A Dfevector.c27 PetscQuadrature quad = NULL; in PetscFEView_Vector_Ascii() local
38 PetscCall(PetscFEGetQuadrature(fe, &quad)); in PetscFEView_Vector_Ascii()
54 if (quad) PetscCall(PetscQuadratureView(quad, v)); in PetscFEView_Vector_Ascii()
303 PetscQuadrature quad; in PetscFECreateVector() local
320 PetscCall(PetscFEGetQuadrature(scalar_fe, &quad)); in PetscFECreateVector()
321 PetscCall(PetscFESetQuadrature(fe_vec, quad)); in PetscFECreateVector()
322 PetscCall(PetscFEGetFaceQuadrature(scalar_fe, &quad)); in PetscFECreateVector()
323 PetscCall(PetscFESetFaceQuadrature(fe_vec, quad)); in PetscFECreateVector()
/petsc/src/dm/impls/moab/
H A Ddmmbfem.cxx110 …nt nverts, const PetscReal *coords, const PetscInt npts, const PetscReal *quad, PetscReal *phypts,… in Compute_Lagrange_Basis_1D_Internal() argument
126 const PetscReal r = quad[j]; in Compute_Lagrange_Basis_1D_Internal()
154 const PetscReal r = quad[j]; in Compute_Lagrange_Basis_1D_Internal()
220 …nt nverts, const PetscReal *coords, const PetscInt npts, const PetscReal *quad, PetscReal *phypts,… in Compute_Lagrange_Basis_2D_Internal() argument
238 const PetscReal r = quad[0 + j * 2]; in Compute_Lagrange_Basis_2D_Internal()
239 const PetscReal s = quad[1 + j * 2]; in Compute_Lagrange_Basis_2D_Internal()
301 const PetscReal r = quad[0 + j * 2]; in Compute_Lagrange_Basis_2D_Internal()
302 const PetscReal s = quad[1 + j * 2]; in Compute_Lagrange_Basis_2D_Internal()
376 …nt nverts, const PetscReal *coords, const PetscInt npts, const PetscReal *quad, PetscReal *phypts,… in Compute_Lagrange_Basis_3D_Internal() argument
397 const PetscReal &r = quad[j * 3 + 0]; in Compute_Lagrange_Basis_3D_Internal()
[all …]
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c18 PetscQuadrature quad = NULL; in PetscFEView_Basic_Ascii() local
25 PetscCall(PetscFEGetQuadrature(fe, &quad)); in PetscFEView_Basic_Ascii()
30 if (quad) PetscCall(PetscQuadratureView(quad, v)); in PetscFEView_Basic_Ascii()
167 PetscQuadrature quad; in PetscFEIntegrate_Basic() local
184 PetscCall(PetscFEGetQuadrature(fe, &quad)); in PetscFEIntegrate_Basic()
203 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in PetscFEIntegrate_Basic()
259 PetscQuadrature quad; in PetscFEIntegrateBd_Basic() local
275 PetscCall(PetscFEGetFaceQuadrature(fe, &quad)); in PetscFEIntegrateBd_Basic()
297 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in PetscFEIntegrateBd_Basic()
397 PetscQuadrature quad; in PetscFEIntegrateResidual_Basic() local
[all …]
/petsc/share/petsc/datafiles/meshes/
H A Dhybrid_3d_cube.msh111 Types: bar(1) tri(2) quad(3) tet(4) hex(5) prism(6) pyramid(7)
/petsc/src/dm/dt/interface/
H A Ddt.c584 static PetscErrorCode PetscQuadratureView_Ascii(PetscQuadrature quad, PetscViewer v) in PetscQuadratureView_Ascii() argument
590 if (quad->Nc > 1) in PetscQuadratureView_Ascii()
591 …" PetscInt_FMT " components\n", DMPolytopeTypes[quad->ct], quad->order, quad->numPoints, quad->dim… in PetscQuadratureView_Ascii()
592 … " points (dim %" PetscInt_FMT ")\n", DMPolytopeTypes[quad->ct], quad->order, quad->numPoints, qua… in PetscQuadratureView_Ascii()
595 for (q = 0; q < quad->numPoints; ++q) { in PetscQuadratureView_Ascii()
598 for (d = 0; d < quad->dim; ++d) { in PetscQuadratureView_Ascii()
600 PetscCall(PetscViewerASCIIPrintf(v, "%+g", (double)quad->points[q * quad->dim + d])); in PetscQuadratureView_Ascii()
603 if (quad->Nc > 1) PetscCall(PetscViewerASCIIPrintf(v, "w%" PetscInt_FMT " (", q)); in PetscQuadratureView_Ascii()
604 for (c = 0; c < quad->Nc; ++c) { in PetscQuadratureView_Ascii()
606 PetscCall(PetscViewerASCIIPrintf(v, "%+g", (double)quad->weights[q * quad->Nc + c])); in PetscQuadratureView_Ascii()
[all …]
/petsc/src/dm/impls/forest/p4est/
H A Dpforest.h2544 p4est_quadrant_t *quad = &fineQuads[i]; in DMPforestGetTransferSF_Point() local
2550 PetscCallP4estReturn(disjoint, p4est_quadrant_disjoint, (quadCoarse, quad)); in DMPforestGetTransferSF_Point()
2554 …if (quadCoarse->level > quad->level || (quadCoarse->level == quad->level && !transferIdent)) { /* … in DMPforestGetTransferSF_Point()
2563 … PetscCallP4est(p4est_quadrant_corner_descendant, (quad, &cornerQuad, j, quadCoarse->level)); in DMPforestGetTransferSF_Point()
2580 PetscCallP4estReturn(disjoint, p4est_quadrant_disjoint, (quadCoarse, quad)); in DMPforestGetTransferSF_Point()
2586 if (quadCoarse->level == quad->level) { /* same quad present in coarse and fine mesh */ in DMPforestGetTransferSF_Point()
2596 PetscInt levelDiff = quad->level - quadCoarse->level; in DMPforestGetTransferSF_Point()
2605 PetscCallP4estReturn(cid, p4est_quadrant_child_id, (quad)); in DMPforestGetTransferSF_Point()
2657 {quad->x, quad->y, in DMPforestGetTransferSF_Point()
2659 quad->z in DMPforestGetTransferSF_Point()
[all …]
/petsc/src/dm/impls/plex/
H A Dplexproject.c1019 PetscQuadrature quad = NULL; in DMProjectLocal_Generic_Plex() local
1029 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, isectIS, &quad)); in DMProjectLocal_Generic_Plex()
1030 if (!quad) { in DMProjectLocal_Generic_Plex()
1032 quad = allPoints; in DMProjectLocal_Generic_Plex()
1035 …PetscDualSpaceGetAllPointsUnion(Nf, sp, isCohesive ? dim - htInc - 1 : dim - htInc, funcs, &quad)); in DMProjectLocal_Generic_Plex()
1047 PetscCall(DMFieldCreateFEGeom(coordField, isectIS, quad, geommode, &fegeom)); in DMProjectLocal_Generic_Plex()
1060 PetscCall(PetscQuadratureDestroy(&quad)); in DMProjectLocal_Generic_Plex()
1066 PetscQuadrature quad = NULL; in DMProjectLocal_Generic_Plex() local
1071 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, pointIS, &quad)); in DMProjectLocal_Generic_Plex()
1072 if (!quad) { in DMProjectLocal_Generic_Plex()
[all …]
/petsc/src/dm/dt/fe/tests/
H A Dex2.c94 PetscErrorCode CellRangeGetFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomM… in CellRangeGetFEGeom() argument
101 PetscCall(PetscObjectGetId((PetscObject)quad, &id)); in CellRangeGetFEGeom()
107 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
113 PetscErrorCode CellRangeRestoreFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBoo… in CellRangeRestoreFEGeom() argument
H A Dex1.c103 PetscErrorCode CellRangeGetFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomM… in CellRangeGetFEGeom() argument
110 PetscCall(PetscObjectGetId((PetscObject)quad, &id)); in CellRangeGetFEGeom()
116 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
122 PetscErrorCode CellRangeRestoreFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBoo… in CellRangeRestoreFEGeom() argument
/petsc/src/dm/impls/swarm/
H A Dswarmpic_plex.c261 PetscQuadrature quad, facequad; in private_DMSwarmInsertPointsUsingCellDM_PLEX() local
268 PetscCall(PetscDTCreateDefaultQuadrature(ct, layout_param, &quad, &facequad)); in private_DMSwarmInsertPointsUsingCellDM_PLEX()
269 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, &xi, NULL)); in private_DMSwarmInsertPointsUsingCellDM_PLEX()
271 PetscCall(PetscQuadratureDestroy(&quad)); in private_DMSwarmInsertPointsUsingCellDM_PLEX()
/petsc/src/dm/impls/swarm/tests/
H A Dex2.c269 PetscQuadrature quad; in CreateParticles_Shape() local
282 PetscCall(PetscFEGetQuadrature(fe, &quad)); in CreateParticles_Shape()
283 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, &qpoints, &qweights)); in CreateParticles_Shape()
445 PetscQuadrature quad = NULL; in InterpolateGradient() local
470 PetscCall(PetscFEGetQuadrature(fe, &quad)); in InterpolateGradient()
475 PetscCall(PetscFVGetQuadrature(fv, &quad)); in InterpolateGradient()
480 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in InterpolateGradient()
497 …PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)… in InterpolateGradient()
/petsc/src/dm/impls/plex/tests/
H A Dex73.c114 PetscQuadrature quad, fquad; in CreateDiscretization() local
127 PetscCall(PetscFEGetQuadrature(feq, &quad)); in CreateDiscretization()
155 PetscCall(PetscFESetQuadrature(fep, quad)); in CreateDiscretization()

123