| /petsc/src/mat/tests/ |
| H A D | ex116.c | 23 PetscScalar *arrayA, *evecs_array, *work, *evals; in main() local
67 PetscCall(PetscMalloc1(n, &evals)); in main()
73 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, &lierr); in main()
90 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
100 …all(PetscPrintf(PETSC_COMM_SELF, "%" PetscInt_FMT " %g\n", (PetscInt)(i + il), (double)evals[i])); in main()
113 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
153 …LAPACKgesvd_("S", "S", &im, &in, arrayA, &im, evals, arrayU, &minMN, arrayVT, &minMN, work, &lwork… in main()
157 …< 10; i++) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%" PetscInt_FMT " %g\n", i, (double)evals[i])); in main()
165 for (i = 0; i < m; i++) arrayU[j * m + i] *= evals[j]; in main()
173 PetscCall(PetscFree(evals)); in main()
|
| H A D | ex118.c | 21 PetscScalar *evecs_array, *D, *E, *evals; in main()
35 evals = E + n; in main()
49 …, &iu, &tol, (PetscReal *)D, (PetscReal *)E, &nevs, &nsplit, (PetscReal *)evals, iblock, isplit, w… in main()
55 …LAPACKstein_(&n, (PetscReal *)D, (PetscReal *)E, &nevs, (PetscReal *)evals, iblock, isplit, evecs_… in main()
61 …nevs; i++) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%" PetscInt_FMT " %g\n", i, (double)evals[i])); in main()
90 PetscCall(CkEigenSolutions(cklvl, T, il - 1, iu - 1, evals, evecs, tols)); in main()
|
| H A D | ex120.c | 16 PetscReal *evals, *rwork; in main() local
135 PetscCall(PetscMalloc1(n, &evals)); in main()
142 LAPACKsyev_("V", "U", &bn, arrayA, &bn, evals, work, &lwork, rwork, &lierr); in main()
163 …LAPACKsyevx_("V", "I", "U", &bn, arrayA, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_arr… in main()
172 LAPACKsygv_(&one, "V", "U", &bn, arrayA, &bn, arrayB, &bn, evals, work, &lwork, rwork, &lierr); in main()
192 …, "U", &bn, arrayA, &bn, arrayB, &bn, &vl, &vu, &il, &iu, &abstol, &nevs, evals, evecs_array, &nn,… in main()
204 …i++) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "%" PetscInt_FMT " %g\n", i + il, (double)evals[i])); in main()
218 PetscCall(CkEigenSolutions(cklvl, A, il - 1, iu - 1, evals, evecs, tols)); in main()
224 PetscCall(PetscFree(evals)); in main()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexmetric.c | 1469 PetscReal *evals; in DMPlexMetricIntersection_Private() local
1480 PetscCall(PetscMalloc4(dim * dim, &evecs, dim * dim, &sqrtM1, dim * dim, &isqrtM1, dim, &evals)); 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()
1519 sqrtj = PetscSqrtReal(evals[j]); 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()
1563 for (i = 0; i < dim; ++i) evals[i] = PetscMax(evals[i], 1.0); in DMPlexMetricIntersection_Private()
1571 …< dim; ++l) M2[i * dim + m] += sqrtM1[j * dim + i] * evecs[j * dim + k] * evals[k] * evecs[l * dim… in DMPlexMetricIntersection_Private()
1579 PetscCall(PetscFree4(evecs, sqrtM1, isqrtM1, evals)); in DMPlexMetricIntersection_Private()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 1417 const PetscScalar *evals; in PCBDDCNedelecSupport() local
1455 PetscCall(VecGetArrayRead(E, &evals)); in PCBDDCNedelecSupport()
1457 if (evals[i] == 0.0) eedgesidxs[cum++] = i + pc->pmat->rmap->rstart; in PCBDDCNedelecSupport()
1458 PetscCall(VecRestoreArrayRead(E, &evals)); in PCBDDCNedelecSupport()
|