boost/odeint
2
0
Fork 0
mirror of https://github.com/boostorg/odeint.git synced 2026-01-26 06:42:23 +00:00
Code Issues Projects Releases Wiki Activity

fixed bug in check_implicit_euler

Browse Source
This commit is contained in:
mariomulansky
2010-09-29 17:39:09 +00:00
parent d67823ec8a
commit 12e81ad32f
1 changed files with 6 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
libs/numeric/odeint/test/check_implicit_euler.cpp
View File

@@ -38,7 +38,7 @@ void jacobi( state_type &x , matrix_type &jacobi , const value_type t )
jacobi( 0 , 0 ) = 1;
jacobi( 0 , 1 ) = 2;
jacobi( 1 , 0 ) = 0;
jacobi( 1 , 1 ) = 3;
jacobi( 1 , 1 ) = 1;
}
BOOST_AUTO_TEST_SUITE( implicit_euler_test )
@@ -53,7 +53,11 @@ BOOST_AUTO_TEST_CASE( test_euler )
stepper.do_step( sys , jacobi , x , 0.0 , 0.1 );
BOOST_CHECK_MESSAGE( abs( x(0) - 0.1 ) < eps , x[0] - 0.1 );
using std::abs;
// compare with analytic solution of above system
BOOST_CHECK_MESSAGE( abs( x(0) - 20.0/81.0 ) < eps , x(0) - 20.0/81.0 );
BOOST_CHECK_MESSAGE( abs( x(1) - 10.0/9.0 ) < eps , x(0) - 10.0/9.0 );
}