Changeset 956

Timestamp:

05/19/10 03:56:03 (15 years ago)

Author:

sarka

Message:

to_setting

Location:

library/bdm

Files:

• base/bdmbase.cpp (modified) (1 diff)
• base/bdmbase.h (modified) (2 diffs)
• base/itpp/base/itfile.h (modified) (1 diff)
• stat/emix.cpp (modified) (5 diffs)
• stat/emix.h (modified) (4 diffs)
• stat/exp_family.cpp (modified) (5 diffs)
• stat/exp_family.h (modified) (27 diffs)
• stat/merger.cpp (modified) (2 diffs)
• stat/merger.h (modified) (2 diffs)

library/bdm/base/bdmbase.cpp

@@ -460 +460 @@
 }
 
+void epdf::from_setting ( const Setting &set ) {
+                root::from_setting( set );
+                UI::get(rv,  set, "rv", UI::optional );
+        }
+
+void epdf::to_setting ( Setting &set ) const {
+                // we do not want to store rvc, therfore, pdf::to_setting( set ) is omitted
+                root::to_setting(set);
+
+                UI::save( &rv, set, "rv" );
+        }
 
 

library/bdm/base/bdmbase.h

@@ -761 +761 @@
         //! \endcode
         //!@}
-        void from_setting ( const Setting &set ) {
-                root::from_setting( set );
-                UI::get(rv,  set, "rv", UI::optional );
-        }
-
-        void to_setting ( Setting &set ) const {
-                // we do not want to store rvc, therfore, pdf::to_setting( set ) is omitted
-                root::to_setting(set);
-
-                UI::save( &rv, set, "rv" );
-        }
+        void from_setting ( const Setting &set );
+        void to_setting ( Setting &set ) const;
 
         vec samplecond ( const vec &cond ) {
@@ -814 +805 @@
 
         void validate() {
+                pdf::validate();
                 iepdf.validate();
                 if ( rv._dsize() < iepdf._rv()._dsize() ) {

library/bdm/base/itpp/base/itfile.h

@@ -367 +367 @@
   //! Remove slack space from the file
   void pack();
+
+  std::string get_fname(){
+                return fname;
+  }
 
 protected:

library/bdm/stat/emix.cpp

@@ -4 +4 @@
 
 void emix_base::validate (){
+        epdf::validate();
         bdm_assert ( no_coms() > 0, "There has to be at least one component." );
 
@@ -116 +117 @@
 
 void emix::from_setting ( const Setting &set ) {
+        emix_base::from_setting(set);
         UI::get ( Coms, set, "pdfs", UI::compulsory );
         UI::get ( w, set, "weights", UI::compulsory );
 }
+void emix::to_setting  (Setting  &set) const {
+        emix_base::to_setting(set);
+        UI::save(Coms, set, "pdfs");
+        UI::save( w, set, "weights");
+}
+
 
 void    emix::validate (){
@@ -179 +187 @@
 }
 
+void mprod::from_setting ( const Setting &set ) {
+                pdf::from_setting(set);
+                Array<shared_ptr<pdf> > temp_array; 
+                UI::get ( temp_array, set, "pdfs", UI::compulsory );
+                set_elements ( temp_array );
+        }
+void    mprod::to_setting  (Setting  &set) const {
+                pdf::to_setting(set);
+                UI::save( pdfs, set, "pdfs");
+        }
+
 void mmix::validate()
-{       
+{       pdf::validate();
         bdm_assert ( Coms.length() > 0, "There has to be at least one component." );
 
@@ -208 +227 @@
 
 void mmix::from_setting ( const Setting &set ) {
+        
+        pdf::from_setting(set);
         UI::get ( Coms, set, "pdfs", UI::compulsory );
 
@@ -215 +236 @@
                 w = 1.0 / len;
         }
+}
+
+void    mmix::to_setting  (Setting  &set) const {
+        pdf::to_setting(set);
+        UI::save( Coms, set, "pdfs");
+        UI::save( w, set, "weights");
 }
 

library/bdm/stat/emix.h

@@ -191 +191 @@
         //!@}
         void from_setting ( const Setting &set );
+        void to_setting  (Setting  &set) const;
 
         void set_rv ( const RV &rv ) {
@@ -262 +263 @@
         //! \endcode
         //!@}
-        void from_setting ( const Setting &set ) {
-                Array<shared_ptr<pdf> > temp_array; 
-                UI::get ( temp_array, set, "pdfs", UI::compulsory );
-                set_elements ( temp_array );
-        }
+        void from_setting ( const Setting &set ) ;
+        void    to_setting  (Setting  &set) const;
+
+
 };
 UIREGISTER ( mprod );
@@ -300 +300 @@
         }
         void validate() {
+                epdf::validate();
                 dls.set_length ( no_factors() );
                 
@@ -387 +388 @@
         //!@}
         void from_setting ( const Setting &set );
-
+        void    to_setting  (Setting  &set) const;
         virtual void validate();
 };

library/bdm/stat/exp_family.cpp

@@ -326 +326 @@
 }
 
+void egiw::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set );
+                UI::get ( dimx, set, "dimx", UI::compulsory );
+                if ( !UI::get ( nu, set, "nu", UI::optional ) ) {
+                        nu = -1;
+                }
+                mat V;
+                if ( !UI::get ( V, set, "V", UI::optional ) ) {
+                        vec dV;
+                        UI::get ( dV, set, "dV", UI::compulsory );
+                        set_parameters ( dimx, ldmat ( dV ), nu );
+
+                } else {
+                        set_parameters ( dimx, V, nu );
+                }
+        }
+
+void egiw::to_setting ( Setting& set ) const {
+                epdf::to_setting ( set );
+                UI::save ( dimx, set, "dimx" );
+                UI::save ( V.to_mat(), set, "V" );
+                UI::save ( nu, set, "nu" );
+        };
+
+void egiw::validate() {
+                eEF::validate();
+                dim = dimx * ( dimx + nPsi );   
+        
+            // check sizes, rvs etc.
+                // also check if RV are meaningful!!!
+                // meaningful =  rv for theta  and rv for r are split!
+        }
 void multiBM::bayes ( const vec &yt, const vec &cond ) {
         if ( frg < 1.0 ) {
@@ -421 +453 @@
 }
 
+void egamma::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set ); // reads rv
+                UI::get ( alpha, set, "alpha", UI::compulsory );
+                UI::get ( beta, set, "beta", UI::compulsory );
+        }
+        
+void egamma::to_setting ( Setting &set ) const
+        {                       
+                epdf::to_setting( set );
+                UI::save( alpha, set, "alpha" );
+                UI::save( beta, set, "beta" );
+        } 
+        
+        
+void egamma::validate() {
+                eEF::validate();
+                bdm_assert ( alpha.length() == beta.length(), "parameters do not match" );
+                dim = alpha.length();
+        }
+
 void mgamma::set_parameters ( double k0, const vec &beta0 ) {
         k = k0;
         iepdf.set_parameters ( k * ones ( beta0.length() ), beta0 );
 }
