Changeset 737 for library/bdm/stat/merger.cpp
- Timestamp:
- 11/25/09 12:14:38 (15 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
library/bdm/stat/merger.cpp
r713 r737 5 5 namespace bdm { 6 6 7 merger_base::merger_base ( const Array<shared_ptr<pdf> > &S ) :8 Npoints(0), DBG(false), dbg_file(0) {7 merger_base::merger_base ( const Array<shared_ptr<pdf> > &S ) : 8 Npoints ( 0 ), DBG ( false ), dbg_file ( 0 ) { 9 9 set_sources ( S ); 10 10 } … … 103 103 //Re-Initialize Mixture model 104 104 Mix.flatten ( &Mix_init ); 105 Mix.bayes_batch ( Smp_ex, empty_vec, w*Npoints );105 Mix.bayes_batch ( Smp_ex, empty_vec, w*Npoints ); 106 106 delete Mpred; 107 107 Mpred = Mix.epredictor ( ); // Allocation => must be deleted at the end!! … … 151 151 if ( pdfs ( i )->dimension() == dim ) { 152 152 // no need for conditioning or marginalization 153 lw_src = pdfs ( i )->evallogcond_mat ( Smp , vec (0));153 lw_src = pdfs ( i )->evallogcond_mat ( Smp , vec ( 0 ) ); 154 154 } else { 155 155 // compute likelihood of marginal on the conditional variable … … 187 187 // Compute likelihood of the partial source 188 188 for ( int k = 0; k < Npoints; k++ ) { 189 lw_src ( k ) += pdfs ( i )->evallogcond ( dls ( i )->pushdown ( Smp ( k ) ), 190 dls ( i )->get_cond ( Smp ( k ) ));189 lw_src ( k ) += pdfs ( i )->evallogcond ( dls ( i )->pushdown ( Smp ( k ) ), 190 dls ( i )->get_cond ( Smp ( k ) ) ); 191 191 } 192 192