Lines Matching refs:fgeom
255 …S ds, PetscInt field, PetscBdPointFn *obj_func, PetscInt Ne, PetscFEGeom *fgeom, const PetscScalar… in PetscFEIntegrateBd_Basic() argument
300 Np = fgeom->numPoints; in PetscFEIntegrateBd_Basic()
301 dE = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
302 isAffine = fgeom->isAffine; in PetscFEIntegrateBd_Basic()
305 const PetscInt face = fgeom->face[e][0]; /* Local face number in cell */ in PetscFEIntegrateBd_Basic()
312 fegeom.dim = fgeom->dim; in PetscFEIntegrateBd_Basic()
313 fegeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
314 cgeom.dim = fgeom->dim; in PetscFEIntegrateBd_Basic()
315 cgeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
318 fegeom.xi = fgeom->xi; in PetscFEIntegrateBd_Basic()
319 fegeom.J = &fgeom->J[e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
320 fegeom.invJ = &fgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
321 fegeom.detJ = &fgeom->detJ[e * Np]; in PetscFEIntegrateBd_Basic()
322 fegeom.n = &fgeom->n[e * Np * dE]; 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()
333 …CoordinatesRefToReal(dE, dim - 1, fegeom.xi, &fgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * (di… in PetscFEIntegrateBd_Basic()
335 fegeom.v = &fgeom->v[(e * Np + q) * dE]; in PetscFEIntegrateBd_Basic()
336 fegeom.J = &fgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
337 fegeom.invJ = &fgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
338 fegeom.detJ = &fgeom->detJ[e * Np + q]; in PetscFEIntegrateBd_Basic()
339 fegeom.n = &fgeom->n[(e * Np + q) * dE]; 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()
496 …PetscDS ds, PetscWeakForm wf, PetscFormKey key, PetscInt Ne, PetscFEGeom *fgeom, const PetscScalar… in PetscFEIntegrateBdResidual_Basic() argument
547 dE = fgeom->dimEmbed; in PetscFEIntegrateBdResidual_Basic()
549 fgeom->dim = dim - 1; in PetscFEIntegrateBdResidual_Basic()
550 …etscCheck(fgeom->dim == qdim, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "FEGeom dim %" PetscInt_FMT "… in PetscFEIntegrateBdResidual_Basic()
553 const PetscInt face = fgeom->face[e][0]; in PetscFEIntegrateBdResidual_Basic()
562 PetscCall(PetscFEGeomGetPoint(fgeom, e, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateBdResidual_Basic()
563 PetscCall(PetscFEGeomGetCellPoint(fgeom, e, q, &cgeom)); in PetscFEIntegrateBdResidual_Basic()
567 if ((fgeom->isAffine && q == 0) || !fgeom->isAffine) { in PetscFEIntegrateBdResidual_Basic()
593 …pdateElementVec_Internal(fe, Tf[field], face, basisReal, basisDerReal, e, fgeom, f0, f1, &elemVec[… in PetscFEIntegrateBdResidual_Basic()
609 …ds, PetscDS dsIn, PetscFormKey key, PetscInt s, PetscInt Ne, PetscFEGeom *fgeom, PetscFEGeom *nbrg… in PetscFEIntegrateHybridResidual_Basic() argument
673 dE = fgeom->dimEmbed; in PetscFEIntegrateHybridResidual_Basic()
674 …etscCheck(fgeom->dim == qdim, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "FEGeom dim %" PetscInt_FMT "… in PetscFEIntegrateHybridResidual_Basic()
678 const PetscInt face[2] = {fgeom->face[e * 2 + 0][0], fgeom->face[e * 2 + 1][2]}; in PetscFEIntegrateHybridResidual_Basic()
679 const PetscInt ornt[2] = {fgeom->face[e * 2 + 0][1], fgeom->face[e * 2 + 1][3]}; in PetscFEIntegrateHybridResidual_Basic()
680 const PetscInt cornt[2] = {fgeom->face[e * 2 + 0][3], fgeom->face[e * 2 + 1][1]}; in PetscFEIntegrateHybridResidual_Basic()
696 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridResidual_Basic()
700 if (debug > 1 && q < fgeom->numPoints) { in PetscFEIntegrateHybridResidual_Basic()
732 …FEUpdateElementVec_Internal(fe, Tf[field], 0, basisReal, basisDerReal, e, fgeom, f0, f1, &elemVec[… in PetscFEIntegrateHybridResidual_Basic()
734 …eElementVec_Hybrid_Internal(fe, Tf[field], 0, s, basisReal, basisDerReal, fgeom, f0, f1, &elemVec[… in PetscFEIntegrateHybridResidual_Basic()
893 …f, PetscFEJacobianType jtype, PetscFormKey key, PetscInt Ne, PetscFEGeom *fgeom, const PetscScalar… in PetscFEIntegrateBdJacobian_Basic() argument
956 Np = fgeom->numPoints; in PetscFEIntegrateBdJacobian_Basic()
957 dE = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
958 isAffine = fgeom->isAffine; in PetscFEIntegrateBdJacobian_Basic()
968 const PetscInt face = fgeom->face[e][0]; in PetscFEIntegrateBdJacobian_Basic()
974 fegeom.dim = fgeom->dim; in PetscFEIntegrateBdJacobian_Basic()
975 fegeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
976 cgeom.dim = fgeom->dim; in PetscFEIntegrateBdJacobian_Basic()
977 cgeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
980 fegeom.xi = fgeom->xi; in PetscFEIntegrateBdJacobian_Basic()
981 fegeom.J = &fgeom->J[e * Np * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
982 fegeom.invJ = &fgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
983 fegeom.detJ = &fgeom->detJ[e * Np]; in PetscFEIntegrateBdJacobian_Basic()
984 fegeom.n = &fgeom->n[e * Np * dE]; 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()
996 …CoordinatesRefToReal(dE, dim - 1, fegeom.xi, &fgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * (di… in PetscFEIntegrateBdJacobian_Basic()
998 fegeom.v = &fgeom->v[(e * Np + q) * dE]; in PetscFEIntegrateBdJacobian_Basic()
999 fegeom.J = &fgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
1000 fegeom.invJ = &fgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrateBdJacobian_Basic()
1001 fegeom.detJ = &fgeom->detJ[e * Np + q]; in PetscFEIntegrateBdJacobian_Basic()
1002 fegeom.n = &fgeom->n[(e * Np + q) * dE]; 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()
1059 …cobianType jtype, PetscFormKey key, PetscInt s, PetscInt Ne, PetscFEGeom *fgeom, PetscFEGeom *nbrg… in PetscFEIntegrateHybridJacobian_Basic() argument
1131 dE = fgeom->dimEmbed; in PetscFEIntegrateHybridJacobian_Basic()
1155 const PetscInt face[2] = {fgeom->face[e * 2 + 0][0], fgeom->face[e * 2 + 1][2]}; in PetscFEIntegrateHybridJacobian_Basic()
1156 const PetscInt ornt[2] = {fgeom->face[e * 2 + 0][1], fgeom->face[e * 2 + 1][3]}; in PetscFEIntegrateHybridJacobian_Basic()
1157 const PetscInt cornt[2] = {fgeom->face[e * 2 + 0][3], fgeom->face[e * 2 + 1][1]}; in PetscFEIntegrateHybridJacobian_Basic()
1167 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridJacobian_Basic()
1171 if (debug > 1 && q < fgeom->numPoints) { in PetscFEIntegrateHybridJacobian_Basic()