+
+
+void mgamma::from_setting ( const Setting &set ) {
+                pdf::from_setting ( set ); // reads rv and rvc
+                vec betatmp; // ugly but necessary
+                UI::get ( betatmp, set, "beta", UI::compulsory );
+                UI::get ( k, set, "k", UI::compulsory );
+                set_parameters ( k, betatmp );
+        }
+
+void mgamma::to_setting  (Setting  &set) const {
+                pdf::to_setting(set);
+                UI::save( _beta, set, "beta");
+                UI::save( k, set, "k");
+
+        }
+
+void mgamma::validate() {
+                pdf_internal<egamma>::validate();
+
+                dim = _beta.length();
+                dimc = _beta.length();
+        }
 
 void eEmp::resample ( RESAMPLING_METHOD method ) {
@@ -536 +611 @@
 
 void migamma_ref::from_setting ( const Setting &set ) {
+        migamma::from_setting(set);
         vec ref;
         UI::get ( ref, set, "ref" , UI::compulsory );
@@ -541 +617 @@
         validate();
 }
-
+void  migamma_ref::to_setting  (Setting  &set) const {
+        migamma::to_setting(set);
+        UI::save ( pow ( refl, 1/(1.0 - l) ), set, "ref");      
+        UI::save(l,set,"l");
+        UI::save(k,set,"k");
+}
 void mlognorm::from_setting ( const Setting &set ) {
+        pdf_internal<elognorm>::from_setting(set);
         vec mu0;
         UI::get ( mu0, set, "mu0", UI::compulsory );
         set_parameters ( mu0.length(), set["k"] );
-        validate();
         condition ( mu0 );
 }
+
+
+void mlognorm::to_setting  (Setting  &set) const {
+        pdf_internal<elognorm>::to_setting(set);
+        UI::save ( exp(mu + sig2), set, "mu0");
+}
+
 
 void mlstudent::condition ( const vec &cond ) {
@@ -586 +674 @@
 }
 
+void eEmp::to_setting ( Setting &set ) const {
+                epdf::to_setting( set );
+                UI::save ( samples, set, "samples" );
+                UI::save ( w, set, "w" );
+        }
+
+void eEmp::from_setting ( const Setting &set ) {
+                epdf::from_setting( set );
+                
+                UI::get( samples, set, "samples", UI::compulsory );
+                UI::get ( w, set, "w", UI::compulsory );
+        }
+
+void    eEmp::validate (){
+          epdf::validate();
+          bdm_assert (samples.length()==w.length(),"samples and weigths are of different lengths");
+          n = w.length();
+          if (n>0)
+                pdf::dim = samples ( 0 ).length();
+        }
+        
+        void eDirich::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set );
+                UI::get ( beta, set, "beta", UI::compulsory );
+        }
+void eDirich::validate() {
+                //check rv
+                eEF::validate();
+                dim = beta.length();
+        }
+
+void eDirich::to_setting ( Setting &set ) const
+        {                       
+                eEF::to_setting( set );
+                UI::save( beta, set, "beta" );
+        } 
+
+void euni::from_setting ( const Setting &set ) {
+                epdf::from_setting ( set ); // reads rv and rvc
+
+                UI::get ( high, set, "high", UI::compulsory );
+                UI::get ( low, set, "low", UI::compulsory );
+                set_parameters ( low, high );
+                
+        }
+        
+void    euni::to_setting  (Setting  &set) const {
+                epdf::to_setting ( set );
+                UI::save ( high, set, "high" );
+                UI::save ( low, set, "low" );
+        }
+        
+void euni::validate() {
+                epdf::validate();
+                bdm_assert ( high.length() == low.length(), "Incompatible high and low vectors" );
+                dim = high.length();
+                bdm_assert ( min ( distance ) > 0.0, "bad support" );
+        }
+
+void mgdirac::from_setting(const Setting& set){
+                        pdf::from_setting(set);
+                        g=UI::build<fnc>(set,"g",UI::compulsory);
+                        validate();
+                }
+void mgdirac::to_setting(Setting &set) const{
+                        pdf::to_setting(set);
+                        UI::save(g.get(), set, "g");
+                }
+void mgdirac::validate() {
+                        pdf::validate();
+                        dim = g->dimension();
+                        dimc = g->dimensionc();
+                }
+
+void mDirich::from_setting ( const Setting &set ) {
+                pdf::from_setting ( set ); // reads rv and rvc
+                if ( _rv()._dsize() > 0 ) {
+                        rvc = _rv().copy_t ( -1 );
+                }
+                vec beta0;
+                if ( !UI::get ( beta0, set, "beta0", UI::optional ) ) {
+                        beta0 = ones ( _rv()._dsize() );
+                }
+                if ( !UI::get ( betac, set, "betac", UI::optional ) ) {
+                        betac = 0.1 * ones ( _rv()._dsize() );
+                }
+                _beta = beta0;
+
+                UI::get ( k, set, "k", UI::compulsory );
+        }
+
+void mDirich::to_setting  (Setting  &set) const {
+                pdf::to_setting(set);
+                UI::save( _beta, set, "beta0");
+                UI::save( betac, set, "betac");
+                UI::save ( k, set, "k" );
+        }
+
+
+void mDirich::validate() {
+                pdf_internal<eDirich>::validate();
+                bdm_assert ( _beta.length() == betac.length(), "beta0 and betac are not compatible" );
+                if ( _rv()._dsize() > 0 ) {
+                        bdm_assert ( ( _rv()._dsize() == dimension() ) , "Size of rv does not match with beta" );
+                }
+                dimc = _beta.length();
+        }
 
 };

