Changeset 180 for tests/merger_test.cpp

Show
Ignore:
Timestamp:
10/15/08 19:08:06 (16 years ago)
Author:
smidl
Message:

Modifications of BDM to reflect changes in basics

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • tests/merger_test.cpp

    r177 r180  
    1313        RV x ( "{x }","1" ); 
    1414 
     15        RV z(x); 
     16         
    1517        enorm<fsqmat> f1 ( x ); 
    1618        enorm<fsqmat> f2 ( x ); 
     
    2022 
    2123        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         
    2529        int Npoints=100; 
    2630        mat x_grid ( 1,Npoints ); 
     
    3539 
    3640        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         
    3747        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); 
    3951        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); 
    4156        vec m3 = M.lognorm_merge(lW); 
     57        M.merge(&g0); 
     58        vec mm3 = M.evalpdflog_m(x_grid); 
    4259         
    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                 
    4465        it_file it("merger_test.it"); 
    4566        it << Name("lf1") << l_f1; 
    4667        it << Name("lf2") << l_f2; 
    47         it << Name("m1") << m1; 
     68/*      it << Name("m1") << m1; 
     69        it << Name("mm1") << mm1;*/ 
    4870        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;*/ 
    5076} 
    51