Lines Matching refs:fe
806 PetscErrorCode PetscFEExpandFaceQuadrature(PetscFE fe, PetscQuadrature fq, PetscQuadrature *efq) in PetscFEExpandFaceQuadrature() argument
817 PetscCall(PetscFEGetDualSpace(fe, &sp)); in PetscFEExpandFaceQuadrature()
1095 static PetscErrorCode PetscFECreatePointTraceDefault_Internal(PetscFE fe, PetscInt refPoint, PetscF… in PetscFECreatePointTraceDefault_Internal() argument
1107 PetscCall(PetscFEGetBasisSpace(fe, &bsp)); in PetscFECreatePointTraceDefault_Internal()
1108 PetscCall(PetscFEGetDualSpace(fe, &dsp)); in PetscFECreatePointTraceDefault_Internal()
1128 PetscCall(PetscFECreate(PetscObjectComm((PetscObject)fe), trFE)); in PetscFECreatePointTraceDefault_Internal()
1130 PetscCall(PetscFEGetNumComponents(fe, &numComp)); in PetscFECreatePointTraceDefault_Internal()
1134 PetscCall(PetscObjectGetName((PetscObject)fe, &name)); in PetscFECreatePointTraceDefault_Internal()
1136 PetscCall(PetscFEGetQuadrature(fe, &fullQuad)); in PetscFECreatePointTraceDefault_Internal()
1148 PETSC_EXTERN PetscErrorCode PetscFECreatePointTrace(PetscFE fe, PetscInt refPoint, PetscFE *trFE) in PetscFECreatePointTrace() argument
1151 PetscValidHeaderSpecific(fe, PETSCFE_CLASSID, 1); in PetscFECreatePointTrace()
1153 if (fe->ops->createpointtrace) PetscUseTypeMethod(fe, createpointtrace, refPoint, trFE); in PetscFECreatePointTrace()
1154 else PetscCall(PetscFECreatePointTraceDefault_Internal(fe, refPoint, trFE)); in PetscFECreatePointTrace()
1158 PetscErrorCode PetscFECreateHeightTrace(PetscFE fe, PetscInt height, PetscFE *trFE) in PetscFECreateHeightTrace() argument
1165 PetscValidHeaderSpecific(fe, PETSCFE_CLASSID, 1); in PetscFECreateHeightTrace()
1168 PetscCall(PetscFEGetDualSpace(fe, &dsp)); in PetscFECreateHeightTrace()
1172 PetscCall(PetscFECreatePointTrace(fe, hStart, trFE)); in PetscFECreateHeightTrace()
1221 PetscErrorCode PetscFEPushforward(PetscFE fe, PetscFEGeom *fegeom, PetscInt Nv, PetscScalar vals[]) in PetscFEPushforward() argument
1224 PetscCall(PetscDualSpacePushforward(fe->dualSpace, fegeom, Nv, fe->numComponents, vals)); in PetscFEPushforward()
1249 PetscErrorCode PetscFEPushforwardGradient(PetscFE fe, PetscFEGeom *fegeom, PetscInt Nv, PetscScalar… in PetscFEPushforwardGradient() argument
1252 PetscCall(PetscDualSpacePushforwardGradient(fe->dualSpace, fegeom, Nv, fe->numComponents, vals)); in PetscFEPushforwardGradient()
1280 PetscErrorCode PetscFEPushforwardHessian(PetscFE fe, PetscFEGeom *fegeom, PetscInt Nv, PetscScalar … in PetscFEPushforwardHessian() argument
1283 PetscCall(PetscDualSpacePushforwardHessian(fe->dualSpace, fegeom, Nv, fe->numComponents, vals)); in PetscFEPushforwardHessian()
1401 PetscFE fe; in PetscFEIntegrate() local
1405 PetscCall(PetscDSGetDiscretization(prob, field, (PetscObject *)&fe)); in PetscFEIntegrate()
1406 …if (fe->ops->integrate) PetscCall((*fe->ops->integrate)(prob, field, Ne, cgeom, coefficients, prob… in PetscFEIntegrate()
1437 PetscFE fe; in PetscFEIntegrateBd() local
1441 PetscCall(PetscDSGetDiscretization(prob, field, (PetscObject *)&fe)); in PetscFEIntegrateBd()
1442 …if (fe->ops->integratebd) PetscCall((*fe->ops->integratebd)(prob, field, obj_func, Ne, geom, coeff… in PetscFEIntegrateBd()
1481 PetscFE fe; in PetscFEIntegrateResidual() local
1485 PetscCall(PetscDSGetDiscretization(ds, key.field, (PetscObject *)&fe)); in PetscFEIntegrateResidual()
1486 …if (fe->ops->integrateresidual) PetscCall((*fe->ops->integrateresidual)(ds, key, Ne, cgeom, coeffi… in PetscFEIntegrateResidual()
1516 PetscFE fe; in PetscFEIntegrateBdResidual() local
1520 PetscCall(PetscDSGetDiscretization(ds, key.field, (PetscObject *)&fe)); in PetscFEIntegrateBdResidual()
1521 …if (fe->ops->integratebdresidual) PetscCall((*fe->ops->integratebdresidual)(ds, wf, key, Ne, fgeom… in PetscFEIntegrateBdResidual()
1553 PetscFE fe; in PetscFEIntegrateHybridResidual() local
1558 PetscCall(PetscDSGetDiscretization(ds, key.field, (PetscObject *)&fe)); in PetscFEIntegrateHybridResidual()
1559 …if (fe->ops->integratehybridresidual) PetscCall((*fe->ops->integratehybridresidual)(ds, dsIn, key,… in PetscFEIntegrateHybridResidual()
1603 PetscFE fe; in PetscFEIntegrateJacobian() local
1610 PetscCall(PetscDSGetDiscretization(rds, key.field / Nf, (PetscObject *)&fe)); in PetscFEIntegrateJacobian()
1611 …if (fe->ops->integratejacobian) PetscCall((*fe->ops->integratejacobian)(rds, cds, jtype, key, Ne, … in PetscFEIntegrateJacobian()
1655 PetscFE fe; in PetscFEIntegrateBdJacobian() local
1661 PetscCall(PetscDSGetDiscretization(ds, key.field / Nf, (PetscObject *)&fe)); in PetscFEIntegrateBdJacobian()
1662 …if (fe->ops->integratebdjacobian) PetscCall((*fe->ops->integratebdjacobian)(ds, wf, jtype, key, Ne… in PetscFEIntegrateBdJacobian()
1708 PetscFE fe; in PetscFEIntegrateHybridJacobian() local
1714 PetscCall(PetscDSGetDiscretization(ds, key.field / Nf, (PetscObject *)&fe)); in PetscFEIntegrateHybridJacobian()
1715 …if (fe->ops->integratehybridjacobian) PetscCall((*fe->ops->integratehybridjacobian)(ds, dsIn, jtyp… in PetscFEIntegrateHybridJacobian()
1736 PetscErrorCode PetscFEGetHeightSubspace(PetscFE fe, PetscInt height, PetscFE *subfe) in PetscFEGetHeightSubspace() argument
1744 PetscValidHeaderSpecific(fe, PETSCFE_CLASSID, 1); in PetscFEGetHeightSubspace()
1747 *subfe = fe; in PetscFEGetHeightSubspace()
1750 PetscCall(PetscFEGetBasisSpace(fe, &P)); in PetscFEGetHeightSubspace()
1751 PetscCall(PetscFEGetDualSpace(fe, &Q)); in PetscFEGetHeightSubspace()
1752 PetscCall(PetscFEGetNumComponents(fe, &Nc)); in PetscFEGetHeightSubspace()
1753 PetscCall(PetscFEGetFaceQuadrature(fe, &subq)); in PetscFEGetHeightSubspace()
1756 if (!fe->subspaces) PetscCall(PetscCalloc1(dim, &fe->subspaces)); in PetscFEGetHeightSubspace()
1758 if (!fe->subspaces[height - 1]) { in PetscFEGetHeightSubspace()
1771 PetscCall(PetscObjectGetName((PetscObject)fe, &name)); in PetscFEGetHeightSubspace()
1774 fe->subspaces[height - 1] = sub; in PetscFEGetHeightSubspace()
1776 *subfe = fe->subspaces[height - 1]; in PetscFEGetHeightSubspace()
1804 PetscErrorCode PetscFERefine(PetscFE fe, PetscFE *feRef) in PetscFERefine() argument
1816 PetscCall(PetscFEGetBasisSpace(fe, &P)); in PetscFERefine()
1817 PetscCall(PetscFEGetDualSpace(fe, &Q)); in PetscFERefine()
1818 PetscCall(PetscFEGetQuadrature(fe, &q)); in PetscFERefine()
1826 PetscCall(DMRefine(K, PetscObjectComm((PetscObject)fe), &Kref)); in PetscFERefine()
1838 PetscCall(PetscFECreate(PetscObjectComm((PetscObject)fe), feRef)); in PetscFERefine()
1842 PetscCall(PetscFEGetNumComponents(fe, &numComp)); in PetscFERefine()
1855 static PetscErrorCode PetscFESetDefaultName_Private(PetscFE fe) in PetscFESetDefaultName_Private() argument
1865 PetscCall(PetscFEGetBasisSpace(fe, &P)); in PetscFESetDefaultName_Private()
1867 PetscCall(PetscFEGetDualSpace(fe, &Q)); in PetscFESetDefaultName_Private()
1890 PetscCall(PetscFESetName(fe, name)); in PetscFESetDefaultName_Private()
2240 PetscErrorCode PetscFELimitDegree(PetscFE fe, PetscInt minDegree, PetscInt maxDegree, PetscFE *newf… in PetscFELimitDegree() argument
2247 PetscCall(PetscFEGetDualSpace(fe, &Q)); in PetscFELimitDegree()
2258 PetscCall(PetscObjectReference((PetscObject)fe)); in PetscFELimitDegree()
2259 *newfe = fe; in PetscFELimitDegree()
2272 PetscCall(PetscFEGetBasisSpace(fe, &P)); in PetscFELimitDegree()
2277 …PetscCall(PetscFECreateLagrangeByCell(PetscObjectComm((PetscObject)fe), dim, Nc, ct, k, PETSC_DETE… in PetscFELimitDegree()
2278 PetscCall(PetscFEGetQuadrature(fe, &q)); in PetscFELimitDegree()
2281 PetscCall(PetscObjectReference((PetscObject)fe)); in PetscFELimitDegree()
2282 *newfe = fe; in PetscFELimitDegree()
2343 PetscErrorCode PetscFESetName(PetscFE fe, const char name[]) in PetscFESetName() argument
2349 PetscCall(PetscFEGetBasisSpace(fe, &P)); in PetscFESetName()
2350 PetscCall(PetscFEGetDualSpace(fe, &Q)); in PetscFESetName()
2351 PetscCall(PetscObjectSetName((PetscObject)fe, name)); in PetscFESetName()
2364 PetscFE fe; in PetscFEEvaluateFieldJets_Internal() local
2376 PetscCall(PetscDSGetDiscretization(ds, f, (PetscObject *)&fe)); in PetscFEEvaluateFieldJets_Internal()
2397 PetscCall(PetscFEPushforwardHessian(fe, fegeom, 1, &u_x[hOffset + fOffset * dE * dE])); in PetscFEEvaluateFieldJets_Internal()
2399 PetscCall(PetscFEPushforward(fe, fegeom, 1, &u[fOffset])); in PetscFEEvaluateFieldJets_Internal()
2400 PetscCall(PetscFEPushforwardGradient(fe, fegeom, 1, &u_x[fOffset * dE])); in PetscFEEvaluateFieldJets_Internal()
2410 PetscCall(PetscFEPushforward(fe, fegeom, 1, &u_t[fOffset])); in PetscFEEvaluateFieldJets_Internal()
2432 PetscFE fe = (PetscFE)ds->disc[f]; in PetscFEEvaluateFieldJets_Hybrid_Internal() local
2458 PetscCall(PetscFEPushforward(fe, isCohesive ? fegeom : &fegeomNbr[s], 1, &u[fOffset])); in PetscFEEvaluateFieldJets_Hybrid_Internal()
2459 …PetscCall(PetscFEPushforwardGradient(fe, isCohesive ? fegeom : &fegeomNbr[s], 1, &u_x[fOffset * dE… in PetscFEEvaluateFieldJets_Hybrid_Internal()
2469 PetscCall(PetscFEPushforward(fe, fegeom, 1, &u_t[fOffset])); in PetscFEEvaluateFieldJets_Hybrid_Internal()
2481 PetscFE fe; in PetscFEEvaluateFaceFields_Internal() local
2486 PetscCall(PetscDSGetDiscretization(prob, field, (PetscObject *)&fe)); in PetscFEEvaluateFaceFields_Internal()
2487 PetscCall(PetscFEGetFaceCentroidTabulation(fe, &Tc)); in PetscFEEvaluateFaceFields_Internal()
2501 PetscErrorCode PetscFEUpdateElementVec_Internal(PetscFE fe, PetscTabulation T, PetscInt r, PetscSca… in PetscFEUpdateElementVec_Internal() argument
2524 PetscCall(PetscFEPushforward(fe, &pgeom, Nb, tmpBasis)); in PetscFEUpdateElementVec_Internal()
2525 PetscCall(PetscFEPushforwardGradient(fe, &pgeom, Nb, tmpBasisDer)); in PetscFEUpdateElementVec_Internal()
2539 PetscErrorCode PetscFEUpdateElementVec_Hybrid_Internal(PetscFE fe, PetscTabulation T, PetscInt r, P… in PetscFEUpdateElementVec_Hybrid_Internal() argument
2557 PetscCall(PetscFEPushforward(fe, fegeom, Nb, tmpBasis)); in PetscFEUpdateElementVec_Hybrid_Internal()
2861 PetscErrorCode PetscFECreateCellGeometry(PetscFE fe, PetscQuadrature quad, PetscFEGeom *cgeom) in PetscFECreateCellGeometry() argument
2869 PetscCall(PetscFEGetDualSpace(fe, &dsp)); in PetscFECreateCellGeometry()
2873 PetscCall(PetscFEGetQuadrature(fe, &quadDef)); in PetscFECreateCellGeometry()
2888 PetscErrorCode PetscFEDestroyCellGeometry(PetscFE fe, PetscFEGeom *cgeom) in PetscFEDestroyCellGeometry() argument