00001 // Copyright (C) 2009 NICTA 00002 // 00003 // Authors: 00004 // - Conrad Sanderson (conradsand at ieee dot org) 00005 // 00006 // This file is part of the Armadillo C++ library. 00007 // It is provided without any warranty of fitness 00008 // for any purpose. You can redistribute this file 00009 // and/or modify it under the terms of the GNU 00010 // Lesser General Public License (LGPL) as published 00011 // by the Free Software Foundation, either version 3 00012 // of the License or (at your option) any later version. 00013 // (see http://www.opensource.org/licenses for more info) 00014 00015 00016 //! \addtogroup fn_qr 00017 //! @{ 00018 00019 00020 //! QR decomposition 00021 template<typename T1> 00022 inline 00023 void 00024 qr(Mat<typename T1::elem_type>& Q, Mat<typename T1::elem_type>& R, const Base<typename T1::elem_type,T1>& X) 00025 { 00026 arma_extra_debug_sigprint(); 00027 00028 typedef typename T1::elem_type eT; 00029 00030 arma_debug_check( (&Q == &R), "qr(): Q and R are the same object"); 00031 00032 const unwrap_check<T1> tmp1(X.get_ref(), Q); 00033 const Mat<eT>& A = tmp1.M; 00034 00035 const unwrap_check< Mat<eT> > tmp2(A, R); 00036 const Mat<eT>& B = tmp2.M; 00037 00038 const bool ok = auxlib::qr(Q, R, B); 00039 00040 if(ok == false) 00041 { 00042 arma_print("qr(): matrix factorisation failed"); 00043 } 00044 00045 } 00046 00047 00048 //! @}