5574   IS              chunkISF, chunkISN;  in DMPlexComputeResidualHybridByKey()  local
5707 …Call(ISCreateGeneral(PETSC_COMM_SELF, 2 * cellChunkSize, neighbors, PETSC_USE_POINTER, &chunkISN));  in DMPlexComputeResidualHybridByKey()
5735     PetscCall(ISGeneralSetIndices(chunkISN, 2 * cellChunkSize, neighbors, PETSC_USE_POINTER));  in DMPlexComputeResidualHybridByKey()
5743         PetscCall(DMFieldCreateDefaultFaceQuadrature(coordField, chunkISN, &affineQuadN));  in DMPlexComputeResidualHybridByKey()
5746 …if (affineQuadN) PetscCall(DMSNESGetFEGeom(coordField, chunkISN, affineQuadN, PETSC_FEGEOM_BASIC, …  in DMPlexComputeResidualHybridByKey()
5750 …if (quadsN[f]) PetscCall(DMSNESGetFEGeom(coordField, chunkISN, quadsN[f], PETSC_FEGEOM_BASIC, &geo…  in DMPlexComputeResidualHybridByKey()
5843   PetscCall(ISDestroy(&chunkISN));  in DMPlexComputeResidualHybridByKey()
6655   IS              chunkISF, chunkISN;  in DMPlexComputeJacobianHybridByKey()  local
6778 …Call(ISCreateGeneral(PETSC_COMM_SELF, 2 * cellChunkSize, neighbors, PETSC_USE_POINTER, &chunkISN));  in DMPlexComputeJacobianHybridByKey()
6816     PetscCall(ISGeneralSetIndices(chunkISN, 2 * cellChunkSize, neighbors, PETSC_USE_POINTER));  in DMPlexComputeJacobianHybridByKey()
6823         PetscCall(DMFieldCreateDefaultFaceQuadrature(coordField, chunkISN, &affineQuadN));  in DMPlexComputeJacobianHybridByKey()
6826 …if (affineQuadN) PetscCall(DMSNESGetFEGeom(coordField, chunkISN, affineQuadN, PETSC_FEGEOM_BASIC, …  in DMPlexComputeJacobianHybridByKey()
6957   PetscCall(ISDestroy(&chunkISN));  in DMPlexComputeJacobianHybridByKey()