1 /*
2 [auto_generated]
3 boost/numeric/odeint/external/eigen/eigen_algebra.hpp
4
5 [begin_description]
6 tba.
7 [end_description]
8
9 Copyright 2013 Christian Shelton
10 Copyright 2013 Karsten Ahnert
11
12 Distributed under the Boost Software License, Version 1.0.
13 (See accompanying file LICENSE_1_0.txt or
14 copy at http://www.boost.org/LICENSE_1_0.txt)
15 */
16
17
18 #ifndef BOOST_NUMERIC_ODEINT_EXTERNAL_EIGEN_EIGEN_ALGEBRA_HPP_INCLUDED
19 #define BOOST_NUMERIC_ODEINT_EXTERNAL_EIGEN_EIGEN_ALGEBRA_HPP_INCLUDED
20
21 #include <Eigen/Dense>
22 #include <boost/numeric/odeint/algebra/vector_space_algebra.hpp>
23
24 // Necessary routines for Eigen matrices to work with vector_space_algebra
25 // from odeint
26 // (that is, it lets odeint treat the eigen matrices correctly, knowing
27 // how to add, multiply, compute the norm, etc)
28 namespace Eigen {
29
30 template<typename D>
31 inline const
32 typename Eigen::CwiseBinaryOp<
33 internal::scalar_sum_op<typename internal::traits<D>::Scalar>,
34 typename DenseBase<D>::ConstantReturnType,
35 const D>
operator +(const typename Eigen::MatrixBase<D> & m,const typename Eigen::internal::traits<D>::Scalar & s)36 operator+(const typename Eigen::MatrixBase<D> &m,
37 const typename Eigen::internal::traits<D>::Scalar &s) {
38 return CwiseBinaryOp<
39 internal::scalar_sum_op<typename internal::traits<D>::Scalar>,
40 typename DenseBase<D>::ConstantReturnType,
41 const D>(DenseBase<D>::Constant(m.rows(), m.cols(), s), m.derived());
42 }
43
44 template<typename D>
45 inline const
46 typename Eigen::CwiseBinaryOp<
47 internal::scalar_sum_op<typename internal::traits<D>::Scalar>,
48 typename DenseBase<D>::ConstantReturnType,
49 const D>
operator +(const typename Eigen::internal::traits<D>::Scalar & s,const typename Eigen::MatrixBase<D> & m)50 operator+(const typename Eigen::internal::traits<D>::Scalar &s,
51 const typename Eigen::MatrixBase<D> &m) {
52 return CwiseBinaryOp<
53 internal::scalar_sum_op<typename internal::traits<D>::Scalar>,
54 typename DenseBase<D>::ConstantReturnType,
55 const D>(DenseBase<D>::Constant(m.rows(), m.cols(), s), m.derived());
56 }
57
58 template<typename D1,typename D2>
59 inline const
60 typename Eigen::CwiseBinaryOp<
61 typename Eigen::internal::scalar_quotient_op<
62 typename Eigen::internal::traits<D1>::Scalar>,
63 const D1, const D2>
operator /(const Eigen::MatrixBase<D1> & x1,const Eigen::MatrixBase<D2> & x2)64 operator/(const Eigen::MatrixBase<D1> &x1, const Eigen::MatrixBase<D2> &x2) {
65 return x1.cwiseQuotient(x2);
66 }
67
68
69 template< typename D >
70 inline const
71 typename Eigen::CwiseUnaryOp<
72 typename Eigen::internal::scalar_abs_op<
73 typename Eigen::internal::traits< D >::Scalar > ,
74 const D >
abs(const Eigen::MatrixBase<D> & m)75 abs( const Eigen::MatrixBase< D > &m ) {
76 return m.cwiseAbs();
77 }
78
79 } // end Eigen namespace
80
81
82 namespace boost {
83 namespace numeric {
84 namespace odeint {
85
86 template<typename B,int S1,int S2,int O, int M1, int M2>
87 struct vector_space_norm_inf< Eigen::Matrix<B,S1,S2,O,M1,M2> >
88 {
89 typedef B result_type;
operator ()boost::numeric::odeint::vector_space_norm_inf90 result_type operator()( const Eigen::Matrix<B,S1,S2,O,M1,M2> &m ) const
91 {
92 return m.template lpNorm<Eigen::Infinity>();
93 }
94 };
95
96 } } } // end boost::numeric::odeint namespace
97
98 #endif // BOOST_NUMERIC_ODEINT_EXTERNAL_EIGEN_EIGEN_ALGEBRA_HPP_INCLUDED
99