Changeset 270 for bdm/estim/merger.h

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

  • bdm/estim/merger.h

    r262 r270  
    4949//!Default constructor 
    5050        merger ( const Array<mpdf*> &S ) : 
    51                         compositepdf ( S ), epdf ( getrv ( false ) ), 
    52                         Mix ( Array<BMEF*> ( 0 ),vec ( 0 ) ), dls ( n ), rvzs ( n ), zdls ( n ), eSmp(rv,0) { 
     51                        compositepdf ( S ), epdf ( ), 
     52                        Mix ( Array<BMEF*> ( 0 ),vec ( 0 ) ), dls ( n ), rvzs ( n ), zdls ( n ), eSmp() { 
    5353                RV ztmp; 
    5454                // Extend rv by rvc! 
    … …  
    9696                const Vec<double> &w = eSmp._w(); 
    9797                const Array<vec> &S = eSmp._samples(); 
    98                 vec tmp=zeros ( rv.count() );  
     98                vec tmp=zeros ( dim);  
    9999                for ( int i=0; i<Ns; i++ ) { 
    100100                        tmp+=w ( i ) *S ( i ); 
    … …  
    110110                cout << sum(w) << "," << w*w <<endl; 
    111111                 
    112                 mat Tmp=zeros(rv.count(), rv.count()); 
     112                mat Tmp=zeros(dim, dim); 
    113113                for ( int i=0; i<Ns; i++ ) { 
    114114                        Tmp+=w ( i ) *outer_product(S ( i ), S(i)); 
    … …  
    120120                const Array<vec> &S = eSmp._samples(); 
    121121                 
    122                 vec tmp=zeros(rv.count()); 
     122                vec tmp=zeros(dim); 
    123123                for ( int i=0; i<Ns; i++ ) { 
    124124                        tmp+=w ( i ) *pow(S ( i ),2);