library/bdm/stat/exp_family.h

@@ -133 +133 @@
                         return std::numeric_limits< double >::max();
                 }
-                void from_setting(const Setting& set){
-                        pdf::from_setting(set);
-                        g=UI::build<fnc>(set,"g",UI::compulsory);
-                        validate();
-                }
-                void to_setting(Setting &set) const{
-                        pdf::to_setting(set);
-                        UI::save(g.get(), set, "g");
-                }
-                void validate() {
-                        dim = g->dimension();
-                        dimc = g->dimensionc();
-                }
+                void from_setting(const Setting& set);
+                void to_setting(Setting &set) const;
+                void validate();
 };
 UIREGISTER(mgdirac);
@@ -186 +176 @@
         void to_setting ( Setting &root ) const ;
         
-        void validate() {
-                bdm_assert ( mu.length() == R.rows(), "mu and R parameters do not match" );
-                dim = mu.length();
-        }
+        void validate();
         //!@}
 
@@ -336 +323 @@
                         }
                         void validate() {
+                                eEF::validate();
                                 dim = H.rows();
                         }
@@ -460 +448 @@
         \sa log_level_enums
         */
-        void from_setting ( const Setting &set ) {
-                epdf::from_setting ( set );
-                UI::get ( dimx, set, "dimx", UI::compulsory );
-                if ( !UI::get ( nu, set, "nu", UI::optional ) ) {
-                        nu = -1;
-                }
-                mat V;
-                if ( !UI::get ( V, set, "V", UI::optional ) ) {
-                        vec dV;
-                        UI::get ( dV, set, "dV", UI::compulsory );
-                        set_parameters ( dimx, ldmat ( dV ), nu );
-                        validate();
-
-                } else {
-                        set_parameters ( dimx, V, nu );
-                        validate();
-                }
-        }
-
-        void to_setting ( Setting& set ) const {
-                epdf::to_setting ( set );
-                UI::save ( dimx, set, "dimx" );
-                UI::save ( V.to_mat(), set, "V" );
-                UI::save ( nu, set, "nu" );
-        };
-
-        void validate() {
-        dim = dimx * ( dimx + nPsi );   
-        
-            // check sizes, rvs etc.
-                // also check if RV are meaningful!!!
-                // meaningful =  rv for theta  and rv for r are split!
-        }
+        void from_setting ( const Setting &set );
+        void to_setting ( Setting& set ) const;
+        void validate();
         void log_register ( bdm::logger& L, const string& prefix );
 
@@ -578 +536 @@
         \endcode
         */
-        void from_setting ( const Setting &set ) {
-                epdf::from_setting ( set );
-                UI::get ( beta, set, "beta", UI::compulsory );
-                validate();
-        }
-        void validate() {
-                //check rv
-                dim = beta.length();
-        }
-
-        void to_setting ( Setting &set ) const
-        {                       
-                eEF::to_setting( set );
-                UI::save( beta, set, "beta" );
-        } 
+        void from_setting ( const Setting &set );
+        void validate();
+        void to_setting ( Setting &set ) const;
 };
 UIREGISTER ( eDirich );
@@ -632 +578 @@
         \endcode
         */
-        void from_setting ( const Setting &set ) {
-                pdf::from_setting ( set ); // reads rv and rvc
-                if ( _rv()._dsize() > 0 ) {
-                        rvc = _rv().copy_t ( -1 );
-                }
-                vec beta0;
-                if ( !UI::get ( beta0, set, "beta0", UI::optional ) ) {
-                        beta0 = ones ( _rv()._dsize() );
-                }
-                if ( !UI::get ( betac, set, "betac", UI::optional ) ) {
-                        betac = 0.1 * ones ( _rv()._dsize() );
-                }
-                _beta = beta0;
-
-                UI::get ( k, set, "k", UI::compulsory );
-                validate();
-        }
-        void validate() {
-                pdf_internal<eDirich>::validate();
-                bdm_assert ( _beta.length() == betac.length(), "beta0 and betac are not compatible" );
-                if ( _rv()._dsize() > 0 ) {
-                        bdm_assert ( ( _rv()._dsize() == dimension() ) , "Size of rv does not match with beta" );
-                }
-                dimc = _beta.length();
-        };
+        void from_setting ( const Setting &set );
+        void    to_setting  (Setting  &set) const;
+        void validate();
 };
 UIREGISTER ( mDirich );
@@ -765 +689 @@
         \endcode
         */
-        void from_setting ( const Setting &set ) {
-                epdf::from_setting ( set ); // reads rv
-                UI::get ( alpha, set, "alpha", UI::compulsory );
-                UI::get ( beta, set, "beta", UI::compulsory );
-                validate();
-        }
-        void validate() {
-                bdm_assert ( alpha.length() == beta.length(), "parameters do not match" );
-                dim = alpha.length();
-        }
+        void from_setting ( const Setting &set );
+        void to_setting ( Setting &set ) const;
+        void validate(); 
 };
 UIREGISTER ( egamma );
@@ -888 +805 @@
          \endcode
          */
-        void from_setting ( const Setting &set ) {
-                epdf::from_setting ( set ); // reads rv and rvc
-
-                UI::get ( high, set, "high", UI::compulsory );
-                UI::get ( low, set, "low", UI::compulsory );
-                set_parameters ( low, high );
-                validate();
-        }
-        void validate() {
-                bdm_assert ( high.length() == low.length(), "Incompatible high and low vectors" );
-                dim = high.length();
-                bdm_assert ( min ( distance ) > 0.0, "bad support" );
-        }
+        void from_setting ( const Setting &set );
+        void    to_setting  (Setting  &set) const;
+        void validate();
 };
 UIREGISTER ( euni );
@@ -922 +829 @@
                 iepdf.set_parameters ( -delta, delta );
         }
-        
+        void    to_setting  (Setting  &set) const {
+                pdf::to_setting ( set ); 
+                UI::save( iepdf.mean(), set, "delta");
+                UI::save(mean, set, "mean");
+        }
         void validate(){
+                pdf_internal<euni>::validate();
                 dimc = mean->dimensionc();
-                iepdf.validate();
+                
         }
 
@@ -1006 +918 @@
                 UI::get ( R0, set, "R", UI::compulsory );
                 set_parameters ( A, mu_const, R0 );
-                validate();
-        };
-        void to_setting (Setting &set) const {
+        }
+
+void to_setting (Setting &set) const {
                 pdf::to_setting(set);
                 UI::save ( A, set, "A");
@@ -1014 +926 @@
                 UI::save ( _R(), set, "R");
         }
-        void validate() {
+
+void validate() {
                 pdf_internal<TEpdf<sq_T> >::validate();
                 if (mu_const.length()==0) { // default in from_setting
@@ -1070 +983 @@
         */
 
-        void from_setting ( const Setting &set ) {
+
+void from_setting ( const Setting &set ) {
                 pdf::from_setting ( set );
                 shared_ptr<fnc> g = UI::build<fnc> ( set, "g", UI::compulsory );
@@ -1084 +998 @@
                 //validate();
         }
-        void validate() {
+        
+
+void to_setting  (Setting  &set) const {
+                UI::save( g,set, "g");
+                UI::save(iepdf._R().to_mat(),set, "R");
+        
+        }
+
+
+
+void validate() {
                 this->iepdf.validate();
                 bdm_assert ( g->dimension() == this->iepdf.dimension(), "incompatible function" );
@@ -1091 +1015 @@
                 this->iepdf.validate();
         }
+
 };
 
@@ -1128 +1053 @@
 
         void validate() {
+                mlnorm<ldmat, enorm>::validate();
                 bdm_assert ( A.rows() == mu_const.length(), "mlstudent: A vs. mu mismatch" );
                 bdm_assert ( _R.rows() == A.rows(), "mlstudent: A vs. R mismatch" );
@@ -1175 +1101 @@
          \endcode
         */
-        void from_setting ( const Setting &set ) {
-                pdf::from_setting ( set ); // reads rv and rvc
-                vec betatmp; // ugly but necessary
-                UI::get ( betatmp, set, "beta", UI::compulsory );
-                UI::get ( k, set, "k", UI::compulsory );
-                set_parameters ( k, betatmp );
-                validate();
-        }
-        void validate() {
-                pdf_internal<egamma>::validate();
-
-                dim = _beta.length();
-                dimc = _beta.length();
-        }
+        void from_setting ( const Setting &set );
+        void    to_setting  (Setting  &set) const;
+        void validate();
 };
 UIREGISTER ( mgamma );
@@ -1236 +1151 @@
 
         void    validate (){
-                iepdf.validate();
+                pdf_internal<eigamma>::validate();
                 dimc = dimension();
 };
@@ -1338 +1253 @@
          */
         void from_setting ( const Setting &set );
-
+        void to_setting  (Setting  &set) const; 
         // TODO dodelat void to_setting( Setting &set ) const;
 };
@@ -1395 +1310 @@
         
         void validate(){
-                iepdf.validate();
+                pdf_internal<elognorm>::validate();
                 dimc = iepdf.dimension();
         }
@@ -1415 +1330 @@
         */
         void from_setting ( const Setting &set );
+        void to_setting  (Setting  &set) const; 
 
         // TODO dodelat void to_setting( Setting &set ) const;
@@ -1449 +1365 @@
         
         virtual void validate (){
+                epdf::validate();
                 dim = p * p;
         }
@@ -1519 +1436 @@
         
         virtual void    validate (){
+                epdf::validate();
                 W.validate();
                 dim = W.dimension();
@@ -1594 +1512 @@
         
         void validate(){
-                iepdf.validate();
+                pdf_internal<eiWishartCh>::validate();
                 dimc = iepdf.dimension();
         }
@@ -1657 +1575 @@
                 samples = Av;
         };
-        virtual void    validate (){
-          bdm_assert (samples.length()==w.length(),"samples and weigths are of different lengths");
-          n = w.length();
-          if (n>0)
-                pdf::dim = samples ( 0 ).length();
-        }
+        virtual void    validate ();
         //! Potentially dangerous, use with care.
         vec& _w()  {
@@ -1709 +1622 @@
         void qbounds ( vec &lb, vec &ub, double perc = 0.95 ) const;
 
-        void to_setting ( Setting &set ) const {
-                epdf::to_setting( set );
-                UI::save ( samples, set, "samples" );
-                UI::save ( w, set, "w" );
-        }
-
-        void from_setting ( const Setting &set ) {
-                epdf::from_setting( set );
-                
-                UI::get( samples, set, "samples", UI::compulsory );
-                UI::get ( w, set, "w", UI::compulsory );
-                validate();
-        }
+        void to_setting ( Setting &set ) const;
+        void from_setting ( const Setting &set );
 
 };
@@ -1745 +1647 @@
         UI::get ( Rtmp, set, "R", UI::compulsory );
         R = Rtmp; // conversion
-        validate();
 }
+
+template<class sq_T>
+void enorm<sq_T>::validate() {
+                eEF::validate();
+                bdm_assert ( mu.length() == R.rows(), "mu and R parameters do not match" );
+                dim = mu.length();
+        }
 
 template<class sq_T>
@@ -1754 +1662 @@
         UI::save ( R.to_mat(), set, "R");
 }
+
+
 
 template<class sq_T>

library/bdm/stat/merger.cpp

@@ -9 +9 @@
         set_sources ( S );
 }
+
 
 void merger_base::set_sources ( const Array<shared_ptr<pdf> > &Sources ) {
@@ -348 +349 @@
 }
 
