Mercurial > hg > segmenter-vamp-plugin
annotate armadillo-2.4.4/include/armadillo_bits/fn_trunc_log.hpp @ 0:8b6102e2a9b0
Armadillo Library
| author | maxzanoni76 <max.zanoni@eecs.qmul.ac.uk> |
|---|---|
| date | Wed, 11 Apr 2012 09:27:06 +0100 |
| parents | |
| children |
| rev | line source |
|---|---|
| max@0 | 1 // Copyright (C) 2008-2010 NICTA (www.nicta.com.au) |
| max@0 | 2 // Copyright (C) 2008-2010 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_trunc_log |
| max@0 | 15 //! @{ |
| max@0 | 16 |
| max@0 | 17 |
| max@0 | 18 |
| max@0 | 19 template<typename eT> |
| max@0 | 20 inline |
| max@0 | 21 static |
| max@0 | 22 typename arma_float_only<eT>::result |
| max@0 | 23 trunc_log(const eT x) |
| max@0 | 24 { |
| max@0 | 25 if(std::numeric_limits<eT>::is_iec559) |
| max@0 | 26 { |
| max@0 | 27 if(x == std::numeric_limits<eT>::infinity()) |
| max@0 | 28 { |
| max@0 | 29 return Math<eT>::log_max(); |
| max@0 | 30 } |
| max@0 | 31 else |
| max@0 | 32 { |
| max@0 | 33 return (x <= eT(0)) ? Math<eT>::log_min() : std::log(x); |
| max@0 | 34 } |
| max@0 | 35 } |
| max@0 | 36 else |
| max@0 | 37 { |
| max@0 | 38 return std::log(x); |
| max@0 | 39 } |
| max@0 | 40 } |
| max@0 | 41 |
| max@0 | 42 |
| max@0 | 43 |
| max@0 | 44 template<typename eT> |
| max@0 | 45 inline |
| max@0 | 46 static |
| max@0 | 47 typename arma_integral_only<eT>::result |
| max@0 | 48 trunc_log(const eT x) |
| max@0 | 49 { |
| max@0 | 50 return eT( trunc_log( double(x) ) ); |
| max@0 | 51 } |
| max@0 | 52 |
| max@0 | 53 |
| max@0 | 54 |
| max@0 | 55 template<typename T> |
| max@0 | 56 inline |
| max@0 | 57 static |
| max@0 | 58 std::complex<T> |
| max@0 | 59 trunc_log(const std::complex<T>& x) |
| max@0 | 60 { |
| max@0 | 61 return std::log(x); |
| max@0 | 62 } |
| max@0 | 63 |
| max@0 | 64 |
| max@0 | 65 |
| max@0 | 66 template<typename T1> |
| max@0 | 67 arma_inline |
| max@0 | 68 const eOp<T1, eop_trunc_log> |
| max@0 | 69 trunc_log(const Base<typename T1::elem_type,T1>& A) |
| max@0 | 70 { |
| max@0 | 71 arma_extra_debug_sigprint(); |
| max@0 | 72 |
| max@0 | 73 return eOp<T1, eop_trunc_log>(A.get_ref()); |
| max@0 | 74 } |
| max@0 | 75 |
| max@0 | 76 |
| max@0 | 77 |
| max@0 | 78 template<typename T1> |
| max@0 | 79 arma_inline |
| max@0 | 80 const eOpCube<T1, eop_trunc_log> |
| max@0 | 81 trunc_log(const BaseCube<typename T1::elem_type,T1>& A) |
| max@0 | 82 { |
| max@0 | 83 arma_extra_debug_sigprint(); |
| max@0 | 84 |
| max@0 | 85 return eOpCube<T1, eop_trunc_log>(A.get_ref()); |
| max@0 | 86 } |
| max@0 | 87 |
| max@0 | 88 |
| max@0 | 89 |
| max@0 | 90 //! @} |