Lines Matching refs:eigs
1056 …PetscReal *eigs, max_eig, l_min = 1.0 / (h_max * h_max), l_max = 1.0 / (h_min * h_min), la_min =… in DMPlexMetricModify_Private() local
1060 PetscCall(PetscMalloc2(dim * dim, &Mpos, dim, &eigs)); in DMPlexMetricModify_Private()
1074 eigs[0] = PetscRealPart(Mpos[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()
1118 eigs[i] = PetscMin(l_max, PetscMax(l_min, PetscAbsReal(eigs[i]))); in DMPlexMetricModify_Private()
1119 max_eig = PetscMax(eigs[i], max_eig); in DMPlexMetricModify_Private()
1125 if (a_max >= 1.0) eigs[i] = PetscMax(eigs[i], max_eig * la_min); in DMPlexMetricModify_Private()
1126 *detMp *= eigs[i]; 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()