Scythe-1.0.3
Matrix scythe::qr_solve ( const Matrix<> &  A,
const Matrix<> &  b 
) [inline]

Solve $Ax=b$ using QR decomposition.

This function solves the system of equations $Ax = b$ using QR decomposition. This function is intended for repeatedly solving systems of equations based on A. That is A stays constant while b varies.

Note:
This function used BLAS/LAPACK support functionality and is only available on machines that provide these libraries. Make sure you enable the SCYTHE_LAPACK preprocessor flag if you wish to use this function. Furthermore, note that the function takes and returns only column-major concrete matrices. Future versions of Scythe will provide a native C++ implementation of this function with support for general matrix templates.
Parameters:
AA Matrix to decompose.
bA Matrix with as many rows as A.
See also:
QRdecomp
qr_solve (const Matrix<>& A, const Matrix<>& b, const QRdecomp& QR)
qr_decomp (const Matrix<>& A)
lu_solve (const Matrix<T,PO1,PS1>&, const Matrix<T,PO2,PS2>&, const Matrix<T,PO3,PS3>&, const Matrix<T,PO4,PS4>&, const Matrix<unsigned int, PO5, PS5>&)
lu_solve (Matrix<T,PO1,PS1>, const Matrix<T,PO2,PS2>&)
chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &)
chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &, const Matrix<T,PO3,PS3> &)
Exceptions:
scythe_null_error(Level 1)
scythe_conformation_error(Level 1)
scythe_type_error(Level 1)
scythe_lapack_internal_error(Level 1)

References scythe::Matrix_base< ORDER, STYLE >::cols(), scythe::Matrix< T_type, ORDER, STYLE >::getArray(), scythe::Matrix_base< ORDER, STYLE >::isNull(), scythe::Matrix_base< ORDER, STYLE >::rows(), SCYTHE_CHECK_10, and scythe::Matrix_base< ORDER, STYLE >::size().