Lines Matching refs:dE
173 …PetscInt dim, dE, Np, numConstants, Nf, NfAux = 0, totDim, totDimAux = 0, cOffset = 0, c… in PetscFEIntegrate_Basic() local
206 dE = cgeom->dimEmbed; 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()
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()
267 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBd_Basic() local
301 dE = fgeom->dimEmbed; 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()
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()
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()
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()
356 for (PetscInt d = 0; d < dE; ++d) { in PetscFEIntegrateBd_Basic()
362 for (PetscInt d = 0; d < dE; ++d) { in PetscFEIntegrateBd_Basic()
405 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateResidual_Basic() local
437 dE = cgeom->dimEmbed; in PetscFEIntegrateResidual_Basic()
444 PetscCall(PetscArrayzero(f1, Nq * T[field]->Nc * dE)); in PetscFEIntegrateResidual_Basic()
455 PetscCall(DMPrintCellMatrix(e, "invJ", dE, dE, fegeom.invJ)); in PetscFEIntegrateResidual_Basic()
460 …for (i = 0; i < n0; ++i) f0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateResidual_Basic()
462 …1_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, numCo… in PetscFEIntegrateResidual_Basic()
464 for (d = 0; d < dE; ++d) f1[(q * T[field]->Nc + c) * dE + d] *= w; in PetscFEIntegrateResidual_Basic()
468 … for (c = 0; c < dE; ++c) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)fegeom.v[c])); in PetscFEIntegrateResidual_Basic()
475 …for (c = 0; c < T[field]->Nc * dE; ++c) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)Pets… in PetscFEIntegrateResidual_Basic()
482 …for (d = 0; d < dE; ++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)PetscRealPart(f1[(q… in PetscFEIntegrateResidual_Basic()
512 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateBdResidual_Basic() local
547 dE = fgeom->dimEmbed; in PetscFEIntegrateBdResidual_Basic()
557 PetscCall(PetscArrayzero(f1, Nq * NcI * dE)); in PetscFEIntegrateBdResidual_Basic()
577 …for (i = 0; i < n0; ++i) f0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateBdResidual_Basic()
579 …_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom… in PetscFEIntegrateBdResidual_Basic()
581 for (d = 0; d < dE; ++d) f1[(q * NcI + c) * dE + d] *= w; in PetscFEIntegrateBdResidual_Basic()
627 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateHybridResidual_Basic() local
673 dE = fgeom->dimEmbed; in PetscFEIntegrateHybridResidual_Basic()
684 PetscCall(PetscArrayzero(f1, Nq * NcS * dE)); in PetscFEIntegrateHybridResidual_Basic()
703 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dE, fegeom.invJ)); in PetscFEIntegrateHybridResidual_Basic()
710 …for (i = 0; i < n0; ++i) f0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateHybridResidual_Basic()
712 …_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom… in PetscFEIntegrateHybridResidual_Basic()
714 for (d = 0; d < dE; ++d) f1[(q * NcS + c) * dE + d] *= w; in PetscFEIntegrateHybridResidual_Basic()
725 …d < dE; ++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, " f1[%" PetscInt_FMT ",%" PetscInt_FMT "] %g"… in PetscFEIntegrateHybridResidual_Basic()
763 PetscInt dE, Np; in PetscFEIntegrateJacobian_Basic() local
818 dE = cgeom->dimEmbed; 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()
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()
854 …for (PetscInt i = 0; i < n0; ++i) g0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_… in PetscFEIntegrateJacobian_Basic()
858 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateJacobian_Basic()
859 …for (PetscInt i = 0; i < n1; ++i) g1_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_… in PetscFEIntegrateJacobian_Basic()
860 for (PetscInt c = 0; c < NcI * NcJ * dE; ++c) g1[c] *= w; in PetscFEIntegrateJacobian_Basic()
863 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateJacobian_Basic()
864 …for (PetscInt i = 0; i < n2; ++i) g2_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_… in PetscFEIntegrateJacobian_Basic()
865 for (PetscInt c = 0; c < NcI * NcJ * dE; ++c) g2[c] *= w; in PetscFEIntegrateJacobian_Basic()
868 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); in PetscFEIntegrateJacobian_Basic()
869 …for (PetscInt i = 0; i < n3; ++i) g3_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_… in PetscFEIntegrateJacobian_Basic()
870 for (PetscInt c = 0; c < NcI * NcJ * dE * dE; ++c) g3[c] *= w; in PetscFEIntegrateJacobian_Basic()
914 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBdJacobian_Basic() local
957 dE = fgeom->dimEmbed; in PetscFEIntegrateBdJacobian_Basic()
961 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
962 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
963 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); 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()
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()
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()
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()
1014 …for (i = 0; i < n0; ++i) g0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateBdJacobian_Basic()
1018 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
1019 …for (i = 0; i < n1; ++i) g1_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateBdJacobian_Basic()
1023 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
1024 …for (i = 0; i < n2; ++i) g2_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateBdJacobian_Basic()
1028 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); in PetscFEIntegrateBdJacobian_Basic()
1029 …for (i = 0; i < n3; ++i) g3_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateBdJacobian_Basic()
1082 PetscInt qNc, Nq, q, dE; in PetscFEIntegrateHybridJacobian_Basic() local
1131 dE = fgeom->dimEmbed; in PetscFEIntegrateHybridJacobian_Basic()
1182 …for (i = 0; i < n0; ++i) g0_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateHybridJacobian_Basic()
1187 …for (i = 0; i < n1; ++i) g1_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateHybridJacobian_Basic()
1192 …for (i = 0; i < n2; ++i) g2_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateHybridJacobian_Basic()
1197 …for (i = 0; i < n3; ++i) g3_func[i](dE, Nf, NfAux, uOff, uOff_x, u, u_t, u_x, aOff, aOff_x, a, NUL… in PetscFEIntegrateHybridJacobian_Basic()