Robuta

https://www.intel.com/content/www/us/en/docs/onemkl/developer-reference-fortran/2023-1/lasq2.html ?lasq2 Computes all the eigenvalues of the symmetric positive definite tridiagonal matrix associated with the quotient difference array z to high relative accuracy....