Lines Matching refs:cellIS
2423 IS cellIS; in DMPlexComputeIntegral_Internal() local
2443 PetscCall(ISCreateStride(PETSC_COMM_SELF, numCells, cStart, 1, &cellIS)); in DMPlexComputeIntegral_Internal()
2461 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeIntegral_Internal()
2463 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &affineQuad)); in DMPlexComputeIntegral_Internal()
2464 …if (affineQuad) PetscCall(DMFieldCreateFEGeom(coordField, cellIS, affineQuad, PETSC_FEGEOM_BASIC, … in DMPlexComputeIntegral_Internal()
2541 …if (!affineQuad) PetscCall(DMFieldCreateFEGeom(coordField, cellIS, q, PETSC_FEGEOM_BASIC, &cgeomFE… in DMPlexComputeIntegral_Internal()
2581 PetscCall(ISDestroy(&cellIS)); in DMPlexComputeIntegral_Internal()
3751 PetscErrorCode DMPlexGetCellFields(DM dm, IS cellIS, Vec locX, PeOp Vec locX_t, PeOp Vec locA, Pets… in DMPlexGetCellFields() argument
3769 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetCellFields()
3814 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetCellFields()
3837 PetscErrorCode DMPlexRestoreCellFields(DM dm, IS cellIS, Vec locX, PeOp Vec locX_t, PeOp Vec locA, … in DMPlexRestoreCellFields() argument
3846 static PetscErrorCode DMPlexGetHybridCellFields(DM dm, IS cellIS, Vec locX, Vec locX_t, Vec locA, P… in DMPlexGetHybridCellFields() argument
3857 PetscValidHeaderSpecific(cellIS, IS_CLASSID, 2); in DMPlexGetHybridCellFields()
3864 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetHybridCellFields()
3971 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetHybridCellFields()
3996 static PetscErrorCode DMPlexGetHybridFields(DM dm, DM dmX[], PetscDS dsX[], IS cellIS, Vec locX[], … in DMPlexGetHybridFields() argument
4009 PetscValidHeaderSpecific(cellIS, IS_CLASSID, 4); in DMPlexGetHybridFields()
4011 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetHybridFields()
4055 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetHybridFields()
4059 static PetscErrorCode DMPlexRestoreHybridFields(DM dm, DM dmX[], PetscDS dsX[], IS cellIS, Vec locX… in DMPlexRestoreHybridFields() argument
4389 PetscErrorCode DMPlexComputeResidual_Patch_Internal(DM dm, PetscSection section, IS cellIS, PetscRe… in DMPlexComputeResidual_Patch_Internal() argument
4440 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeResidual_Patch_Internal()
4442 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &affineQuad)); in DMPlexComputeResidual_Patch_Internal()
4443 …if (affineQuad) PetscCall(DMSNESGetFEGeom(coordField, cellIS, affineQuad, PETSC_FEGEOM_BASIC, &aff… in DMPlexComputeResidual_Patch_Internal()
4460 PetscCall(DMSNESGetFEGeom(coordField, cellIS, quads[f], PETSC_FEGEOM_BASIC, &geoms[f])); in DMPlexComputeResidual_Patch_Internal()
4466 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidual_Patch_Internal()
4594 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidual_Patch_Internal()
4598 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuad, PETSC_FALSE, &affineGeom)); in DMPlexComputeResidual_Patch_Internal()
4602 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, quads[f], PETSC_FALSE, &geoms[f])); in DMPlexComputeResidual_Patch_Internal()
4618 …tch_Internal(DM dm, PetscSection section, PetscSection globalSection, IS cellIS, PetscReal t, Pets… in DMPlexComputeJacobian_Patch_Internal() argument
4637 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeJacobian_Patch_Internal()
4638 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobian_Patch_Internal()
4704 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeJacobian_Patch_Internal()
4705 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &qGeom)); in DMPlexComputeJacobian_Patch_Internal()
4713 PetscCall(DMSNESGetFEGeom(coordField, cellIS, qGeom, PETSC_FEGEOM_BASIC, &cgeomFEM)); in DMPlexComputeJacobian_Patch_Internal()
4788 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, qGeom, PETSC_FALSE, &cgeomFEM)); in DMPlexComputeJacobian_Patch_Internal()
5175 PetscErrorCode DMPlexComputeResidualByKey(DM dm, PetscFormKey key, IS cellIS, PetscReal time, Vec l… in DMPlexComputeResidualByKey() argument
5202 if (!cellIS) goto end; in DMPlexComputeResidualByKey()
5203 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidualByKey()
5241 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeResidualByKey()
5243 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &affineQuad)); in DMPlexComputeResidualByKey()
5244 …if (affineQuad) PetscCall(DMSNESGetFEGeom(coordField, cellIS, affineQuad, PETSC_FEGEOM_BASIC, &aff… in DMPlexComputeResidualByKey()
5261 PetscCall(DMSNESGetFEGeom(coordField, cellIS, quads[f], PETSC_FEGEOM_BASIC, &geoms[f])); in DMPlexComputeResidualByKey()
5461 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidualByKey()
5467 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuad, PETSC_FALSE, &affineGeom)); in DMPlexComputeResidualByKey()
5471 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, quads[f], PETSC_FALSE, &geoms[f])); in DMPlexComputeResidualByKey()
5555 PetscErrorCode DMPlexComputeResidualHybridByKey(DM dm, PetscFormKey key[], IS cellIS, PetscReal tim… in DMPlexComputeResidualHybridByKey() argument
5587 if (!cellIS) goto end; in DMPlexComputeResidualHybridByKey()
5588 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidualHybridByKey()
5589 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeResidualHybridByKey()
5675 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeResidualHybridByKey()
5710 PetscCall(DMPlexGetHybridCellFields(dm, cellIS, locX, locX_t, locA[2], &u, &u_t, &a[2])); in DMPlexComputeResidualHybridByKey()
5711 PetscCall(DMPlexGetHybridFields(dm, dmAux, dsAux, cellIS, locA, PETSC_TRUE, a)); in DMPlexComputeResidualHybridByKey()
5712 PetscCall(DMPlexGetHybridFields(dm, dmScale, dsScale, cellIS, locS, PETSC_TRUE, s)); in DMPlexComputeResidualHybridByKey()
5835 PetscCall(DMPlexRestoreCellFields(dm, cellIS, locX, locX_t, locA[2], &u, &u_t, &a[2])); in DMPlexComputeResidualHybridByKey()
5836 PetscCall(DMPlexRestoreHybridFields(dm, dmAux, dsAux, cellIS, locA, PETSC_TRUE, a)); in DMPlexComputeResidualHybridByKey()
5837 PetscCall(DMPlexRestoreHybridFields(dm, dmScale, dsScale, cellIS, locS, PETSC_TRUE, s)); in DMPlexComputeResidualHybridByKey()
5844 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeResidualHybridByKey()
5846 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuadF, PETSC_FALSE, &affineGeomF)); in DMPlexComputeResidualHybridByKey()
5848 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuadN, PETSC_FALSE, &affineGeomN)); in DMPlexComputeResidualHybridByKey()
5852 …if (geomsF) PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, quadsF[f], PETSC_FALSE, &geomsF[f])); in DMPlexComputeResidualHybridByKey()
5854 …if (geomsN) PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, quadsN[f], PETSC_FALSE, &geomsN[f])); in DMPlexComputeResidualHybridByKey()
6172 PetscErrorCode DMPlexComputeJacobianByKey(DM dm, PetscFormKey key, IS cellIS, PetscReal t, PetscRea… in DMPlexComputeJacobianByKey() argument
6202 if (!cellIS) goto end; in DMPlexComputeJacobianByKey()
6203 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianByKey()
6204 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeJacobianByKey()
6266 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeJacobianByKey()
6267 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &qGeom)); in DMPlexComputeJacobianByKey()
6273 PetscCall(DMSNESGetFEGeom(coordField, cellIS, qGeom, PETSC_FEGEOM_BASIC, &cgeomFEM)); in DMPlexComputeJacobianByKey()
6300 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, qGeom, PETSC_FALSE, &cgeomFEM)); in DMPlexComputeJacobianByKey()
6362 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianByKey()
6385 PetscErrorCode DMPlexComputeJacobianByKeyGeneral(DM dmr, DM dmc, PetscFormKey key, IS cellIS, Petsc… in DMPlexComputeJacobianByKeyGeneral() argument
6421 if (!cellIS) goto end; in DMPlexComputeJacobianByKeyGeneral()
6422 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianByKeyGeneral()
6423 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeJacobianByKeyGeneral()
6489 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeJacobianByKeyGeneral()
6490 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &qGeom)); in DMPlexComputeJacobianByKeyGeneral()
6496 PetscCall(DMSNESGetFEGeom(coordField, cellIS, qGeom, PETSC_FEGEOM_BASIC, &cgeomFEM)); in DMPlexComputeJacobianByKeyGeneral()
6523 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, qGeom, PETSC_FALSE, &cgeomFEM)); in DMPlexComputeJacobianByKeyGeneral()
6585 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianByKeyGeneral()
6631 PetscErrorCode DMPlexComputeJacobianHybridByKey(DM dm, PetscFormKey key[], IS cellIS, PetscReal t, … in DMPlexComputeJacobianHybridByKey() argument
6669 if (!cellIS) goto end; in DMPlexComputeJacobianHybridByKey()
6670 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianHybridByKey()
6671 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeJacobianHybridByKey()
6759 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeJacobianHybridByKey()
6781 PetscCall(DMPlexGetHybridCellFields(dm, cellIS, locX, locX_t, locA[2], &u, &u_t, &a[2])); in DMPlexComputeJacobianHybridByKey()
6782 PetscCall(DMPlexGetHybridFields(dm, dmAux, dsAux, cellIS, locA, PETSC_TRUE, a)); in DMPlexComputeJacobianHybridByKey()
6783 PetscCall(DMPlexGetHybridFields(dm, dmScale, dsScale, cellIS, locS, PETSC_TRUE, s)); in DMPlexComputeJacobianHybridByKey()
6947 PetscCall(DMPlexRestoreCellFields(dm, cellIS, locX, locX_t, locA[2], &u, &u_t, &a[2])); in DMPlexComputeJacobianHybridByKey()
6948 PetscCall(DMPlexRestoreHybridFields(dm, dmAux, dsAux, cellIS, locA, PETSC_TRUE, a)); in DMPlexComputeJacobianHybridByKey()
6958 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianHybridByKey()
6960 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuadF, PETSC_FALSE, &affineGeomF)); in DMPlexComputeJacobianHybridByKey()
6962 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, affineQuadN, PETSC_FALSE, &affineGeomN)); in DMPlexComputeJacobianHybridByKey()
6967 …if (geomsF) PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, quadsF[f], PETSC_FALSE, &geomsF[f])); in DMPlexComputeJacobianHybridByKey()
7002 PetscErrorCode DMPlexComputeJacobianActionByKey(DM dm, PetscFormKey key, IS cellIS, PetscReal t, Pe… in DMPlexComputeJacobianActionByKey() argument
7022 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMPlexComputeJacobianActionByKey()
7023 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianActionByKey()
7086 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeJacobianActionByKey()
7087 if (maxDegree <= 1) PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &qGeom)); in DMPlexComputeJacobianActionByKey()
7093 PetscCall(DMSNESGetFEGeom(coordField, cellIS, qGeom, PETSC_FEGEOM_BASIC, &cgeomFEM)); in DMPlexComputeJacobianActionByKey()
7114 PetscCall(DMSNESRestoreFEGeom(coordField, cellIS, qGeom, PETSC_FALSE, &cgeomFEM)); in DMPlexComputeJacobianActionByKey()
7141 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexComputeJacobianActionByKey()