Searched refs:subdiag (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1418 …DTSymmetricTridiagonalEigensolve(PetscInt n, PetscReal diag[], PetscReal subdiag[], PetscReal eigs… in PetscDTSymmetricTridiagonalEigensolve() argument
1446 …PetscCallBLAS("LAPACKstegr", LAPACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &… in PetscDTSymmetricTridiagonalEigensolve()
1452 …PetscCallBLAS("LAPACKstegr", LAPACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &… in PetscDTSymmetricTridiagonalEigensolve()
1462 PetscCallBLAS("LAPACKsteqr", LAPACKsteqr_("I", &bn, diag, subdiag, V, &ldz, work, &info)); in PetscDTSymmetricTridiagonalEigensolve()
1657 PetscReal *diag, *subdiag; in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal() local
1660 PetscCall(PetscMalloc2(npoints, &diag, npoints, &subdiag)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1662 PetscCall(PetscDTJacobiMatrix_Internal(npoints, a, b, diag, subdiag)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1663 for (i = 0; i < npoints - 1; i++) subdiag[i] = PetscSqrtReal(subdiag[i]); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1664 PetscCall(PetscDTSymmetricTridiagonalEigensolve(npoints, diag, subdiag, x, V)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
1667 PetscCall(PetscFree2(diag, subdiag)); in PetscDTGaussJacobiQuadrature_GolubWelsch_Internal()
|