Changeset 180 for tests/merger_test.cpp
- Timestamp:
- 10/15/08 19:08:06 (16 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
tests/merger_test.cpp
r177 r180 13 13 RV x ( "{x }","1" ); 14 14 15 RV z(x); 16 15 17 enorm<fsqmat> f1 ( x ); 16 18 enorm<fsqmat> f2 ( x ); … … 20 22 21 23 Array<mpdf* > A ( 2 ); 22 A ( 0 ) =&mepdf(f1); 23 A ( 1 ) =&mepdf(f2); 24 24 mepdf A1(f1); 25 mepdf A2(f2); 26 A ( 0 ) =&A1; 27 A ( 1 ) =&A2; 28 25 29 int Npoints=100; 26 30 mat x_grid ( 1,Npoints ); … … 35 39 36 40 merger M ( A ); 41 enorm<fsqmat> g0(x); 42 g0.set_parameters(vec("0.0"),mat("100.0")); 43 44 M.set_parameters(1.2,1000,3); 45 M.merge(&g0); 46 37 47 vec m2 = M.lognorm_merge(lW); 38 M.set_parameters(31.0/30.0,10,100); 48 vec mm2 = M.evalpdflog_m(x_grid); 49 50 /* M.set_parameters(1.0,10,3); 39 51 vec m1 = M.lognorm_merge(lW); 40 M.set_parameters(3.0,10,100); 52 M.merge(&g0); 53 vec mm1 = M.evalpdflog_m(x_grid); 54 55 M.set_parameters(3.0,10,3); 41 56 vec m3 = M.lognorm_merge(lW); 57 M.merge(&g0); 58 vec mm3 = M.evalpdflog_m(x_grid); 42 59 43 60 M.set_parameters(300.0,10,3); 61 vec m300 = M.lognorm_merge(lW); 62 M.merge(&g0); 63 vec mm300 = M.evalpdflog_m(x_grid);*/ 64 44 65 it_file it("merger_test.it"); 45 66 it << Name("lf1") << l_f1; 46 67 it << Name("lf2") << l_f2; 47 it << Name("m1") << m1; 68 /* it << Name("m1") << m1; 69 it << Name("mm1") << mm1;*/ 48 70 it << Name("m2") << m2; 49 it << Name("m3") << m3; 71 it << Name("mm2") << mm2; 72 /* it << Name("m3") << m3; 73 it << Name("mm3") << mm3; 74 it << Name("m300") << m300; 75 it << Name("mm300") << mm300;*/ 50 76 } 51