Lines Matching refs:q
60 PetscInt Nc, Nq, q, k, c; in PetscFESetUp_Basic() local
70 for (q = 0; q < Nq; ++q) { in PetscFESetUp_Basic()
71 …for (c = 0; c < Nc; ++c) fem->invV[j * pdim + k] += Bf[(q * pdim + k) * Nc + c] * weights[q * Nc +… in PetscFESetUp_Basic()
176 PetscInt qNc, Nq, q; in PetscFEIntegrate_Basic() local
221 for (q = 0; q < Nq; ++q) { 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()
234 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrate_Basic()
235 if (debug > 1 && q < Np) { in PetscFEIntegrate_Basic()
241 if (debug) PetscCall(PetscPrintf(PETSC_COMM_SELF, " quad point %" PetscInt_FMT "\n", q)); 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()
267 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBd_Basic() local
328 for (q = 0; q < Nq; ++q) { in PetscFEIntegrateBd_Basic()
333 …RefToReal(dE, dim - 1, fegeom.xi, &fgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * (dim - 1)], x); 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()
346 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateBd_Basic()
347 if (debug > 1 && q < Np) { in PetscFEIntegrateBd_Basic()
353 if (debug > 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, " quad point %" PetscInt_FMT "\n", 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()
405 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateResidual_Basic() local
445 for (q = 0; q < Nq; ++q) { in PetscFEIntegrateResidual_Basic()
449 PetscCall(PetscFEGeomGetPoint(cgeom, e, q, &quadPoints[q * cgeom->dim], &fegeom)); in PetscFEIntegrateResidual_Basic()
451 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateResidual_Basic()
452 if (debug > 1 && q < cgeom->numPoints) { 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 …u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, numConstants, constants, &f0[q * T[field]->Nc]); in PetscFEIntegrateResidual_Basic()
461 for (c = 0; c < T[field]->Nc; ++c) f0[q * T[field]->Nc + c] *= w; in PetscFEIntegrateResidual_Basic()
462 … u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, numConstants, constants, &f1[q * T[field]->Nc * dE]… in PetscFEIntegrateResidual_Basic()
464 for (d = 0; d < dE; ++d) f1[(q * T[field]->Nc + c) * dE + d] *= w; in PetscFEIntegrateResidual_Basic()
467 …etscPrintf(PETSC_COMM_SELF, " quad point %" PetscInt_FMT " wt %g x:", q, (double)quadWeights[q])); in PetscFEIntegrateResidual_Basic()
478 …++c) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)PetscRealPart(f0[q * T[field]->Nc + c])… in PetscFEIntegrateResidual_Basic()
482 …+d) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)PetscRealPart(f1[(q * T[field]->Nc + c) … in PetscFEIntegrateResidual_Basic()
512 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateBdResidual_Basic() local
558 for (q = 0; q < Nq; ++q) { 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()
565 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateBdResidual_Basic()
567 if ((fgeom->isAffine && q == 0) || !fgeom->isAffine) { in PetscFEIntegrateBdResidual_Basic()
575 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, Tf, &cgeom, &coefficients[cOffset], P… in PetscFEIntegrateBdResidual_Basic()
576 …if (dsAux) PetscCall(PetscFEEvaluateFieldJets_Internal(dsAux, NfAux, auxOnBd ? 0 : face, q, TfAux,… in PetscFEIntegrateBdResidual_Basic()
577 …_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f0[q * NcI]); in PetscFEIntegrateBdResidual_Basic()
578 for (c = 0; c < NcI; ++c) f0[q * NcI + c] *= w; in PetscFEIntegrateBdResidual_Basic()
579 …_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f1[q * NcI * dE]); in PetscFEIntegrateBdResidual_Basic()
581 for (d = 0; d < dE; ++d) f1[(q * NcI + c) * dE + d] *= w; in PetscFEIntegrateBdResidual_Basic()
583 …all(PetscPrintf(PETSC_COMM_SELF, " elem %" PetscInt_FMT " quad point %" PetscInt_FMT "\n", e, q)); in PetscFEIntegrateBdResidual_Basic()
585 …rintf(PETSC_COMM_SELF, " f0[%" PetscInt_FMT "] %g\n", c, (double)PetscRealPart(f0[q * NcI + c]))); in PetscFEIntegrateBdResidual_Basic()
587 …1[%" PetscInt_FMT ",%" PetscInt_FMT "] %g", c, d, (double)PetscRealPart(f1[(q * NcI + c) * dim + d… in PetscFEIntegrateBdResidual_Basic()
627 PetscInt qdim, qNc, Nq, q, dE; in PetscFEIntegrateHybridResidual_Basic() local
689 for (q = 0; q < Nq; ++q) { in PetscFEIntegrateHybridResidual_Basic()
694 …rmuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[0], ornt[0]), field, q, &qpt[0])); in PetscFEIntegrateHybridResidual_Basic()
695 …rmuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[1], ornt[1]), field, q, &qpt[1])); in PetscFEIntegrateHybridResidual_Basic()
696 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridResidual_Basic()
697 PetscCall(PetscFEGeomGetPoint(nbrgeom, e * 2, q, NULL, &fegeomN[0])); in PetscFEIntegrateHybridResidual_Basic()
698 PetscCall(PetscFEGeomGetPoint(nbrgeom, e * 2 + 1, q, NULL, &fegeomN[1])); in PetscFEIntegrateHybridResidual_Basic()
699 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateHybridResidual_Basic()
700 if (debug > 1 && q < fgeom->numPoints) { in PetscFEIntegrateHybridResidual_Basic()
706 …OMM_SELF, " quad point %" PetscInt_FMT " weight %g detJ %g\n", q, (double)quadWeights[q], (double… in PetscFEIntegrateHybridResidual_Basic()
708 …PetscCall(PetscFEEvaluateFieldJets_Hybrid_Internal(dsIn, Nf, 0, q, Tf, face, qpt, TfIn, &fegeom, f… in PetscFEIntegrateHybridResidual_Basic()
709 …FEEvaluateFieldJets_Internal(dsAux, NfAux, auxOnBd ? 0 : face[s], auxOnBd ? q : qpt[s], TfAux, &fe… in PetscFEIntegrateHybridResidual_Basic()
710 …_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f0[q * NcS]); in PetscFEIntegrateHybridResidual_Basic()
711 for (c = 0; c < NcS; ++c) f0[q * NcS + c] *= w; in PetscFEIntegrateHybridResidual_Basic()
712 …_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f1[q * NcS * dE]); in PetscFEIntegrateHybridResidual_Basic()
714 for (d = 0; d < dE; ++d) f1[(q * NcS + c) * dE + d] *= w; in PetscFEIntegrateHybridResidual_Basic()
716 …" quad point %" PetscInt_FMT " field %" PetscInt_FMT " side %" PetscInt_FMT "\n", e, q, field, s)); in PetscFEIntegrateHybridResidual_Basic()
723 …rintf(PETSC_COMM_SELF, " f0[%" PetscInt_FMT "] %g\n", c, (double)PetscRealPart(f0[q * NcS + c]))); in PetscFEIntegrateHybridResidual_Basic()
725 …1[%" PetscInt_FMT ",%" PetscInt_FMT "] %g", c, d, (double)PetscRealPart(f1[(q * NcS + c) * dE + d]… in PetscFEIntegrateHybridResidual_Basic()
836 for (PetscInt q = 0; q < Nq; ++q) { in PetscFEIntegrateJacobian_Basic() local
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()
848 …OMM_SELF, " quad point %" PetscInt_FMT " weight %g detJ %g\n", q, (double)quadWeights[q], (double… 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()
873 …PetscCall(PetscFEUpdateElementMat_Internal(feI, feJ, 0, q, rT[fieldI], basisReal, basisDerReal, cT… in PetscFEIntegrateJacobian_Basic()
914 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBdJacobian_Basic() local
990 for (q = 0; q < Nq; ++q) { in PetscFEIntegrateBdJacobian_Basic()
994 if (debug) PetscCall(PetscPrintf(PETSC_COMM_SELF, " quad point %" PetscInt_FMT "\n", q)); in PetscFEIntegrateBdJacobian_Basic()
996 …RefToReal(dE, dim - 1, fegeom.xi, &fgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * (dim - 1)], x); 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()
1009 w = fegeom.detJ[0] * quadWeights[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 …PetscCall(PetscFEUpdateElementMat_Internal(feI, feJ, face, q, T[fieldI], basisReal, basisDerReal, … in PetscFEIntegrateBdJacobian_Basic()
1082 PetscInt qNc, Nq, q, dE; in PetscFEIntegrateHybridJacobian_Basic() local
1160 for (q = 0; q < Nq; ++q) { in PetscFEIntegrateHybridJacobian_Basic()
1165 …muteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[0], ornt[0]), fieldI, q, &qpt[0])); in PetscFEIntegrateHybridJacobian_Basic()
1166 …muteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, ornt[1], cornt[1]), fieldI, q, &qpt[1])); in PetscFEIntegrateHybridJacobian_Basic()
1167 PetscCall(PetscFEGeomGetPoint(fgeom, e * 2, q, &quadPoints[q * fgeom->dim], &fegeom)); in PetscFEIntegrateHybridJacobian_Basic()
1168 PetscCall(PetscFEGeomGetPoint(nbrgeom, e * 2, q, NULL, &fegeomN[0])); in PetscFEIntegrateHybridJacobian_Basic()
1169 PetscCall(PetscFEGeomGetPoint(nbrgeom, e * 2 + 1, q, NULL, &fegeomN[1])); in PetscFEIntegrateHybridJacobian_Basic()
1170 w = fegeom.detJ[0] * quadWeights[q]; in PetscFEIntegrateHybridJacobian_Basic()
1171 if (debug > 1 && q < fgeom->numPoints) { in PetscFEIntegrateHybridJacobian_Basic()
1177 if (debug) PetscCall(PetscPrintf(PETSC_COMM_SELF, " quad point %" PetscInt_FMT "\n", q)); in PetscFEIntegrateHybridJacobian_Basic()
1178 …if (coefficients) PetscCall(PetscFEEvaluateFieldJets_Hybrid_Internal(dsIn, Nf, 0, q, T, face, qpt,… in PetscFEIntegrateHybridJacobian_Basic()
1179 …FEEvaluateFieldJets_Internal(dsAux, NfAux, auxOnBd ? 0 : face[s], auxOnBd ? q : qpt[s], TAux, &feg… in PetscFEIntegrateHybridJacobian_Basic()
1204 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 0, 0, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1206 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 0, 0, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1207 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 1, 1, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1212 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 0, 0, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1213 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 1, 1, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1215 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 0, 0, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1216 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 0, 1, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1217 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 1, 0, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1218 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, 1, 1, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()
1221 …tMat_Hybrid_Internal(feI, isCohesiveFieldI, feJ, isCohesiveFieldJ, 0, s, s, q, T[fieldI], basisRea… in PetscFEIntegrateHybridJacobian_Basic()