fn_qr.hpp

Go to the documentation of this file.
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 //! @}