Chris@49
|
1 // Copyright (C) 2009-2010 NICTA (www.nicta.com.au)
|
Chris@49
|
2 // Copyright (C) 2009-2010 Conrad Sanderson
|
Chris@49
|
3 // Copyright (C) 2009-2010 Dimitrios Bouzas
|
Chris@49
|
4 //
|
Chris@49
|
5 // This Source Code Form is subject to the terms of the Mozilla Public
|
Chris@49
|
6 // License, v. 2.0. If a copy of the MPL was not distributed with this
|
Chris@49
|
7 // file, You can obtain one at http://mozilla.org/MPL/2.0/.
|
Chris@49
|
8
|
Chris@49
|
9
|
Chris@49
|
10
|
Chris@49
|
11 //! \addtogroup fn_eps
|
Chris@49
|
12 //! @{
|
Chris@49
|
13
|
Chris@49
|
14
|
Chris@49
|
15
|
Chris@49
|
16 //! \brief
|
Chris@49
|
17 //! eps version for non-complex matrices and vectors
|
Chris@49
|
18 template<typename T1>
|
Chris@49
|
19 inline
|
Chris@49
|
20 const eOp<T1, eop_eps>
|
Chris@49
|
21 eps(const Base<typename T1::elem_type, T1>& X, const typename arma_not_cx<typename T1::elem_type>::result* junk = 0)
|
Chris@49
|
22 {
|
Chris@49
|
23 arma_extra_debug_sigprint();
|
Chris@49
|
24 arma_ignore(junk);
|
Chris@49
|
25
|
Chris@49
|
26 return eOp<T1, eop_eps>(X.get_ref());
|
Chris@49
|
27 }
|
Chris@49
|
28
|
Chris@49
|
29
|
Chris@49
|
30
|
Chris@49
|
31 //! \brief
|
Chris@49
|
32 //! eps version for complex matrices and vectors
|
Chris@49
|
33 template<typename T1>
|
Chris@49
|
34 inline
|
Chris@49
|
35 Mat< typename T1::pod_type >
|
Chris@49
|
36 eps(const Base< std::complex<typename T1::pod_type>, T1>& X, const typename arma_cx_only<typename T1::elem_type>::result* junk = 0)
|
Chris@49
|
37 {
|
Chris@49
|
38 arma_extra_debug_sigprint();
|
Chris@49
|
39
|
Chris@49
|
40 arma_ignore(junk);
|
Chris@49
|
41
|
Chris@49
|
42 typedef typename T1::pod_type T;
|
Chris@49
|
43 typedef typename T1::elem_type eT;
|
Chris@49
|
44
|
Chris@49
|
45 const unwrap<T1> tmp(X.get_ref());
|
Chris@49
|
46 const Mat<eT>& A = tmp.M;
|
Chris@49
|
47
|
Chris@49
|
48 Mat<T> out(A.n_rows, A.n_cols);
|
Chris@49
|
49
|
Chris@49
|
50 T* out_mem = out.memptr();
|
Chris@49
|
51 const eT* A_mem = A.memptr();
|
Chris@49
|
52 const uword n_elem = A.n_elem;
|
Chris@49
|
53
|
Chris@49
|
54 for(uword i=0; i<n_elem; ++i)
|
Chris@49
|
55 {
|
Chris@49
|
56 out_mem[i] = eop_aux::direct_eps( A_mem[i] );
|
Chris@49
|
57 }
|
Chris@49
|
58
|
Chris@49
|
59
|
Chris@49
|
60 return out;
|
Chris@49
|
61 }
|
Chris@49
|
62
|
Chris@49
|
63
|
Chris@49
|
64
|
Chris@49
|
65 template<typename eT>
|
Chris@49
|
66 arma_inline
|
Chris@49
|
67 arma_warn_unused
|
Chris@49
|
68 typename arma_integral_only<eT>::result
|
Chris@49
|
69 eps(const eT& x)
|
Chris@49
|
70 {
|
Chris@49
|
71 arma_ignore(x);
|
Chris@49
|
72
|
Chris@49
|
73 return eT(0);
|
Chris@49
|
74 }
|
Chris@49
|
75
|
Chris@49
|
76
|
Chris@49
|
77
|
Chris@49
|
78 template<typename eT>
|
Chris@49
|
79 arma_inline
|
Chris@49
|
80 arma_warn_unused
|
Chris@49
|
81 typename arma_real_only<eT>::result
|
Chris@49
|
82 eps(const eT& x)
|
Chris@49
|
83 {
|
Chris@49
|
84 return eop_aux::direct_eps(x);
|
Chris@49
|
85 }
|
Chris@49
|
86
|
Chris@49
|
87
|
Chris@49
|
88
|
Chris@49
|
89 template<typename T>
|
Chris@49
|
90 arma_inline
|
Chris@49
|
91 arma_warn_unused
|
Chris@49
|
92 typename arma_real_only<T>::result
|
Chris@49
|
93 eps(const std::complex<T>& x)
|
Chris@49
|
94 {
|
Chris@49
|
95 return eop_aux::direct_eps(x);
|
Chris@49
|
96 }
|
Chris@49
|
97
|
Chris@49
|
98
|
Chris@49
|
99
|
Chris@49
|
100 //! @}
|