1034 static PetscErrorCode LAPACKsyevFail(PetscInt dim, PetscScalar Mpos[])  in LAPACKsyevFail()  argument
1044 … 1) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%15.8e, ", (double)PetscAbsScalar(Mpos[i * dim + j])));  in LAPACKsyevFail()
1045 … else PetscCall(PetscPrintf(PETSC_COMM_SELF, "%15.8e", (double)PetscAbsScalar(Mpos[i * dim + j])));  in LAPACKsyevFail()
1057   PetscScalar *Mpos;  in DMPlexMetricModify_Private()  local
1060   PetscCall(PetscMalloc2(dim * dim, &Mpos, dim, &eigs));  in DMPlexMetricModify_Private()
1064     Mpos[i * dim + i] = Mp[i * dim + i];  in DMPlexMetricModify_Private()
1066       Mpos[i * dim + j] = 0.5 * (Mp[i * dim + j] + Mp[j * dim + i]);  in DMPlexMetricModify_Private()
1067       Mpos[j * dim + i] = Mpos[i * dim + j];  in DMPlexMetricModify_Private()
1074     eigs[0] = PetscRealPart(Mpos[0]);  in DMPlexMetricModify_Private()
1075     Mpos[0] = 1.0;  in DMPlexMetricModify_Private()
1093 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, Mpos, &nb, eigs, work, &lwork, rwork, &lier…  in DMPlexMetricModify_Private()
1097 …    PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, Mpos, &nb, eigs, work, &lwork, &lierr));  in DMPlexMetricModify_Private()
1101           Mpos[i * dim + i] = Mp[i * dim + i];  in DMPlexMetricModify_Private()
1103             Mpos[i * dim + j] = 0.5 * (Mp[i * dim + j] + Mp[j * dim + i]);  in DMPlexMetricModify_Private()
1104             Mpos[j * dim + i] = Mpos[i * dim + j];  in DMPlexMetricModify_Private()
1107         PetscCall(LAPACKsyevFail(dim, Mpos));  in DMPlexMetricModify_Private()
1133       for (k = 0; k < dim; ++k) Mp[i * dim + j] += Mpos[k * dim + i] * eigs[k] * Mpos[k * dim + j];  in DMPlexMetricModify_Private()
1136   PetscCall(PetscFree2(Mpos, eigs));  in DMPlexMetricModify_Private()