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
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;
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 }
00094
00095 }
00096
00097 }
00098
00099
00100 #endif // ! MLN_FUN_INTERNAL_X2X_LINEAR_IMPL_HH