Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 2008, 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #ifndef MLN_NORM_L1_HH 00027 # define MLN_NORM_L1_HH 00028 00033 00034 # include <mln/math/abs.hh> 00035 # include <mln/algebra/vec.hh> 00036 00037 00038 namespace mln 00039 { 00040 00041 namespace norm 00042 { 00043 00046 template <unsigned n, typename C> 00047 mln_sum_product(C,C) l1(const C (&vec)[n]); 00048 00049 template <unsigned n, typename C> 00050 mln_sum_product(C,C) l1(const algebra::vec<n,C>& vec); 00052 00055 template <unsigned n, typename C> 00056 mln_sum_product(C,C) l1_distance(const C (&vec1)[n], const C (&vec2)[n]); 00057 00058 template <unsigned n, typename C> 00059 mln_sum_product(C,C) l1_distance(const algebra::vec<n,C>& vec1, 00060 const algebra::vec<n,C>& vec2); 00062 00063 00064 # ifndef MLN_INCLUDE_ONLY 00065 00066 namespace impl 00067 { 00068 template <unsigned n, typename C, typename V> 00069 inline 00070 mln_sum_product(C,C) 00071 l1_(const V& vec) 00072 { 00073 typedef mln_sum_product(C,C) M; 00074 M m = 0; 00075 for (unsigned i = 0; i < n; ++i) 00076 { 00077 M v_i = static_cast<M>(mln::math::abs(vec[i])); 00078 m = static_cast<M>(m + v_i); 00079 } 00080 return m; 00081 } 00082 00083 template <unsigned n, typename C, typename V> 00084 inline 00085 mln_sum_product(C,C) 00086 l1_distance_(const V& vec1, const V& vec2) 00087 { 00088 typedef mln_sum_product(C,C) D; 00089 D d = 0; 00090 for (unsigned i = 0; i < n; ++i) 00091 { 00092 D v1_v2 = static_cast<D>(mln::math::abs(vec1[i] - vec2[i])); 00093 d = static_cast<D>(d + v1_v2); 00094 } 00095 return d; 00096 } 00097 00098 } // end of namespace mln::norm::impl 00099 00100 00101 /*----------. 00102 | Facades. | 00103 `----------*/ 00104 00105 template <unsigned n, typename C> 00106 inline 00107 mln_sum_product(C,C) 00108 l1(const C (&vec)[n]) 00109 { 00110 return impl::l1_<n, C>(vec); 00111 } 00112 00113 template <unsigned n, typename C> 00114 inline 00115 mln_sum_product(C,C) 00116 l1(const algebra::vec<n,C>& vec) 00117 { 00118 return impl::l1_<n, C>(vec); 00119 } 00120 00121 template <unsigned n, typename C> 00122 inline 00123 mln_sum_product(C,C) 00124 l1_distance(const C (&vec1)[n], const C (&vec2)[n]) 00125 { 00126 return impl::l1_distance_<n, C>(vec1, vec2); 00127 } 00128 00129 template <unsigned n, typename C> 00130 inline 00131 mln_sum_product(C,C) 00132 l1_distance(const algebra::vec<n,C>& vec1, const algebra::vec<n,C>& vec2) 00133 { 00134 return impl::l1_distance_<n, C>(vec1, vec2); 00135 } 00136 00137 # endif // ! MLN_INCLUDE_ONLY 00138 00139 } // end of namespace mln::norm 00140 00141 } // end of namespace mln 00142 00143 00144 #endif // ! MLN_NORM_L1_HH