Lines Matching refs:cgeom
162 …scFEIntegrate_Basic(PetscDS ds, PetscInt field, PetscInt Ne, PetscFEGeom *cgeom, const PetscScalar… in PetscFEIntegrate_Basic() argument
205 Np = cgeom->numPoints; in PetscFEIntegrate_Basic()
206 dE = cgeom->dimEmbed; in PetscFEIntegrate_Basic()
207 isAffine = cgeom->isAffine; in PetscFEIntegrate_Basic()
211 fegeom.dim = cgeom->dim; in PetscFEIntegrate_Basic()
212 fegeom.dimEmbed = cgeom->dimEmbed; in PetscFEIntegrate_Basic()
216 fegeom.xi = cgeom->xi; in PetscFEIntegrate_Basic()
217 fegeom.J = &cgeom->J[e * Np * dE * dE]; in PetscFEIntegrate_Basic()
218 fegeom.invJ = &cgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrate_Basic()
219 fegeom.detJ = &cgeom->detJ[e * Np]; in PetscFEIntegrate_Basic()
226 …CoordinatesRefToReal(dE, dim, fegeom.xi, &cgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * dim], x… in PetscFEIntegrate_Basic()
228 fegeom.v = &cgeom->v[(e * Np + q) * dE]; in PetscFEIntegrate_Basic()
229 fegeom.J = &cgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
230 fegeom.invJ = &cgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
231 fegeom.detJ = &cgeom->detJ[e * Np + q]; in PetscFEIntegrate_Basic()
304 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBd_Basic() local
314 cgeom.dim = fgeom->dim; in PetscFEIntegrateBd_Basic()
315 cgeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
324 cgeom.J = &fgeom->suppJ[0][e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
325 cgeom.invJ = &fgeom->suppInvJ[0][e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
326 cgeom.detJ = &fgeom->suppDetJ[0][e * Np]; in PetscFEIntegrateBd_Basic()
341 cgeom.J = &fgeom->suppJ[0][(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
342 cgeom.invJ = &fgeom->suppInvJ[0][(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
343 cgeom.detJ = &fgeom->suppDetJ[0][e * Np + q]; in PetscFEIntegrateBd_Basic()
376 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, Tf, &cgeom, &coefficients[cOffset], N… in PetscFEIntegrateBd_Basic()
377 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, face, q, TfAux, &cgeom, &coef… in PetscFEIntegrateBd_Basic()
389 …ateResidual_Basic(PetscDS ds, PetscFormKey key, PetscInt Ne, PetscFEGeom *cgeom, const PetscScalar… in PetscFEIntegrateResidual_Basic() argument
437 dE = cgeom->dimEmbed; in PetscFEIntegrateResidual_Basic()
438 …etscCheck(cgeom->dim == qdim, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "FEGeom dim %" PetscInt_FMT "… in PetscFEIntegrateResidual_Basic()
449 PetscCall(PetscFEGeomGetPoint(cgeom, e, q, &quadPoints[q * cgeom->dim], &fegeom)); in PetscFEIntegrateResidual_Basic()
452 if (debug > 1 && q < cgeom->numPoints) { in PetscFEIntegrateResidual_Basic()
489 …PetscCall(PetscFEUpdateElementVec_Internal(fe, T[field], 0, basisReal, basisDerReal, e, cgeom, f0,… in PetscFEIntegrateResidual_Basic()
552 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBdResidual_Basic() local
563 PetscCall(PetscFEGeomGetCellPoint(fgeom, e, q, &cgeom)); in PetscFEIntegrateBdResidual_Basic()
575 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, Tf, &cgeom, &coefficients[cOffset], P… in PetscFEIntegrateBdResidual_Basic()
576 …FEEvaluateFieldJets_Internal(dsAux, NfAux, auxOnBd ? 0 : face, q, TfAux, &cgeom, &coefficientsAux[… in PetscFEIntegrateBdResidual_Basic()
743 …s, PetscFEJacobianType jtype, PetscFormKey key, PetscInt Ne, PetscFEGeom *cgeom, const PetscScalar… in PetscFEIntegrateJacobian_Basic() argument
817 Np = cgeom->numPoints; in PetscFEIntegrateJacobian_Basic()
818 dE = cgeom->dimEmbed; in PetscFEIntegrateJacobian_Basic()
819 isAffine = cgeom->isAffine; in PetscFEIntegrateJacobian_Basic()
826 fegeom.dim = cgeom->dim; in PetscFEIntegrateJacobian_Basic()
827 fegeom.dimEmbed = cgeom->dimEmbed; in PetscFEIntegrateJacobian_Basic()
831 fegeom.xi = cgeom->xi; in PetscFEIntegrateJacobian_Basic()
832 fegeom.J = &cgeom->J[e * Np * dE * dE]; in PetscFEIntegrateJacobian_Basic()
833 fegeom.invJ = &cgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrateJacobian_Basic()
834 fegeom.detJ = &cgeom->detJ[e * Np]; in PetscFEIntegrateJacobian_Basic()
840 …CoordinatesRefToReal(dE, dim, fegeom.xi, &cgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * dim], x… in PetscFEIntegrateJacobian_Basic()
842 fegeom.v = &cgeom->v[(e * Np + q) * dE]; in PetscFEIntegrateJacobian_Basic()
843 fegeom.J = &cgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrateJacobian_Basic()
844 fegeom.invJ = &cgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrateJacobian_Basic()
845 fegeom.detJ = &cgeom->detJ[e * Np + q]; in PetscFEIntegrateJacobian_Basic()
967 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBdJacobian_Basic() local
976 cgeom.dim = fgeom->dim; in PetscFEIntegrateBdJacobian_Basic()
977 cgeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
986 cgeom.J = &fgeom->suppJ[0][e * Np * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
987 cgeom.invJ = &fgeom->suppInvJ[0][e * Np * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
988 cgeom.detJ = &fgeom->suppDetJ[0][e * Np]; in PetscFEIntegrateBdJacobian_Basic()
1004 cgeom.J = &fgeom->suppJ[0][(e * Np + q) * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
1005 cgeom.invJ = &fgeom->suppInvJ[0][(e * Np + q) * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
1006 cgeom.detJ = &fgeom->suppDetJ[0][e * Np + q]; in PetscFEIntegrateBdJacobian_Basic()
1010 …if (coefficients) PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, T, &cgeom, &coeffic… in PetscFEIntegrateBdJacobian_Basic()
1011 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, face, q, TAux, &cgeom, &coeff… in PetscFEIntegrateBdJacobian_Basic()
1033 …q, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &cgeom, g0, g1, g2, g3, t… in PetscFEIntegrateBdJacobian_Basic()