Changeset 270 for tests/merger_test.cpp

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

Changes in the very root classes!
* rv and rvc are no longer compulsory,
* samplecond does not return ll
* BM has drv

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • tests/merger_test.cpp

    r262 r270  
    1212 
    1313        RNG_randomize(); 
    14          
     14 
    1515        RV x ( "{x }","1" ); 
    1616 
    17         RV z(x); 
    18          
    19         enorm<fsqmat> f1 ( x ); 
    20         enorm<fsqmat> f2 ( x ); 
     17        RV z ( x ); 
     18 
     19        enorm<fsqmat> f1; f1.set_rv ( x ); 
     20        enorm<fsqmat> f2; f2.set_rv ( x ); 
    2121 
    2222        f1.set_parameters ( "-5",mat ( "2" ) ); 
    … …  
    2424 
    2525        Array<mpdf* > A ( 2 ); 
    26         mepdf A1(&f1); 
    27         mepdf A2(&f2); 
     26        mepdf A1 ( &f1 ); 
     27        mepdf A2 ( &f2 ); 
    2828        A ( 0 ) =&A1; 
    2929        A ( 1 ) =&A2; 
    30          
     30 
    3131        int Npoints=100; 
    3232        mat x_grid ( 1,Npoints ); 
    … …  
    3636 
    3737        vec l_f2=f2.evallog_m ( x_grid ); 
    38         mat lW(2,Npoints); 
    39         lW.set_row(0, l_f1); 
    40         lW.set_row(1, l_f2 ); 
     38        mat lW ( 2,Npoints ); 
     39        lW.set_row ( 0, l_f1 ); 
     40        lW.set_row ( 1, l_f2 ); 
    4141 
    4242        merger M ( A ); 
    43         enorm<fsqmat> g0(x); 
    44         g0.set_parameters(vec("0.0"),mat("100.0")); 
    45          
    46         M.set_parameters(1.2,200,3); 
    47         M.merge(&g0); 
    48          
    49         vec m2 = M.lognorm_merge(lW); 
    50         vec mm2 = M.evallog_m(x_grid); 
    51          
    52         it_file it("merger_test.it"); 
    53         it << Name("x_grid") << x_grid; 
    54         it << Name("lf1") << l_f1; 
    55         it << Name("lf2") << l_f2; 
    56         it << Name("m2") << m2; 
    57         it << Name("mm2") << mm2; 
     43        enorm<fsqmat> g0; g0.set_rv ( x ); 
     44        g0.set_parameters ( vec ( "0.0" ),mat ( "100.0" ) ); 
     45 
     46        M.set_parameters ( 1.2,200,3 ); 
     47        M.merge ( &g0 ); 
     48 
     49        vec m2 = M.lognorm_merge ( lW ); 
     50        vec mm2 = M.evallog_m ( x_grid ); 
     51 
     52        it_file it ( "merger_test.it" ); 
     53        it << Name ( "x_grid" ) << x_grid; 
     54        it << Name ( "lf1" ) << l_f1; 
     55        it << Name ( "lf2" ) << l_f2; 
     56        it << Name ( "m2" ) << m2; 
     57        it << Name ( "mm2" ) << mm2; 
    5858}