I need to compute the `k`

smallest eigenvalues of a symmetric positive definite real matrix, and their corresponding eigenvectors. Can I use a Julia package for this, or is there something under LinearAlgebra that helps?

`LinearAlgebra.eigen`

takes `vl,vu`

parameters, to return eigenvalues within `[vl,vu]`

, but I do not know *a priori* the magnitude of the eigenvalues, so I wouldn’t know what `vl,vu`

to use to get the `k`

smallest ones.