Lines Matching refs:fegeom
209 PetscFEGeom fegeom; in PetscFEIntegrate_Basic() local
211 fegeom.dim = cgeom->dim; in PetscFEIntegrate_Basic()
212 fegeom.dimEmbed = cgeom->dimEmbed; in PetscFEIntegrate_Basic()
213 fegeom.xi = NULL; in PetscFEIntegrate_Basic()
215 fegeom.v = x; 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()
220 } else fegeom.xi = NULL; 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()
233 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrate_Basic()
234 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrate_Basic()
236 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrate_Basic()
238 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrate_Basic()
242 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, 0, q, T, &fegeom, &coefficients[cOffset], NULL… in PetscFEIntegrate_Basic()
243 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, 0, q, TAux, &fegeom, &coeffic… in PetscFEIntegrate_Basic()
244 …obj_func(dim, Nf, NfAux, uOff, uOff_x, u, NULL, u_x, aOff, aOff_x, a, NULL, a_x, 0.0, fegeom.v, nu… in PetscFEIntegrate_Basic()
304 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBd_Basic() local
306 fegeom.n = NULL; in PetscFEIntegrateBd_Basic()
307 fegeom.v = NULL; in PetscFEIntegrateBd_Basic()
308 fegeom.xi = NULL; in PetscFEIntegrateBd_Basic()
309 fegeom.J = NULL; in PetscFEIntegrateBd_Basic()
310 fegeom.invJ = NULL; in PetscFEIntegrateBd_Basic()
311 fegeom.detJ = NULL; in PetscFEIntegrateBd_Basic()
312 fegeom.dim = fgeom->dim; in PetscFEIntegrateBd_Basic()
313 fegeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
317 fegeom.v = x; 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()
327 } else fegeom.xi = NULL; 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()
345 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBd_Basic()
346 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateBd_Basic()
348 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrateBd_Basic()
350 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrateBd_Basic()
358 PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)fegeom.v[d])); in PetscFEIntegrateBd_Basic()
364 PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)fegeom.n[d])); in PetscFEIntegrateBd_Basic()
378 …NfAux, uOff, uOff_x, u, NULL, u_x, aOff, aOff_x, a, NULL, a_x, 0.0, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBd_Basic()
440 PetscFEGeom fegeom; in PetscFEIntegrateResidual_Basic() local
442 fegeom.v = x; /* workspace */ in PetscFEIntegrateResidual_Basic()
449 PetscCall(PetscFEGeomGetPoint(cgeom, e, q, &quadPoints[q * cgeom->dim], &fegeom)); in PetscFEIntegrateResidual_Basic()
450 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateResidual_Basic()
451 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateResidual_Basic()
453 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrateResidual_Basic()
455 PetscCall(DMPrintCellMatrix(e, "invJ", dE, dE, fegeom.invJ)); in PetscFEIntegrateResidual_Basic()
458 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, 0, q, T, &fegeom, &coefficients[cOffset], Pets… in PetscFEIntegrateResidual_Basic()
459 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, 0, q, TAux, &fegeom, &coeffic… in PetscFEIntegrateResidual_Basic()
460 …dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, numConstants, c… in PetscFEIntegrateResidual_Basic()
462 …dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, numConstants, c… in PetscFEIntegrateResidual_Basic()
468 … for (c = 0; c < dE; ++c) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)fegeom.v[c])); in PetscFEIntegrateResidual_Basic()
552 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBdResidual_Basic() local
555 fegeom.v = x; /* Workspace */ in PetscFEIntegrateBdResidual_Basic()
562 PetscCall(PetscFEGeomGetPoint(fgeom, e, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateBdResidual_Basic()
564 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBdResidual_Basic()
565 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateBdResidual_Basic()
568 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrateBdResidual_Basic()
570 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrateBdResidual_Basic()
571 PetscCall(DMPrintCellVector(e, "n", dim, fegeom.n)); in PetscFEIntegrateBdResidual_Basic()
577 …f, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdResidual_Basic()
579 …f, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdResidual_Basic()
677 PetscFEGeom fegeom, fegeomN[2]; in PetscFEIntegrateHybridResidual_Basic() local
682 fegeom.v = x; /* Workspace */ in PetscFEIntegrateHybridResidual_Basic()
696 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridResidual_Basic()
699 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateHybridResidual_Basic()
701 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrateHybridResidual_Basic()
703 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dE, fegeom.invJ)); in PetscFEIntegrateHybridResidual_Basic()
706 … point %" PetscInt_FMT " weight %g detJ %g\n", q, (double)quadWeights[q], (double)fegeom.detJ[0])); in PetscFEIntegrateHybridResidual_Basic()
708 …EEvaluateFieldJets_Hybrid_Internal(dsIn, Nf, 0, q, Tf, face, qpt, TfIn, &fegeom, fegeomN, &coeffic… in PetscFEIntegrateHybridResidual_Basic()
709 …ernal(dsAux, NfAux, auxOnBd ? 0 : face[s], auxOnBd ? q : qpt[s], TfAux, &fegeom, &coefficientsAux[… in PetscFEIntegrateHybridResidual_Basic()
710 …f, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridResidual_Basic()
712 …f, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridResidual_Basic()
824 PetscFEGeom fegeom; in PetscFEIntegrateJacobian_Basic() local
826 fegeom.dim = cgeom->dim; in PetscFEIntegrateJacobian_Basic()
827 fegeom.dimEmbed = cgeom->dimEmbed; in PetscFEIntegrateJacobian_Basic()
828 fegeom.xi = NULL; in PetscFEIntegrateJacobian_Basic()
830 fegeom.v = x; 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()
835 } else fegeom.xi = NULL; 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()
847 PetscCall(PetscDSSetCellParameters(rds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateJacobian_Basic()
848 … point %" PetscInt_FMT " weight %g detJ %g\n", q, (double)quadWeights[q], (double)fegeom.detJ[0])); in PetscFEIntegrateJacobian_Basic()
849 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateJacobian_Basic()
850 …if (coefficients) PetscCall(PetscFEEvaluateFieldJets_Internal(rds, Nf, 0, q, rT, &fegeom, &coeffic… in PetscFEIntegrateJacobian_Basic()
851 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, 0, q, TAux, &fegeom, &coeffic… in PetscFEIntegrateJacobian_Basic()
854 …Aux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, c… in PetscFEIntegrateJacobian_Basic()
859 …Aux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, c… in PetscFEIntegrateJacobian_Basic()
864 …Aux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, c… in PetscFEIntegrateJacobian_Basic()
869 …Aux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, c… in PetscFEIntegrateJacobian_Basic()
873 …rT[fieldI], basisReal, basisDerReal, cT[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, c… in PetscFEIntegrateJacobian_Basic()
967 PetscFEGeom fegeom, cgeom; in PetscFEIntegrateBdJacobian_Basic() local
969 fegeom.n = NULL; in PetscFEIntegrateBdJacobian_Basic()
970 fegeom.v = NULL; in PetscFEIntegrateBdJacobian_Basic()
971 fegeom.xi = NULL; in PetscFEIntegrateBdJacobian_Basic()
972 fegeom.J = NULL; in PetscFEIntegrateBdJacobian_Basic()
973 fegeom.detJ = NULL; in PetscFEIntegrateBdJacobian_Basic()
974 fegeom.dim = fgeom->dim; in PetscFEIntegrateBdJacobian_Basic()
975 fegeom.dimEmbed = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
979 fegeom.v = x; 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()
989 } else fegeom.xi = NULL; 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()
1008 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBdJacobian_Basic()
1009 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateBdJacobian_Basic()
1014 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdJacobian_Basic()
1019 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdJacobian_Basic()
1024 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdJacobian_Basic()
1029 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateBdJacobian_Basic()
1154 PetscFEGeom fegeom, fegeomN[2]; in PetscFEIntegrateHybridJacobian_Basic() local
1159 fegeom.v = x; /* Workspace */ in PetscFEIntegrateHybridJacobian_Basic()
1167 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridJacobian_Basic()
1170 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateHybridJacobian_Basic()
1172 PetscCall(PetscPrintf(PETSC_COMM_SELF, " detJ: %g\n", (double)fegeom.detJ[0])); in PetscFEIntegrateHybridJacobian_Basic()
1174 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrateHybridJacobian_Basic()
1178 …FEEvaluateFieldJets_Hybrid_Internal(dsIn, Nf, 0, q, T, face, qpt, TfIn, &fegeom, fegeomN, &coeffic… in PetscFEIntegrateHybridJacobian_Basic()
1179 …ternal(dsAux, NfAux, auxOnBd ? 0 : face[s], auxOnBd ? q : qpt[s], TAux, &fegeom, &coefficientsAux[… in PetscFEIntegrateHybridJacobian_Basic()
1182 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridJacobian_Basic()
1187 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridJacobian_Basic()
1192 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridJacobian_Basic()
1197 …uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, fegeom.n, numConstan… in PetscFEIntegrateHybridJacobian_Basic()
1204 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, e… in PetscFEIntegrateHybridJacobian_Basic()
1206 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, e… in PetscFEIntegrateHybridJacobian_Basic()
1207 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, &g0[NcI * NcJ], &… in PetscFEIntegrateHybridJacobian_Basic()
1212 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, e… in PetscFEIntegrateHybridJacobian_Basic()
1213 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, &g0[NcI * NcJ], &… in PetscFEIntegrateHybridJacobian_Basic()
1215 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, e… in PetscFEIntegrateHybridJacobian_Basic()
1216 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, &g0[NcI * NcJ], &… in PetscFEIntegrateHybridJacobian_Basic()
1217 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, &g0[NcI * NcJ * 2… in PetscFEIntegrateHybridJacobian_Basic()
1218 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, &g0[NcI * NcJ * 3… in PetscFEIntegrateHybridJacobian_Basic()
1221 …, T[fieldI], basisReal, basisDerReal, T[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, e… in PetscFEIntegrateHybridJacobian_Basic()