+void merger_mix::from_setting ( const Setting& set ) {
+                merger_base::from_setting ( set );
+                Ncoms=DFLT_Ncoms;
+                UI::get( Ncoms, set, "ncoms", UI::optional );
+                effss_coef=DFLT_effss_coef;
+                UI::get (effss_coef , set,  "effss_coef", UI::optional);
+                stop_niter=10;
+                UI::get ( stop_niter, set,"stop_niter", UI::optional );
+        
+        }
+
+void    merger_mix::to_setting  (Setting  &set) const  {
+                merger_base::to_setting(set);
+                UI::save( Ncoms, set, "ncoms");
+                UI::save (effss_coef , set,  "effss_coef");
+                UI::save ( stop_niter, set,"stop_niter");
+
+        }
+
+void merger_mix::validate() {
+                merger_base::validate();
+                bdm_assert(Ncoms>0,"Ncoms too small");
+        }
+
+void merger_base::from_setting ( const Setting& set ) {
+                // get support
+                // find which method to use
+                epdf::from_setting (set);
+                string meth_str;
+                UI::get<string> ( meth_str, set, "method", UI::compulsory );
+                if ( !strcmp ( meth_str.c_str(), "arithmetic" ) )
+                        set_method ( ARITHMETIC );
+                else {
+                        if ( !strcmp ( meth_str.c_str(), "geometric" ) )
+                                set_method ( GEOMETRIC );
+                        else if ( !strcmp ( meth_str.c_str(), "lognormal" ) ) {
+                                set_method ( LOGNORMAL );
+                                set.lookupValue ( "beta", beta );
+                        }
+                }
+                string dbg_file;
+                if ( UI::get ( dbg_file, set, "dbg_file" ) )
+                        set_debug_file ( dbg_file );
+                //validate() - not used
+        }
+        
+        void merger_base::to_setting  (Setting  &set) const {
+                epdf::to_setting(set);
+                
+                UI::save( METHOD, set, "method");
+                UI::save ( dbg_file->get_fname(), set, "dbg_file" );
+                UI::save (beta, set, "beta" );
+        }
+
+        void merger_base::validate() {
+//              bdm_assert ( eSmp._w().length() > 0, "Empty support, use set_support()." );
+//              bdm_assert ( dim == eSmp._samples() ( 0 ).length(), "Support points and rv are not compatible!" );
+                epdf::validate();
+                bdm_assert ( isnamed(), "mergers must be named" );
+        }
 // DEFAULTS FOR MERGER_BASE
 const MERGER_METHOD merger_base::DFLT_METHOD = LOGNORMAL;

library/bdm/stat/merger.h

@@ -166 +166 @@
 
         //! load from setting
-        void from_setting ( const Setting& set ) {
-                // get support
-                // find which method to use
-                string meth_str;
-                UI::get<string> ( meth_str, set, "method", UI::compulsory );
-                if ( !strcmp ( meth_str.c_str(), "arithmetic" ) )
-                        set_method ( ARITHMETIC );
-                else {
-                        if ( !strcmp ( meth_str.c_str(), "geometric" ) )
-                                set_method ( GEOMETRIC );
-                        else if ( !strcmp ( meth_str.c_str(), "lognormal" ) ) {
-                                set_method ( LOGNORMAL );
-                                set.lookupValue ( "beta", beta );
-                        }
-                }
-                string dbg_file;
-                if ( UI::get ( dbg_file, set, "dbg_file" ) )
-                        set_debug_file ( dbg_file );
-                //validate() - not used
-        }
-
-        void validate() {
-//              bdm_assert ( eSmp._w().length() > 0, "Empty support, use set_support()." );
-//              bdm_assert ( dim == eSmp._samples() ( 0 ).length(), "Support points and rv are not compatible!" );
-                bdm_assert ( isnamed(), "mergers must be named" );
-        }
+        void from_setting ( const Setting& set );
+        
+        void to_setting  (Setting  &set) const ;
+
+        void validate() ;
         //!@}
 };
@@ -268 +247 @@
         }
         //! from_settings
-        void from_setting ( const Setting& set ) {
-                merger_base::from_setting ( set );
-                Ncoms=DFLT_Ncoms;
-                UI::get( Ncoms, set, "ncoms", UI::optional );
-                effss_coef=DFLT_effss_coef;
-                UI::get (effss_coef , set,  "effss_coef", UI::optional);
-                stop_niter=10;
-                UI::get ( stop_niter, set,"stop_niter", UI::optional );
-                
-                validate();
-        }
-        void validate() {
-                bdm_assert(Ncoms>0,"Ncoms too small");
-        }
+        void from_setting ( const Setting& set );
+        void    to_setting  (Setting  &set) const;
+        void validate();
 
         //! @}