1470   PetscScalar *evecs, *sqrtM1, *isqrtM1;  in DMPlexMetricIntersection_Private()  local
1480   PetscCall(PetscMalloc4(dim * dim, &evecs, dim * dim, &sqrtM1, dim * dim, &isqrtM1, dim, &evals));  in DMPlexMetricIntersection_Private()
1482     for (j = 0; j < dim; ++j) evecs[i * dim + j] = M1[i * dim + j];  in DMPlexMetricIntersection_Private()
1501 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li…  in DMPlexMetricIntersection_Private()
1505 …  PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, &lierr));  in DMPlexMetricIntersection_Private()
1520             sqrtM1[i * dim + k] += evecs[j * dim + i] * sqrtj * evecs[j * dim + k];  in DMPlexMetricIntersection_Private()
1521             isqrtM1[i * dim + k] += evecs[j * dim + i] * (1.0 / sqrtj) * evecs[j * dim + k];  in DMPlexMetricIntersection_Private()
1529           evecs[i * dim + l] = 0.0;  in DMPlexMetricIntersection_Private()
1531 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k…  in DMPlexMetricIntersection_Private()
1542 …PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, rwork, &li…  in DMPlexMetricIntersection_Private()
1546 …  PetscCallBLAS("LAPACKsyev", LAPACKsyev_("V", "U", &nb, evecs, &nb, evals, work, &lwork, &lierr));  in DMPlexMetricIntersection_Private()
1551             evecs[i * dim + l] = 0.0;  in DMPlexMetricIntersection_Private()
1553 …for (k = 0; k < dim; ++k) evecs[i * dim + l] += isqrtM1[j * dim + i] * M2[j * dim + k] * isqrtM1[k…  in DMPlexMetricIntersection_Private()
1557         PetscCall(LAPACKsyevFail(dim, evecs));  in DMPlexMetricIntersection_Private()
1571 …0; l < dim; ++l) M2[i * dim + m] += sqrtM1[j * dim + i] * evecs[j * dim + k] * evals[k] * evecs[l …  in DMPlexMetricIntersection_Private()
1579   PetscCall(PetscFree4(evecs, sqrtM1, isqrtM1, evals));  in DMPlexMetricIntersection_Private()