annotate armadillo-2.4.4/include/armadillo_bits/fn_qr.hpp @ 18:8d046a9d36aa slimline

Back out rev 13:ac07c60aa798. Like an idiot, I committed a whole pile of unrelated changes in the guise of a single typo fix. Will re-commit in stages
author Chris Cannam
date Thu, 10 May 2012 10:45:44 +0100
parents 8b6102e2a9b0
children
rev   line source
max@0 1 // Copyright (C) 2009-2012 NICTA (www.nicta.com.au)
max@0 2 // Copyright (C) 2009-2012 Conrad Sanderson
max@0 3 //
max@0 4 // This file is part of the Armadillo C++ library.
max@0 5 // It is provided without any warranty of fitness
max@0 6 // for any purpose. You can redistribute this file
max@0 7 // and/or modify it under the terms of the GNU
max@0 8 // Lesser General Public License (LGPL) as published
max@0 9 // by the Free Software Foundation, either version 3
max@0 10 // of the License or (at your option) any later version.
max@0 11 // (see http://www.opensource.org/licenses for more info)
max@0 12
max@0 13
max@0 14 //! \addtogroup fn_qr
max@0 15 //! @{
max@0 16
max@0 17
max@0 18
max@0 19 //! QR decomposition
max@0 20 template<typename T1>
max@0 21 inline
max@0 22 bool
max@0 23 qr
max@0 24 (
max@0 25 Mat<typename T1::elem_type>& Q,
max@0 26 Mat<typename T1::elem_type>& R,
max@0 27 const Base<typename T1::elem_type,T1>& X,
max@0 28 const typename arma_blas_type_only<typename T1::elem_type>::result* junk = 0
max@0 29 )
max@0 30 {
max@0 31 arma_extra_debug_sigprint();
max@0 32 arma_ignore(junk);
max@0 33
max@0 34 arma_debug_check( (&Q == &R), "qr(): Q and R are the same object");
max@0 35
max@0 36 const bool status = auxlib::qr(Q, R, X);
max@0 37
max@0 38 if(status == false)
max@0 39 {
max@0 40 Q.reset();
max@0 41 R.reset();
max@0 42 arma_bad("qr(): failed to converge", false);
max@0 43 }
max@0 44
max@0 45 return status;
max@0 46 }
max@0 47
max@0 48
max@0 49
max@0 50 //! @}