Uses the LU factorization to compute the solution to the system of linear equations with a square matrix A and multiple right-hand sides, and provides error bounds on the solution.
Computes selected eigenvalues and, optionally, eigenvectors of a Hermitian matrix.
Computes an LU factorization of a general band matrix, using partial pivoting with row interchanges. The routine is called by p?dttrs.
Reduces the first nb columns of a general rectangular matrix A so that elements below the k-th subdiagonal are zero, by an orthogonal/unitary transformation, and returns auxiliary matrices that are needed to apply the transformation to the unreduced part of A.
Scales a symmetric/Hermitian matrix, using scaling factors computed by p?poequ .
Sorts the numbers in an array and the corresponding vectors in increasing order.
Multiplies a general matrix by the orthogonal/unitary matrix from an LQ factorization determined by p?gelqf (unblocked algorithm).
Computes scaled eigenvector corresponding to given eigenvalue.
Computes an LU factorization of a general tridiagonal matrix with no pivoting (local blocked algorithm).
Copies a submatrix from one trapezoidal matrix to another.