Changeset 395 for library/bdm/stat/emix.h

Timestamp:

06/22/09 13:17:49 (15 years ago)

Author:

smidl

Message:

merging works for merger_mx

Files:

• library/bdm/stat/emix.h (modified) (3 diffs)

library/bdm/stat/emix.h

@@ -235 +235 @@
         mprod (){};
         mprod (Array<mpdf*> mFacs ){set_elements( mFacs );};
-        void set_elements(Array<mpdf*> mFacs ) {
+        void set_elements(Array<mpdf*> mFacs , bool own=false) {
                 
-                set_elements(mFacs);
-                
+                compositepdf::set_elements(mFacs,own);
+                dls.set_size(mFacs.length());
+                epdfs.set_size(mFacs.length());
+                                
                 ep=&dummy;
                 RV rv=getrv ( true );
-                set_rv ( rv );dummy.set_parameters ( rv._dsize() );
+                set_rv ( rv );
+                dummy.set_parameters ( rv._dsize() );
                 setrvc ( ep->_rv(),rvc );
                 // rv and rvc established = > we can link them with mpdfs
@@ -270 +273 @@
                 return res;
         }
+        vec evallogcond_m(const mat &Dt, const vec &cond) {
+                vec tmp(Dt.cols());
+                for(int i=0;i<Dt.cols(); i++){
+                        tmp(i) = evallogcond(Dt.get_col(i),cond);
+                }
+                return tmp;
+        };
+        vec evallogcond_m(const Array<vec> &Dt, const vec &cond) {
+                vec tmp(Dt.length());
+                for(int i=0;i<Dt.length(); i++){
+                        tmp(i) = evallogcond(Dt(i),cond);
+                }
+                return tmp;             
+        };
+
+
         //TODO smarter...
         vec samplecond ( const vec &cond ) {
@@ -294 +313 @@
 
         ~mprod() {};
-};
+        //! Load from structure with elements:
+        //!  \code
+        //! { class='mprod';
+        //!   mpdfs = (..., ...);     // list of mpdfs in the order of chain rule
+        //! }
+        //! \endcode
+        //!@}
+        void from_setting(const Setting &set){
+                Array<mpdf*> Atmp; //temporary Array
+                UI::get(Atmp,set, "mpdfs");
+                set_elements(Atmp,true);
+        }
+        
+};
+UIREGISTER(mprod);
 
 //! Product of independent epdfs. For dependent pdfs, use mprod.