Template Function mcf_qr_decompose

• Defined in File matrix_continued_fraction.hpp

Function Documentation

template<typename T>
void mcf_qr_decompose(const matrix<T> &C, matrix<T> &Q, matrix<T> &R)

QR decomposition of a p×p matrix via modified Gram-Schmidt: C = Q*R, Q unitary, R upper-triangular.

Used by matrix_continued_fraction::convert_B_format() to move an off-diagonal Lanczos block from Hermitian to upper-triangular form.