Changeset 124 for bdm/stat/emix.cpp

Timestamp:

06/20/08 11:16:12 (16 years ago)

Author:

smidl

Message:

zmeny v mmix

Files:

• bdm/stat/emix.cpp (modified) (1 diff)

bdm/stat/emix.cpp

@@ -3 +3 @@
 using namespace itpp;
 
-void emix::set_parameters ( const vec &w0, const Array<epdf*> &Coms0 ){
+void emix::set_parameters ( const vec &w0, const Array<epdf*> &Coms0 ) {
         w = w0;
         int i;
-        for(i=0;i<w.length();i++){it_assert_debug(rv.equal(Coms0(i)->_rv()),"RVs do not match!"); }
+        for ( i=0;i<w.length();i++ ) {
+                epdf* Mp=Coms0 ( i );
+                it_assert_debug ( rv.equal ( Coms0 ( i )->_rv() ),"RVs do not match!" );
+        }
         Coms = Coms0;
 }
 
-vec emix::sample() const{
-        //Sample which component 
+vec emix::sample() const {
+        //Sample which component
         vec cumDist = cumsum ( w );
         double u0 = UniRNG.sample();
 
         int i=0;
-        while ((w(i)<u0) && (i<(w.length()-1))){i++;}
-        
-        return Coms(i)->sample();
+        while ( ( w ( i ) <u0 ) && ( i< ( w.length()-1 ) ) ) {i++;}
+
+        return Coms ( i )->sample();
 }