1478     PetscScalar   *elemMatrix, *elMat;  in CreateStaticData()  local
1522             PetscScalar *valuesOrig = elMat = elemMatrix;  in CreateStaticData()
1523             PetscCall(PetscArrayzero(elMat, totDim * totDim));  in CreateStaticData()
1524             elMat[(fieldA * Nb + q) * totDim + fieldA * Nb + q] = 1;  in CreateStaticData()
1525 …lex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numindices, &indices, NULL, &elMat));  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()
1551 …     pointMaps[maps[grid].num_reduced][jj].scale += PetscRealPart(elMat[f * numindices + ff + ii]);  in CreateStaticData()
1578 …                  for (f = 0; f < numindices; ++f) tmp += PetscRealPart(elMat[d * numindices + f]);  in CreateStaticData()
1590 …lex[grid], section[grid], globsection[grid], ej, PETSC_TRUE, &numindices, &indices, NULL, &elMat));  in CreateStaticData()
1591 …if (elMat != valuesOrig) PetscCall(DMRestoreWorkArray(ctx->plex[grid], numindices * numindices, MP…  in CreateStaticData()