00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #ifndef MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH
00027 # define MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH
00028
00033
00034 # include <mln/core/concept/function.hh>
00035 # include <mln/algebra/h_mat.hh>
00036 # include <mln/algebra/h_vec.hh>
00037
00038 namespace mln
00039 {
00040
00041 namespace fun
00042 {
00043
00044 namespace internal
00045 {
00046
00047 template <typename V, typename E>
00048 struct x2x_linear_impl_
00049 {
00050 static const unsigned dim = V::dim;
00051
00052 typedef V argument;
00053 typedef V result;
00054 typedef typename V::coord coord;
00055 typedef algebra::h_mat<dim, coord> matrix;
00056
00057 V operator()(const V& x) const
00058 {
00059 algebra::h_vec<dim, coord> 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 E>
00076 inline
00077 x2x_linear_impl_<V,E>::x2x_linear_impl_()
00078 {
00079 }
00080
00081
00082 template <typename V, typename E>
00083 inline
00084 const typename x2x_linear_impl_<V,E>::matrix&
00085 x2x_linear_impl_<V,E>::mat() const
00086 {
00087 return m_;
00088 }
00089
00090
00091 # endif // ! MLN_INCLUDE_ONLY
00092
00093 }
00094
00095 }
00096
00097 }
00098
00099
00100 #endif // ! MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH