Changeset 384 for library/bdm/stat/mixtures.h

Timestamp:

06/19/09 10:17:25 (15 years ago)

Author:

mido

Message:

possibly broken?

Location:

library/bdm

Files:

• . (moved) (moved from bdm)
• stat/mixtures.h (moved) (moved from bdm/stat/emix.h) (5 diffs)

library/bdm/stat/mixtures.h

@@ -11 +11 @@
 */
 
-#ifndef MX_H
-#define MX_H
+#ifndef STAT_MIXTURES_H
+#define STAT_MIXTURES_H
 
 #define LOG2  0.69314718055995  
 
-#include "libEF.h"
-//#include <std>
+#include "exp_family.h"
 
 namespace bdm {
@@ -234 +233 @@
         /*!\brief Constructor from list of mFacs,
         */
-        mprod ( Array<mpdf*> mFacs ) : compositepdf ( mFacs ), mpdf (), epdfs ( n ), dls ( n ) {
+        mprod (){};
+        mprod (Array<mpdf*> mFacs ){set_elements( mFacs );};
+        void set_elements(Array<mpdf*> mFacs ) {
+                
+                set_elements(mFacs);
+                
                 ep=&dummy;
                 RV rv=getrv ( true );
@@ -240 +244 @@
                 setrvc ( ep->_rv(),rvc );
                 // rv and rvc established = > we can link them with mpdfs
-                for ( int i = 0;i < n;i++ ) {
+                for ( int i = 0;i < mpdfs.length(); i++ ) {
                         dls ( i ) = new datalink_m2m;
                         dls(i)->set_connection( mpdfs ( i )->_rv(), mpdfs ( i )->_rvc(), _rv(), _rvc() );
                 }
 
-                for ( int i=0;i<n;i++ ) {
+                for ( int i=0; i<mpdfs.length(); i++ ) {
                         epdfs ( i ) =& ( mpdfs ( i )->_epdf() );
                 }
@@ -253 +257 @@
                 int i;
                 double res = 0.0;
-                for ( i = n - 1;i >= 0;i-- ) {
+                for ( i = mpdfs.length() - 1;i >= 0;i-- ) {
                         /*                      if ( mpdfs(i)->_rvc().count() >0) {
                                                         mpdfs ( i )->condition ( dls ( i )->get_cond ( val,cond ) );
@@ -272 +276 @@
                 vec smpi;
                 // Hard assumption here!!! We are going backwards, to assure that samples that are needed from smp are already generated!
-                for ( int i = ( n - 1 );i >= 0;i-- ) {
+                for ( int i = ( mpdfs.length() - 1 );i >= 0;i-- ) {
                         if ( mpdfs ( i )->dimensionc() ) {
                                 mpdfs ( i )->condition ( dls ( i )->get_cond ( smp ,cond ) ); // smp is val here!!