Lines Matching refs:ds
162 PETSC_INTERN PetscErrorCode PetscFEIntegrate_Basic(PetscDS ds, PetscInt field, PetscInt Ne, PetscFE… in PetscFEIntegrate_Basic() argument
164 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrate_Basic()
179 PetscCall(PetscDSGetObjective(ds, field, &obj_func)); in PetscFEIntegrate_Basic()
181 PetscCall(PetscDSGetDiscretization(ds, field, (PetscObject *)&fe)); in PetscFEIntegrate_Basic()
185 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrate_Basic()
186 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrate_Basic()
187 PetscCall(PetscDSGetComponentOffsets(ds, &uOff)); in PetscFEIntegrate_Basic()
188 PetscCall(PetscDSGetComponentDerivativeOffsets(ds, &uOff_x)); in PetscFEIntegrate_Basic()
189 PetscCall(PetscDSGetTabulation(ds, &T)); in PetscFEIntegrate_Basic()
190 PetscCall(PetscDSGetEvaluationArrays(ds, &u, NULL, &u_x)); in PetscFEIntegrate_Basic()
191 PetscCall(PetscDSGetWorkspace(ds, &x, NULL, NULL, NULL, NULL)); in PetscFEIntegrate_Basic()
192 PetscCall(PetscDSSetIntegrationParameters(ds, field, PETSC_DETERMINE)); in PetscFEIntegrate_Basic()
193 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrate_Basic()
233 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrate_Basic()
242 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, 0, q, T, &fegeom, &coefficients[cOffset], NULL… in PetscFEIntegrate_Basic()
255 PETSC_INTERN PetscErrorCode PetscFEIntegrateBd_Basic(PetscDS ds, PetscInt field, PetscBdPointFn *ob… in PetscFEIntegrateBd_Basic() argument
257 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateBd_Basic()
272 PetscCall(PetscDSGetDiscretization(ds, field, (PetscObject *)&fe)); in PetscFEIntegrateBd_Basic()
276 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateBd_Basic()
277 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateBd_Basic()
278 PetscCall(PetscDSGetComponentOffsets(ds, &uOff)); in PetscFEIntegrateBd_Basic()
279 PetscCall(PetscDSGetComponentDerivativeOffsets(ds, &uOff_x)); in PetscFEIntegrateBd_Basic()
280 PetscCall(PetscDSGetEvaluationArrays(ds, &u, NULL, &u_x)); in PetscFEIntegrateBd_Basic()
281 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, NULL, NULL)); in PetscFEIntegrateBd_Basic()
282 PetscCall(PetscDSGetFaceTabulation(ds, &Tf)); in PetscFEIntegrateBd_Basic()
283 PetscCall(PetscDSSetIntegrationParameters(ds, field, PETSC_DETERMINE)); in PetscFEIntegrateBd_Basic()
284 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateBd_Basic()
345 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBd_Basic()
376 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, Tf, &cgeom, &coefficients[cOffset], N… in PetscFEIntegrateBd_Basic()
389 PetscErrorCode PetscFEIntegrateResidual_Basic(PetscDS ds, PetscFormKey key, PetscInt Ne, PetscFEGeo… in PetscFEIntegrateResidual_Basic() argument
391 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateResidual_Basic()
408 PetscCall(PetscDSGetDiscretization(ds, field, (PetscObject *)&fe)); in PetscFEIntegrateResidual_Basic()
412 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateResidual_Basic()
413 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateResidual_Basic()
414 PetscCall(PetscDSGetComponentOffsets(ds, &uOff)); in PetscFEIntegrateResidual_Basic()
415 PetscCall(PetscDSGetComponentDerivativeOffsets(ds, &uOff_x)); in PetscFEIntegrateResidual_Basic()
416 PetscCall(PetscDSGetFieldOffset(ds, field, &fOffset)); in PetscFEIntegrateResidual_Basic()
417 PetscCall(PetscDSGetWeakForm(ds, &wf)); in PetscFEIntegrateResidual_Basic()
420 PetscCall(PetscDSGetEvaluationArrays(ds, &u, coefficients_t ? &u_t : NULL, &u_x)); in PetscFEIntegrateResidual_Basic()
421 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, NULL, NULL)); in PetscFEIntegrateResidual_Basic()
422 PetscCall(PetscDSGetWeakFormArrays(ds, &f0, &f1, NULL, NULL, NULL, NULL)); in PetscFEIntegrateResidual_Basic()
423 PetscCall(PetscDSGetTabulation(ds, &T)); in PetscFEIntegrateResidual_Basic()
424 PetscCall(PetscDSSetIntegrationParameters(ds, field, PETSC_DETERMINE)); in PetscFEIntegrateResidual_Basic()
425 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateResidual_Basic()
450 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateResidual_Basic()
458 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, 0, q, T, &fegeom, &coefficients[cOffset], Pets… in PetscFEIntegrateResidual_Basic()
496 PetscErrorCode PetscFEIntegrateBdResidual_Basic(PetscDS ds, PetscWeakForm wf, PetscFormKey key, Pet… in PetscFEIntegrateBdResidual_Basic() argument
498 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateBdResidual_Basic()
515 PetscCall(PetscDSGetDiscretization(ds, field, (PetscObject *)&fe)); in PetscFEIntegrateBdResidual_Basic()
519 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateBdResidual_Basic()
520 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateBdResidual_Basic()
521 PetscCall(PetscDSGetComponentOffsets(ds, &uOff)); in PetscFEIntegrateBdResidual_Basic()
522 PetscCall(PetscDSGetComponentDerivativeOffsets(ds, &uOff_x)); in PetscFEIntegrateBdResidual_Basic()
523 PetscCall(PetscDSGetFieldOffset(ds, field, &fOffset)); in PetscFEIntegrateBdResidual_Basic()
526 PetscCall(PetscDSGetEvaluationArrays(ds, &u, coefficients_t ? &u_t : NULL, &u_x)); in PetscFEIntegrateBdResidual_Basic()
527 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, NULL, NULL)); in PetscFEIntegrateBdResidual_Basic()
528 PetscCall(PetscDSGetWeakFormArrays(ds, &f0, &f1, NULL, NULL, NULL, NULL)); in PetscFEIntegrateBdResidual_Basic()
529 PetscCall(PetscDSGetFaceTabulation(ds, &Tf)); in PetscFEIntegrateBdResidual_Basic()
530 PetscCall(PetscDSSetIntegrationParameters(ds, field, PETSC_DETERMINE)); in PetscFEIntegrateBdResidual_Basic()
531 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateBdResidual_Basic()
564 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBdResidual_Basic()
575 …PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, Tf, &cgeom, &coefficients[cOffset], P… in PetscFEIntegrateBdResidual_Basic()
609 PETSC_INTERN PetscErrorCode PetscFEIntegrateHybridResidual_Basic(PetscDS ds, PetscDS dsIn, PetscFor… in PetscFEIntegrateHybridResidual_Basic() argument
611 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateHybridResidual_Basic()
631 PetscCall(PetscDSGetDiscretization(ds, field, (PetscObject *)&fe)); in PetscFEIntegrateHybridResidual_Basic()
634 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateHybridResidual_Basic()
635 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateHybridResidual_Basic()
639 PetscCall(PetscDSGetFieldOffsetCohesive(ds, field, &fOffset)); in PetscFEIntegrateHybridResidual_Basic()
640 PetscCall(PetscDSGetWeakForm(ds, &wf)); in PetscFEIntegrateHybridResidual_Basic()
643 PetscCall(PetscDSGetEvaluationArrays(ds, &u, coefficients_t ? &u_t : NULL, &u_x)); in PetscFEIntegrateHybridResidual_Basic()
644 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, NULL, NULL)); in PetscFEIntegrateHybridResidual_Basic()
645 PetscCall(PetscDSGetWeakFormArrays(ds, &f0, &f1, NULL, NULL, NULL, NULL)); in PetscFEIntegrateHybridResidual_Basic()
647 PetscCall(PetscDSGetTabulation(ds, &Tf)); in PetscFEIntegrateHybridResidual_Basic()
649 PetscCall(PetscDSSetIntegrationParameters(ds, field, PETSC_DETERMINE)); in PetscFEIntegrateHybridResidual_Basic()
650 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateHybridResidual_Basic()
663 PetscCall(PetscDSGetCohesive(ds, field, &isCohesiveField)); in PetscFEIntegrateHybridResidual_Basic()
694 …PetscCall(PetscDSPermuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[0], ornt[0]), … in PetscFEIntegrateHybridResidual_Basic()
695 …PetscCall(PetscDSPermuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[1], ornt[1]), … in PetscFEIntegrateHybridResidual_Basic()
893 PETSC_INTERN PetscErrorCode PetscFEIntegrateBdJacobian_Basic(PetscDS ds, PetscWeakForm wf, PetscFEJ… in PetscFEIntegrateBdJacobian_Basic() argument
895 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateBdJacobian_Basic()
917 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateBdJacobian_Basic()
920 PetscCall(PetscDSGetDiscretization(ds, fieldI, (PetscObject *)&feI)); in PetscFEIntegrateBdJacobian_Basic()
921 PetscCall(PetscDSGetDiscretization(ds, fieldJ, (PetscObject *)&feJ)); in PetscFEIntegrateBdJacobian_Basic()
925 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateBdJacobian_Basic()
926 PetscCall(PetscDSGetComponentOffsets(ds, &uOff)); in PetscFEIntegrateBdJacobian_Basic()
927 PetscCall(PetscDSGetComponentDerivativeOffsets(ds, &uOff_x)); in PetscFEIntegrateBdJacobian_Basic()
928 PetscCall(PetscDSGetFieldOffset(ds, fieldI, &offsetI)); in PetscFEIntegrateBdJacobian_Basic()
929 PetscCall(PetscDSGetFieldOffset(ds, fieldJ, &offsetJ)); in PetscFEIntegrateBdJacobian_Basic()
941 PetscCall(PetscDSGetEvaluationArrays(ds, &u, coefficients_t ? &u_t : NULL, &u_x)); in PetscFEIntegrateBdJacobian_Basic()
942 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, &testReal, &testDerReal)); in PetscFEIntegrateBdJacobian_Basic()
943 PetscCall(PetscDSGetWeakFormArrays(ds, NULL, NULL, &g0, &g1, &g2, &g3)); in PetscFEIntegrateBdJacobian_Basic()
944 PetscCall(PetscDSGetFaceTabulation(ds, &T)); in PetscFEIntegrateBdJacobian_Basic()
945 PetscCall(PetscDSSetIntegrationParameters(ds, fieldI, fieldJ)); in PetscFEIntegrateBdJacobian_Basic()
946 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateBdJacobian_Basic()
1008 PetscCall(PetscDSSetCellParameters(ds, fegeom.detJ[0] * cellScale)); in PetscFEIntegrateBdJacobian_Basic()
1010 …if (coefficients) PetscCall(PetscFEEvaluateFieldJets_Internal(ds, Nf, face, q, T, &cgeom, &coeffic… in PetscFEIntegrateBdJacobian_Basic()
1059 PETSC_INTERN PetscErrorCode PetscFEIntegrateHybridJacobian_Basic(PetscDS ds, PetscDS dsIn, PetscFEJ… in PetscFEIntegrateHybridJacobian_Basic() argument
1061 const PetscInt debug = ds->printIntegrate; in PetscFEIntegrateHybridJacobian_Basic()
1085 PetscCall(PetscDSGetNumFields(ds, &Nf)); in PetscFEIntegrateHybridJacobian_Basic()
1089 PetscCall(PetscDSGetDiscretization(ds, fieldI, (PetscObject *)&feI)); in PetscFEIntegrateHybridJacobian_Basic()
1090 PetscCall(PetscDSGetDiscretization(ds, fieldJ, (PetscObject *)&feJ)); in PetscFEIntegrateHybridJacobian_Basic()
1093 PetscCall(PetscDSGetTotalDimension(ds, &totDim)); in PetscFEIntegrateHybridJacobian_Basic()
1094 …PetscCall(PetscDSGetComponentOffsetsCohesive(ds, 0, &uOff)); // Change 0 to s for one-sided offsets in PetscFEIntegrateHybridJacobian_Basic()
1095 PetscCall(PetscDSGetComponentDerivativeOffsetsCohesive(ds, s, &uOff_x)); in PetscFEIntegrateHybridJacobian_Basic()
1096 PetscCall(PetscDSGetWeakForm(ds, &wf)); in PetscFEIntegrateHybridJacobian_Basic()
1108 PetscCall(PetscDSGetEvaluationArrays(ds, &u, coefficients_t ? &u_t : NULL, &u_x)); in PetscFEIntegrateHybridJacobian_Basic()
1109 PetscCall(PetscDSGetWorkspace(ds, &x, &basisReal, &basisDerReal, &testReal, &testDerReal)); in PetscFEIntegrateHybridJacobian_Basic()
1110 PetscCall(PetscDSGetWeakFormArrays(ds, NULL, NULL, &g0, &g1, &g2, &g3)); in PetscFEIntegrateHybridJacobian_Basic()
1111 PetscCall(PetscDSGetTabulation(ds, &T)); in PetscFEIntegrateHybridJacobian_Basic()
1113 PetscCall(PetscDSGetFieldOffsetCohesive(ds, fieldI, &offsetI)); in PetscFEIntegrateHybridJacobian_Basic()
1114 PetscCall(PetscDSGetFieldOffsetCohesive(ds, fieldJ, &offsetJ)); in PetscFEIntegrateHybridJacobian_Basic()
1115 PetscCall(PetscDSSetIntegrationParameters(ds, fieldI, fieldJ)); in PetscFEIntegrateHybridJacobian_Basic()
1116 PetscCall(PetscDSGetConstants(ds, &numConstants, &constants)); in PetscFEIntegrateHybridJacobian_Basic()
1129 PetscCall(PetscDSGetCohesive(ds, fieldI, &isCohesiveFieldI)); in PetscFEIntegrateHybridJacobian_Basic()
1130 PetscCall(PetscDSGetCohesive(ds, fieldJ, &isCohesiveFieldJ)); in PetscFEIntegrateHybridJacobian_Basic()
1165 …PetscCall(PetscDSPermuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, cornt[0], ornt[0]), … in PetscFEIntegrateHybridJacobian_Basic()
1166 …PetscCall(PetscDSPermuteQuadPoint(ds, DMPolytopeTypeComposeOrientationInv(ct, ornt[1], cornt[1]), … in PetscFEIntegrateHybridJacobian_Basic()