|
Solve given a QR decomposition.
This function solves the system of equations using the results of a QR decomposition. This function requires the actual QR decomposition to be performed ahead of time; by qr_decomp() for example.
This function is intended for repeatedly solving systems of equations based on A. That is A stays constant while b varies.
- Note:
- This function requires BLAS/LAPACK 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 this 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:
-
A | A Matrix to decompose. |
b | A Matrix with as many rows as A. |
QR | A QRdecomp object containing the result of the QR decomposition of A. |
- See also:
- QRdecomp
-
qr_solve (const Matrix<>& A, const Matrix<>& b)
-
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:
-
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().
|