Changeset 477 for library/tests/merger_test.cpp
- Timestamp:
- 08/05/09 14:40:03 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/tests/merger_test.cpp
r386 r477 13 13 RNG_randomize(); 14 14 15 RV x ( "{x }", "1" );15 RV x ( "{x }", "1" ); 16 16 17 17 RV z ( x ); 18 18 19 enorm<fsqmat> f1; f1.set_rv ( x ); 20 enorm<fsqmat> f2; f2.set_rv ( x ); 19 enorm<fsqmat> f1; 20 f1.set_rv ( x ); 21 enorm<fsqmat> f2; 22 f2.set_rv ( x ); 21 23 22 f1.set_parameters ( "-5", mat ( "2" ) );23 f2.set_parameters ( "5", mat ( "10" ) );24 f1.set_parameters ( "-5", mat ( "2" ) ); 25 f2.set_parameters ( "5", mat ( "10" ) ); 24 26 25 27 Array<mpdf* > A ( 2 ); 26 28 mepdf A1 ( &f1 ); 27 29 mepdf A2 ( &f2 ); 28 A ( 0 ) = &A1;29 A ( 1 ) = &A2;30 A ( 0 ) = &A1; 31 A ( 1 ) = &A2; 30 32 31 int Npoints =100;32 mat x_grid ( 1, Npoints );33 x_grid.set_row ( 0, linspace ( -10.0, 10.0 ) );33 int Npoints = 100; 34 mat x_grid ( 1, Npoints ); 35 x_grid.set_row ( 0, linspace ( -10.0, 10.0 ) ); 34 36 35 vec l_f1 =f1.evallog_m ( x_grid );36 vec l_f2 =f2.evallog_m ( x_grid );37 mat lW ( 2, Npoints );37 vec l_f1 = f1.evallog_m ( x_grid ); 38 vec l_f2 = f2.evallog_m ( x_grid ); 39 mat lW ( 2, Npoints ); 38 40 lW.set_row ( 0, l_f1 ); 39 41 lW.set_row ( 1, l_f2 ); 40 42 41 43 merger_base M ( A ); 42 enorm<fsqmat> g0; g0.set_rv ( x ); 43 g0.set_parameters ( vec ( "0.0" ),mat ( "100.0" ) ); 44 enorm<fsqmat> g0; 45 g0.set_rv ( x ); 46 g0.set_parameters ( vec ( "0.0" ), mat ( "100.0" ) ); 44 47 45 M.set_method ( 46 M.set_support (g0,200);48 M.set_method ( LOGNORMAL, 1.2 ); 49 M.set_support ( g0, 200 ); 47 50 M.merge (); 48 51