1516       for (PetscInt ej = cStart, eidx = 0; ej < cEnd; ++ej, ++eidx, glb_elem_idx++) {  in CreateStaticData()  local
1525 …PetscCall(DMPlexGetClosureIndices(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRU…  in CreateStaticData()
1590 …exRestoreClosureIndices(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numind…  in CreateStaticData()
1604       for (PetscInt ej = 0; ej < maps[grid].num_reduced; ++ej) {  in CreateStaticData()  local
1606           maps[grid].c_maps[ej][q].scale = pointMaps[ej][q].scale;  in CreateStaticData()
1607           maps[grid].c_maps[ej][q].gid   = pointMaps[ej][q].gid;  in CreateStaticData()
1626         for (PetscInt ej = 0; ej < numCells[grid]; ++ej, glb_elem_idx++) {  in CreateStaticData()  local
1628 …const LandauIdx *const Idxs             = &maps[grid].gIdx[ej][0][0]; // just use field-0 maps, Th…  in CreateStaticData()
1653           for (PetscInt ej = 0; ej < numCells[grid]; ++ej, glb_elem_idx++) {  in CreateStaticData()  local
1657               const LandauIdx *const Idxs = &maps[grid].gIdx[ej][fieldA][0];  in CreateStaticData()
1726       PetscInt     cStart, cEnd, ej;  in CreateStaticData()  local
1753       for (ej = 0; ej < numCells[grid]; ++ej, ++outer_ej) {  in CreateStaticData()
1757 …PetscCall(DMPlexComputeCellGeometryFEM(ctx->plex[grid], ej + cStart, quad, vj, Jdummy, invJe, detJ…  in CreateStaticData()
1758 …rgy_tensor_trick) PetscCall(DMPlexVecGetClosure(dmEnergy, e_section, v2_2, ej + cStart, NULL, &coe…  in CreateStaticData()
1784 …SqrtReal(xx[gidx] * xx[gidx] + yy[gidx] * yy[gidx] + zz[gidx] * zz[gidx]), ej, qj, dg2_c2, (double…  in CreateStaticData()
1805 …tensor_trick) PetscCall(DMPlexVecRestoreClosure(dmEnergy, e_section, v2_2, ej + cStart, NULL, &coe…  in CreateStaticData()