Changeset 181 for bdm/stat/emix.h

Timestamp:

10/15/08 19:11:17 (16 years ago)

Author:

smidl

Message:

Regenerated doc

Files:

• bdm/stat/emix.h (modified) (5 diffs)

bdm/stat/emix.h

@@ -40 +40 @@
 public:
         //!Default constructor
-        emix (const RV &rv ) : epdf ( rv ) {};
+        emix ( const RV &rv ) : epdf ( rv ) {};
         //! Set weights \c w and components \c Coms , Coms are not copied!
         void set_parameters ( const vec &w, const Array<epdf*> &Coms, bool copy=true );
@@ -60 +60 @@
         //! returns a pointer to the internal mean value. Use with Care!
         vec& _w() {return w;}
-        virtual ~emix(){if (destroyComs){for(int i=0;i<Coms.length();i++){delete Coms(i);}}}
+        virtual ~emix() {if ( destroyComs ) {for ( int i=0;i<Coms.length();i++ ) {delete Coms ( i );}}}
         //! Auxiliary function for taking ownership of the Coms()
-        void ownComs(){destroyComs=true;}
+        void ownComs() {destroyComs=true;}
 };
 
@@ -75 +75 @@
 class mprod: public compositepdf, public mpdf {
 protected:
-        // pointers to epdfs
+        //! pointers to epdfs - shortcut to mpdfs()._epdf()
         Array<epdf*> epdfs;
         //! Indeces of rvc in common rvc
-        Array<ivec> rvcinds;
+        Array<ivec> rvcsinrvc;
+        //! Indeces of common rvc in rvcs
+        Array<ivec> rvcinrvcs;
 public:
         /*!\brief Constructor from list of mFacs,
-        Additional parameter overlap is left for future use. Do not set to true for mprod.
         */
-        mprod ( Array<mpdf*> mFacs): compositepdf(mFacs), mpdf(getrv(true),RV()), epdfs(n), rvcinds(n)
-        {       setrvc(rv,rvc);
-                setrvcinrv(rvc,rvcinds);
-                setindices(rv); 
-        for(int i=0;i<n;i++){epdfs(i)=&(mpdfs(i)->_epdf());}
-        };
-
-        double evalpdflog ( const vec &val ) const {
+        mprod ( Array<mpdf*> mFacs ) : compositepdf ( mFacs ), mpdf ( getrv ( true ),RV() ), epdfs ( n ), rvcsinrvc ( n ) {
+                setrvc ( rv,rvc );
+                setindices ( rv );
+                for ( int i = 0;i < n;i++ ) {
+                        // find ith rvc in common rv
+                        rvcsinrvc ( i ) = mpdfs ( i )->_rvc().dataind ( rvc );
+                        //and vice-versa
+                        rvcinrvcs ( i ) = rvc.dataind ( mpdfs ( i )->_rvc() );
+                }
+
+                for ( int i=0;i<n;i++ ) {
+                        epdfs ( i ) =& ( mpdfs ( i )->_epdf() );
+                }
+        };
+
+        double evalcond ( const vec &val, const vec &cond ) const {
                 int i;
                 double res = 0.0;
+                vec condi;
                 for ( i = n - 1;i > 0;i++ ) {
-                        if ( rvcinds ( i ).length() > 0 )
-                                {mpdfs ( i )->condition ( val ( rvcinds ( i ) ) );}
+                        if ( ( rvcsinrvc ( i ).length() > 0 ) || ( rvcsinrv ( i ).length() > 0 ) ) {
+                                condi = zeros ( rvcsinrvc ( i ).length() +rvcsinrv ( i ).length() );
+                                //copy part of the rvc into condi
+                                set_subvector ( condi, rvcinrvcs ( i ), get_vec ( cond,rvcsinrvc ( i ) ) );
+                                //copy part of the rv into condi
+                                set_subvector ( condi, rvinrvcs ( i ), get_vec ( val,rvcsinrv ( i ) ) );
+                                mpdfs ( i )->condition ( condi );
+                        }
                         // add logarithms
                         res += epdfs ( i )->evalpdflog ( val ( rvsinrv ( i ) ) );
@@ -107 +123 @@
                 ll = 0;
                 for ( int i = ( n - 1 );i >= 0;i-- ) {
-                        if (( rvcinds ( i ).length() > 0 )||( rvcsinrv ( i ).length() > 0 )) {
-                                condi = zeros ( rvcsinrv(i).length() + rvcinds(i).length() );
+                        if ( ( rvcsinrvc ( i ).length() > 0 ) || ( rvcsinrv ( i ).length() > 0 ) ) {
+                                condi = zeros ( rvcsinrv ( i ).length() + rvcsinrvc ( i ).length() );
                                 // copy data in condition
-                                set_subvector ( condi,rvcinds ( i ), cond );
+                                set_subvector ( condi,rvcsinrvc ( i ), cond );
                                 // copy data in already generated sample
-                                set_subvector ( condi,rvinrvcs ( i ), get_vec(smp,rvcsinrv(i)) );
+                                set_subvector ( condi,rvinrvcs ( i ), get_vec ( smp,rvcsinrv ( i ) ) );
 
                                 mpdfs ( i )->condition ( condi );
@@ -169 +185 @@
                 double tmp=0;
                 for ( int i=0;i<epdfs.length();i++ ) {
-                        tmp+=epdfs(i)->evalpdflog(val(rvinds(i)));
+                        tmp+=epdfs ( i )->evalpdflog ( val ( rvinds ( i ) ) );
                 }
                 return tmp;
         }
         //!access function
-        const epdf* operator () (int i) const {it_assert_debug(i<epdfs.length(),"wrong index");return epdfs(i);}
+        const epdf* operator () ( int i ) const {it_assert_debug ( i<epdfs.length(),"wrong index" );return epdfs ( i );}
 };