Changeset 488 for library/bdm/stat/emix.cpp

Timestamp:

08/08/09 13:43:13 (15 years ago)

Author:

smidl

Message:

changes in mpdf -> compile OK, broken tests!

Files:

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

library/bdm/stat/emix.cpp

@@ -199 +199 @@
         return aprgiw;
 };
+
+vec mmix::samplecond(const vec &cond) {
+        //Sample which component
+        vec cumDist = cumsum ( w );
+        double u0;
+#pragma omp critical
+        u0 = UniRNG.sample();
+
+        int i = 0;
+        while ( ( cumDist ( i ) < u0 ) && ( i < ( w.length() - 1 ) ) ) {
+                i++;
+        }
+
+        return Coms ( i )->samplecond(cond);
+}
 
 }