Changeset 280 for bdm/stat/libEF.h

Timestamp:

02/24/09 14:12:46 (15 years ago)

Author:

smidl

Message:

progress...

Files:

• bdm/stat/libEF.h (modified) (6 diffs)

bdm/stat/libEF.h

@@ -116 +116 @@
         //!@{
 
-        enorm ( ):eEF ( ), mu ( ),R ( ) {};
+        enorm ( ) :eEF ( ), mu ( ),R ( ) {};
         enorm ( const vec &mu,const sq_T &R ) {set_parameters ( mu,R );}
         void set_parameters ( const vec &mu,const sq_T &R );
@@ -366 +366 @@
         //!@{
         eigamma ( ) :eEF ( ), eg(),alpha ( eg._alpha() ), beta ( eg._beta() ) {};
-        eigamma ( const vec &a, const vec &b ):eEF ( ), eg(),alpha ( eg._alpha() ), beta ( eg._beta() ) {eg.set_parameters ( a,b );};
+        eigamma ( const vec &a, const vec &b ) :eEF ( ), eg(),alpha ( eg._alpha() ), beta ( eg._beta() ) {eg.set_parameters ( a,b );};
         void set_parameters ( const vec &a, const vec &b ) {eg.set_parameters ( a,b );};
         //!@}
@@ -476 +476 @@
         template<class sq_M>
         friend std::ostream &operator<< ( std::ostream &os,  mlnorm<sq_M> &ml );
+};
+
+//! Mpdf with general function for mean value
+template<class sq_T>
+class mgnorm : public mEF {
+protected:
+        //! Internal epdf that arise by conditioning on \c rvc
+        enorm<sq_T> epdf;
+        vec &mu;
+        fnc* g;
+public:
+        //!default constructor
+        mgnorm() :mu ( epdf._mu() ) {ep=&epdf;}
+        //!set mean function
+        void set_parameters ( fnc* g0, const sq_T &R0 ) {g=g0; epdf.set_parameters ( zeros ( g->_dimy() ), R0 );}
+        void condition ( const vec &cond ) {mu=g->eval ( cond );};
 };
 
@@ -663 +679 @@
         Array<vec> samples;
 public:
-        //! \name Constructors 
-        //!@{
-        eEmp ( ) :epdf ( ),w (  ),samples ( ) {};
-        eEmp (const eEmp &e ): epdf(e), w(e.w), samples(e.samples) {};
-        //!@}
-        
+        //! \name Constructors
+        //!@{
+        eEmp ( ) :epdf ( ),w ( ),samples ( ) {};
+        eEmp ( const eEmp &e ) : epdf ( e ), w ( e.w ), samples ( e.samples ) {};
+        //!@}
+
         //! Set samples and weights
         void set_parameters ( const vec &w0, const epdf* pdf0 );
@@ -817 +833 @@
 template<class sq_T>
 epdf* enorm<sq_T>::marginal ( const RV &rvn ) const {
-        it_assert_debug(isnamed(), "rv description is not assigned");
+        it_assert_debug ( isnamed(), "rv description is not assigned" );
         ivec irvn = rvn.dataind ( rv );
 
         sq_T Rn ( R,irvn ); //select rows and columns of R
-        
+
         enorm<sq_T>* tmp = new enorm<sq_T>;
         tmp->set_rv ( rvn );
@@ -856 +872 @@
 
         mlnorm<sq_T>* tmp=new mlnorm<sq_T> ( );
-        tmp->set_rv(rvn); tmp->set_rvc(rvc);
+        tmp->set_rv ( rvn ); tmp->set_rvc ( rvc );
         tmp->set_parameters ( A,mu1-A*mu2,R_n );
         return tmp;