annotate armadillo-3.900.4/include/armadillo_bits/fn_eps.hpp @ 84:55a047986812 tip

Update library URI so as not to be document-local
author Chris Cannam
date Wed, 22 Apr 2020 14:21:57 +0100
parents 1ec0e2823891
children
rev   line source
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 //! @}