[Boost-commit] svn:boost r72764 - in sandbox/odeint/branches/karsten: boost/numeric/odeint/util libs/numeric/odeint/ideas/state_wrapper

Subject: [Boost-commit] svn:boost r72764 - in sandbox/odeint/branches/karsten: boost/numeric/odeint/util libs/numeric/odeint/ideas/state_wrapper
From: mario.mulansky_at_[hidden]
Date: 2011-06-26 10:03:30

Author: mariomulansky
Date: 2011-06-26 10:03:29 EDT (Sun, 26 Jun 2011)
New Revision: 72764
URL: http://svn.boost.org/trac/boost/changeset/72764

Log:
added resizing to wrapper in ideas
Added:
   sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_array.cpp   (contents, props changed)
Binary files modified: 
   sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/Jamfile
Text files modified: 
   sandbox/odeint/branches/karsten/boost/numeric/odeint/util/is_resizeable.hpp                 |     6 ++++--                                  
   sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/explicit_euler.hpp  |    39 +++++++++++++++++++++++++++++++++++---- 
   sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_gsl_vector.cpp |    13 +++++++++++++                           
   3 files changed, 52 insertions(+), 6 deletions(-)

Modified: sandbox/odeint/branches/karsten/boost/numeric/odeint/util/is_resizeable.hpp
==============================================================================
--- sandbox/odeint/branches/karsten/boost/numeric/odeint/util/is_resizeable.hpp	(original)
+++ sandbox/odeint/branches/karsten/boost/numeric/odeint/util/is_resizeable.hpp	2011-06-26 10:03:29 EDT (Sun, 26 Jun 2011)
@@ -23,7 +23,8 @@
 template< class Container >
 struct is_resizeable
 {
-	struct type : public boost::false_type { };
+	//struct type : public boost::false_type { };
+	typedef boost::false_type type;
         const static bool value = type::value;
 };
 
@@ -33,7 +34,8 @@
 template< class V, class A >
 struct is_resizeable< std::vector< V , A  > >
 {
-	struct type : public boost::true_type { };
+	//struct type : public boost::true_type { };
+    typedef boost::true_type type;
         const static bool value = type::value;
 };
 

Modified: sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/Jamfile
==============================================================================
Binary files. No diff available.

Modified: sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/explicit_euler.hpp
==============================================================================
--- sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/explicit_euler.hpp	(original)
+++ sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/explicit_euler.hpp	2011-06-26 10:03:29 EDT (Sun, 26 Jun 2011)
@@ -4,10 +4,40 @@
 
 #include <boost/numeric/odeint/algebra/range_algebra.hpp>
 #include <boost/numeric/odeint/algebra/default_operations.hpp>
-#include <boost/numeric/odeint/util/resize.hpp>
+#include <boost/numeric/odeint/util/is_resizeable.hpp>
+#include <boost/type_traits/integral_constant.hpp>
+
+
+template< class V , typename resizeable = typename boost::numeric::odeint::is_resizeable< V >::type>
+struct state_wrapper;
+
+//two standard implementations, with and without resizing depending on is_resizeable< StateType >
+
+template< class V >
+struct state_wrapper< V , boost::true_type > // with resizing
+{
+    typedef typename V::value_type value_type;
+
+    V m_v;
+
+    state_wrapper() : m_v() { }
+
+    template< class StateIn >
+    bool same_size( const StateIn &x ) const
+    {
+        return ( boost::size( m_v ) == boost::size( x ) );
+    }
+
+    template< class StateIn >
+    void resize( const StateIn &x )
+    {
+        m_v.resize( x.size() );
+    }
+};
+
 
 template< class V >
-struct state_wrapper
+struct state_wrapper< V , boost::false_type > // without resizing
 {
     typedef typename V::value_type value_type;
 
@@ -15,6 +45,7 @@
 
     state_wrapper() : m_v() { }
 
+    //no resize method
 };
 
 
@@ -45,9 +76,9 @@
     template< class State >
     bool adjust_size( const State &x )
     {
-        if( boost::numeric::odeint::same_size( x , m_dxdt.m_v ) )
+        if( m_dxdt.same_size( x ) )
         {
-            boost::numeric::odeint::resize( x , m_dxdt.m_v );
+            m_dxdt.resize( x );
             return true;
         } else
             return false;

Added: sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_array.cpp
==============================================================================
--- (empty file)
+++ sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_array.cpp	2011-06-26 10:03:29 EDT (Sun, 26 Jun 2011)
@@ -0,0 +1,34 @@
+#include <iostream>
+#include <boost/array.hpp>
+
+#include "explicit_euler.hpp"
+#include "size_adjuster.hpp"
+
+using namespace std;
+
+typedef boost::array< double , 3 > state_type;
+
+const double sigma = 10.0;
+const double R = 28.0;
+const double b = 8.0 / 3.0;
+
+void lorenz( const state_type &x , state_type &dxdt , double t )
+{
+    dxdt[0] = sigma * ( x[1] - x[0] );
+    dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
+    dxdt[2] = x[0]*x[1] - b * x[2];
+}
+
+int main()
+{
+    explicit_euler< state_type , initially_resizer > euler;
+    //you would use never_resizer for arrays, but just to show that no resizing is called
+    //even with initially_resizer I use this one here.
+
+    state_type x;
+    x[0] = 1.0; x[1] = 1.0; x[2] = 2.0;
+
+    euler.do_step( lorenz , x , 0.0 , 0.1 );
+
+    cout << x[0] << "  " << x[1] << "  " << x[2] << endl;
+}

Modified: sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_gsl_vector.cpp
==============================================================================
--- sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_gsl_vector.cpp	(original)
+++ sandbox/odeint/branches/karsten/libs/numeric/odeint/ideas/state_wrapper/test_gsl_vector.cpp	2011-06-26 10:03:29 EDT (Sun, 26 Jun 2011)
@@ -37,6 +37,19 @@
         gsl_vector_free( m_v );
     }
 
+    bool same_size( const gsl_vector *x )
+    {
+        return ( m_v->size == x->size );
+    }
+
+    void resize( const gsl_vector *x )
+    {
+        if( x->size == 0 ) return;
+
+        gsl_vector_free( m_v );
+        m_v = gsl_vector_alloc( x->size );
+    }
+
 };
 
 int main() {