|
max@0
|
1 // Copyright (C) 2008-2011 NICTA (www.nicta.com.au)
|
|
max@0
|
2 // Copyright (C) 2008-2011 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_inv
|
|
max@0
|
15 //! @{
|
|
max@0
|
16
|
|
max@0
|
17
|
|
max@0
|
18
|
|
max@0
|
19 //! delayed matrix inverse (general matrices)
|
|
max@0
|
20 template<typename T1>
|
|
max@0
|
21 arma_inline
|
|
max@0
|
22 const Op<T1, op_inv>
|
|
max@0
|
23 inv
|
|
max@0
|
24 (
|
|
max@0
|
25 const Base<typename T1::elem_type,T1>& X,
|
|
max@0
|
26 const bool slow = false,
|
|
max@0
|
27 const typename arma_blas_type_only<typename T1::elem_type>::result* junk = 0
|
|
max@0
|
28 )
|
|
max@0
|
29 {
|
|
max@0
|
30 arma_extra_debug_sigprint();
|
|
max@0
|
31 arma_ignore(junk);
|
|
max@0
|
32
|
|
max@0
|
33 return Op<T1, op_inv>(X.get_ref(), ((slow == false) ? 0 : 1), 0);
|
|
max@0
|
34 }
|
|
max@0
|
35
|
|
max@0
|
36
|
|
max@0
|
37
|
|
max@0
|
38 //! remove the inverse operation if applied twice consecutively
|
|
max@0
|
39 template<typename T1>
|
|
max@0
|
40 arma_inline
|
|
max@0
|
41 const T1&
|
|
max@0
|
42 inv(const Op<T1, op_inv>& X, const bool slow = false)
|
|
max@0
|
43 {
|
|
max@0
|
44 arma_extra_debug_sigprint();
|
|
max@0
|
45 arma_ignore(slow);
|
|
max@0
|
46
|
|
max@0
|
47 return X.m;
|
|
max@0
|
48 }
|
|
max@0
|
49
|
|
max@0
|
50
|
|
max@0
|
51
|
|
max@0
|
52 //! delayed matrix inverse (triangular matrices)
|
|
max@0
|
53 template<typename T1>
|
|
max@0
|
54 arma_inline
|
|
max@0
|
55 const Op<T1, op_inv_tr>
|
|
max@0
|
56 inv
|
|
max@0
|
57 (
|
|
max@0
|
58 const Op<T1, op_trimat>& X,
|
|
max@0
|
59 const bool slow = false,
|
|
max@0
|
60 const typename arma_blas_type_only<typename T1::elem_type>::result* junk = 0
|
|
max@0
|
61 )
|
|
max@0
|
62 {
|
|
max@0
|
63 arma_extra_debug_sigprint();
|
|
max@0
|
64 arma_ignore(slow);
|
|
max@0
|
65 arma_ignore(junk);
|
|
max@0
|
66
|
|
max@0
|
67 return Op<T1, op_inv_tr>(X.m, X.aux_uword_a, 0);
|
|
max@0
|
68 }
|
|
max@0
|
69
|
|
max@0
|
70
|
|
max@0
|
71
|
|
max@0
|
72 //! delayed matrix inverse (symmetric positive definite matrices)
|
|
max@0
|
73 template<typename T1>
|
|
max@0
|
74 arma_inline
|
|
max@0
|
75 const Op<T1, op_inv_sympd>
|
|
max@0
|
76 inv
|
|
max@0
|
77 (
|
|
max@0
|
78 const Op<T1, op_sympd>& X,
|
|
max@0
|
79 const bool slow = false,
|
|
max@0
|
80 const typename arma_blas_type_only<typename T1::elem_type>::result* junk = 0
|
|
max@0
|
81 )
|
|
max@0
|
82 {
|
|
max@0
|
83 arma_extra_debug_sigprint();
|
|
max@0
|
84 arma_ignore(slow);
|
|
max@0
|
85 arma_ignore(junk);
|
|
max@0
|
86
|
|
max@0
|
87 return Op<T1, op_inv_sympd>(X.m, 0, 0);
|
|
max@0
|
88 }
|
|
max@0
|
89
|
|
max@0
|
90
|
|
max@0
|
91
|
|
max@0
|
92 template<typename T1>
|
|
max@0
|
93 inline
|
|
max@0
|
94 bool
|
|
max@0
|
95 inv
|
|
max@0
|
96 (
|
|
max@0
|
97 Mat<typename T1::elem_type>& out,
|
|
max@0
|
98 const Base<typename T1::elem_type,T1>& X,
|
|
max@0
|
99 const bool slow = false,
|
|
max@0
|
100 const typename arma_blas_type_only<typename T1::elem_type>::result* junk = 0
|
|
max@0
|
101 )
|
|
max@0
|
102 {
|
|
max@0
|
103 arma_extra_debug_sigprint();
|
|
max@0
|
104 arma_ignore(junk);
|
|
max@0
|
105
|
|
max@0
|
106 try
|
|
max@0
|
107 {
|
|
max@0
|
108 out = inv(X,slow);
|
|
max@0
|
109 }
|
|
max@0
|
110 catch(std::runtime_error&)
|
|
max@0
|
111 {
|
|
max@0
|
112 return false;
|
|
max@0
|
113 }
|
|
max@0
|
114
|
|
max@0
|
115 return true;
|
|
max@0
|
116 }
|
|
max@0
|
117
|
|
max@0
|
118
|
|
max@0
|
119
|
|
max@0
|
120 //! @}
|
