Show
Ignore:
Timestamp:
08/05/09 14:40:03 (15 years ago)
Author:
mido
Message:

panove, vite, jak jsem peclivej na upravu kodu.. snad se vam bude libit:) konfigurace je v souboru /system/astylerc

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • library/tests/merger_iter_test.cpp

    r386 r477  
    1414 
    1515        RNG_randomize(); 
    16          
    17         RV x ( "{x }","1" ); 
    18         RV y ( "{y }","1" ); 
    1916 
    20         RV xy=x; xy.add(y); 
    21          
    22         enorm<fsqmat> f1;f1.set_rv ( xy ); 
    23         enorm<fsqmat> f2;f2.set_rv ( xy ); 
    24         enorm<fsqmat> f3;f3.set_rv(y); 
    25                  
    26         f1.set_parameters ( "4 3",mat ( "0.4 0.3; 0.3 0.4" ) ); 
    27         f2.set_parameters ( "1 3",mat ( "0.3 -0.2; -0.2 0.3" ) ); 
    28         f3.set_parameters ( "2",mat("0.4") ); 
    29          
     17        RV x ( "{x }", "1" ); 
     18        RV y ( "{y }", "1" ); 
     19 
     20        RV xy = x; 
     21        xy.add ( y ); 
     22 
     23        enorm<fsqmat> f1; 
     24        f1.set_rv ( xy ); 
     25        enorm<fsqmat> f2; 
     26        f2.set_rv ( xy ); 
     27        enorm<fsqmat> f3; 
     28        f3.set_rv ( y ); 
     29 
     30        f1.set_parameters ( "4 3", mat ( "0.4 0.3; 0.3 0.4" ) ); 
     31        f2.set_parameters ( "1 3", mat ( "0.3 -0.2; -0.2 0.3" ) ); 
     32        f3.set_parameters ( "2", mat ( "0.4" ) ); 
     33 
    3034        Array<mpdf* > A ( 3 ); 
    31         mepdf A1(&f1); 
    32         mepdf A2(&f2); 
    33         mepdf A3(&f3); 
    34         A ( 0 ) =&A1; 
    35         A ( 1 ) =&A2; 
    36         A ( 2 ) =&A3; 
    37          
    38         int Npoints=100; 
    39         mat x_grid(1,Npoints); 
    40         x_grid.set_row(0,linspace ( -2.0, 4.0, Npoints )); 
    41         mat y_grid(1,Npoints); 
    42         y_grid.set_row(0,linspace ( -2.0, 4.0, Npoints )); 
    43          
    44         mat Grid (2,Npoints*Npoints); 
    45         Grid.set_submatrix ( 0,0, kron(x_grid,ones(1,Npoints)) ); 
    46         Grid.set_submatrix ( 1,0, kron(ones(1,Npoints), y_grid) ); 
    47          
     35        mepdf A1 ( &f1 ); 
     36        mepdf A2 ( &f2 ); 
     37        mepdf A3 ( &f3 ); 
     38        A ( 0 ) = &A1; 
     39        A ( 1 ) = &A2; 
     40        A ( 2 ) = &A3; 
     41 
     42        int Npoints = 100; 
     43        mat x_grid ( 1, Npoints ); 
     44        x_grid.set_row ( 0, linspace ( -2.0, 4.0, Npoints ) ); 
     45        mat y_grid ( 1, Npoints ); 
     46        y_grid.set_row ( 0, linspace ( -2.0, 4.0, Npoints ) ); 
     47 
     48        mat Grid ( 2, Npoints*Npoints ); 
     49        Grid.set_submatrix ( 0, 0, kron ( x_grid, ones ( 1, Npoints ) ) ); 
     50        Grid.set_submatrix ( 1, 0, kron ( ones ( 1, Npoints ), y_grid ) ); 
     51 
    4852        merger_mix M ( A ); 
    49         enorm<fsqmat> g0; g0.set_rv(xy); 
    50         g0.set_parameters(vec("4 4"),mat("1 0; 0 1")); 
    51          
    52         M.set_parameters(5); 
    53         M.set_method(LOGNORMAL,1.2); 
    54         M.set_support(g0,400); 
     53        enorm<fsqmat> g0; 
     54        g0.set_rv ( xy ); 
     55        g0.set_parameters ( vec ( "4 4" ), mat ( "1 0; 0 1" ) ); 
     56 
     57        M.set_parameters ( 5 ); 
     58        M.set_method ( LOGNORMAL, 1.2 ); 
     59        M.set_support ( g0, 400 ); 
    5560        M.merge(); 
    56          
     61 
    5762        MixEF &MM = M._Mix(); 
    5863        epdf* MP = MM.epredictor();//xy 
    59          
    60         vec Res1 = M.evallog_m(Grid); 
    61         mat Res2 = ((emix*)MP)->evallog_M(Grid); 
    62          
    63         it_file it("merger_iter_test.it"); 
    64         it << Name("Npoints") << Npoints; 
    65         it << Name("Grid") << Grid; 
    66         it << Name("Res1") << Res1; 
    67         it << Name("Res2") << Res2; 
     64 
     65        vec Res1 = M.evallog_m ( Grid ); 
     66        mat Res2 = ( ( emix* ) MP )->evallog_M ( Grid ); 
     67 
     68        it_file it ( "merger_iter_test.it" ); 
     69        it << Name ( "Npoints" ) << Npoints; 
     70        it << Name ( "Grid" ) << Grid; 
     71        it << Name ( "Res1" ) << Res1; 
     72        it << Name ( "Res2" ) << Res2; 
    6873}