5574   IS              chunkISF, chunkISN;  in DMPlexComputeResidualHybridByKey()  local
5706 …PetscCall(ISCreateGeneral(PETSC_COMM_SELF, 2 * cellChunkSize, faces, PETSC_USE_POINTER, &chunkISF)…  in DMPlexComputeResidualHybridByKey()
5734     PetscCall(ISGeneralSetIndices(chunkISF, 2 * cellChunkSize, faces, PETSC_USE_POINTER));  in DMPlexComputeResidualHybridByKey()
5738 …   if (!affineQuadF) PetscCall(DMFieldCreateDefaultQuadrature(coordField, chunkISF, &affineQuadF));  in DMPlexComputeResidualHybridByKey()
5739 …if (affineQuadF) PetscCall(DMSNESGetFEGeom(coordField, chunkISF, affineQuadF, PETSC_FEGEOM_COHESIV…  in DMPlexComputeResidualHybridByKey()
5749 …if (quadsF[f]) PetscCall(DMSNESGetFEGeom(coordField, chunkISF, quadsF[f], PETSC_FEGEOM_COHESIVE, &…  in DMPlexComputeResidualHybridByKey()
5842   PetscCall(ISDestroy(&chunkISF));  in DMPlexComputeResidualHybridByKey()
6655   IS              chunkISF, chunkISN;  in DMPlexComputeJacobianHybridByKey()  local
6777 …PetscCall(ISCreateGeneral(PETSC_COMM_SELF, 2 * cellChunkSize, faces, PETSC_USE_POINTER, &chunkISF)…  in DMPlexComputeJacobianHybridByKey()
6815     PetscCall(ISGeneralSetIndices(chunkISF, 2 * cellChunkSize, faces, PETSC_USE_POINTER));  in DMPlexComputeJacobianHybridByKey()
6818 …   if (!affineQuadF) PetscCall(DMFieldCreateDefaultQuadrature(coordField, chunkISF, &affineQuadF));  in DMPlexComputeJacobianHybridByKey()
6819 …if (affineQuadF) PetscCall(DMSNESGetFEGeom(coordField, chunkISF, affineQuadF, PETSC_FEGEOM_COHESIV…  in DMPlexComputeJacobianHybridByKey()
6830 …if (quadsF[f]) PetscCall(DMSNESGetFEGeom(coordField, chunkISF, quadsF[f], PETSC_FEGEOM_COHESIVE, &…  in DMPlexComputeJacobianHybridByKey()
6956   PetscCall(ISDestroy(&chunkISF));  in DMPlexComputeJacobianHybridByKey()