1521             PetscInt     numindices, *indices;  in CreateStaticData()  local
1525 …es(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numindices, &indices, NULL,…  in CreateStaticData()
1527 …PetscCheck(numindices == Nb, ctx->comm, PETSC_ERR_ARG_WRONG, "numindices != Nb numindices=%" Petsc…  in CreateStaticData()
1528 …           for (PetscInt q = 0; q < numindices; ++q) maps[grid].gIdx[eidx][fieldA][q] = indices[q];  in CreateStaticData()
1531 …for (PetscInt f = 0; f < numindices; ++f) { // look for a non-zero on the diagonal (is this too co…  in CreateStaticData()
1532                 if (PetscAbs(PetscRealPart(elMat[f * numindices + f])) > PETSC_MACHINE_EPSILON) {  in CreateStaticData()
1534 …if (PetscAbs(PetscRealPart(elMat[f * numindices + f] - 1.)) < PETSC_MACHINE_EPSILON) { // normal v…  in CreateStaticData()
1550 …if (ff + ii < numindices) {                                                                   // 3…  in CreateStaticData()
1551 …     pointMaps[maps[grid].num_reduced][jj].scale += PetscRealPart(elMat[f * numindices + ff + ii]);  in CreateStaticData()
1565 …   } while (++jj < maps[grid].num_face && ++f < numindices); // jj is incremented if we hit the end  in CreateStaticData()
1576                       for (d = 0, tmp = 0; d < numindices; ++d) {  in CreateStaticData()
1578 …                  for (f = 0; f < numindices; ++f) tmp += PetscRealPart(elMat[d * numindices + f]);  in CreateStaticData()
1590 …es(ctx->plex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numindices, &indices, NULL,…  in CreateStaticData()
1591 …Mat != valuesOrig) PetscCall(DMRestoreWorkArray(ctx->plex[grid], numindices * numindices, MPIU_SCA…  in CreateStaticData()