Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 2008, 2009 EPITA Research and Development 00002 // Laboratory (LRDE) 00003 // 00004 // This file is part of Olena. 00005 // 00006 // Olena is free software: you can redistribute it and/or modify it under 00007 // the terms of the GNU General Public License as published by the Free 00008 // Software Foundation, version 2 of the License. 00009 // 00010 // Olena is distributed in the hope that it will be useful, 00011 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00012 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00013 // General Public License for more details. 00014 // 00015 // You should have received a copy of the GNU General Public License 00016 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00017 // 00018 // As a special exception, you may use this file as part of a free 00019 // software project without restriction. Specifically, if other files 00020 // instantiate templates or use macros or inline functions from this 00021 // file, or you compile this file and link it with other files to produce 00022 // an executable, this file does not by itself cause the resulting 00023 // executable to be covered by the GNU General Public License. This 00024 // exception does not however invalidate any other reasons why the 00025 // executable file might be covered by the GNU General Public License. 00026 00027 #ifndef MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH 00028 # define MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH 00029 00034 00035 # include <mln/core/concept/function.hh> 00036 # include <mln/algebra/h_mat.hh> 00037 # include <mln/algebra/h_vec.hh> 00038 00039 namespace mln 00040 { 00041 00042 namespace fun 00043 { 00044 00045 namespace internal 00046 { 00047 00048 template <typename V, typename C, typename E> 00049 struct x2x_linear_impl_ 00050 { 00051 static const unsigned dim = V::dim; // To please g++ when comparing dim in metal bexpr. 00052 00053 typedef V argument; 00054 typedef V result; 00055 typedef algebra::h_mat<dim, C> matrix; 00056 00057 V operator()(const V& x) const 00058 { 00059 algebra::h_vec<dim, C> tmp = m_ * x.to_h_vec(); 00060 return tmp.to_vec(); 00061 } 00062 00063 const matrix& mat() const; 00064 00065 protected: 00066 x2x_linear_impl_(); 00067 00068 matrix m_; 00069 }; 00070 00071 00072 00073 # ifndef MLN_INCLUDE_ONLY 00074 00075 template <typename V, typename C, typename E> 00076 inline 00077 x2x_linear_impl_<V,C,E>::x2x_linear_impl_() 00078 { 00079 } 00080 00081 00082 template <typename V, typename C, typename E> 00083 inline 00084 const typename x2x_linear_impl_<V,C,E>::matrix& 00085 x2x_linear_impl_<V,C,E>::mat() const 00086 { 00087 return m_; 00088 } 00089 00090 00091 # endif // ! MLN_INCLUDE_ONLY 00092 00093 } // end of namespace mln::fun::internal 00094 00095 } // end of namespace mln::fun 00096 00097 } // end of namespace mln 00098 00099 00100 #endif // ! MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH