Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 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_FUN_STAT_MAHALANOBIS_HH 00027 # define MLN_FUN_STAT_MAHALANOBIS_HH 00028 00032 00033 # include <cmath> 00034 # include <mln/core/concept/function.hh> 00035 # include <mln/algebra/vec.hh> 00036 # include <mln/algebra/mat.hh> 00037 00038 00039 namespace mln 00040 { 00041 00042 namespace fun 00043 { 00044 00045 namespace stat 00046 { 00047 00048 template <typename V> 00049 struct mahalanobis 00050 : public Function_v2v< mahalanobis<V> >, 00051 private metal::equal< V, algebra::vec<V::dim,float> >::check_t 00052 { 00053 enum { n = V::dim }; 00054 typedef float result; 00055 00056 mahalanobis(const algebra::mat<n,n,float>& var, 00057 const algebra::vec<n,float>& mean); 00058 // Tech. note: using n (instead of V::dim) above helps g++-2.95. 00059 00060 float operator()(const V& v) const; 00061 00062 typedef algebra::vec<n,float> mean_t; 00063 00064 mean_t mean() const; 00065 00066 protected: 00067 algebra::mat<n,n,float> var_1_; 00068 algebra::vec<n,float> mean_; 00069 }; 00070 00071 00072 # ifndef MLN_INCLUDE_ONLY 00073 00074 template <typename V> 00075 inline 00076 mahalanobis<V>::mahalanobis(const algebra::mat<n,n,float>& var, 00077 const algebra::vec<n,float>& mean) 00078 { 00079 var_1_ = var._1(); 00080 mean_ = mean; 00081 } 00082 00083 template <typename V> 00084 inline 00085 float 00086 mahalanobis<V>::operator()(const V& v) const 00087 { 00088 return std::sqrt((v - mean_).t() * var_1_ * (v - mean_)); 00089 } 00090 00091 template <typename V> 00092 inline 00093 typename mahalanobis<V>::mean_t 00094 mahalanobis<V>::mean() const 00095 { 00096 return mean_; 00097 } 00098 00099 # endif // ! MLN_INCLUDE_ONLY 00100 00101 } // end of namespace mln::fun::stat 00102 00103 } // end of namespace mln::fun 00104 00105 } // end of namespace mln 00106 00107 00108 #endif // ! MLN_FUN_STAT_MAHALANOBIS_HH