Changeset 569 for library/bdm/stat/discrete.h

Timestamp:

08/20/09 00:54:47 (15 years ago)

Author:

smidl

Message:

new object discrete_support, merger adapted to accept this input as well

Files:

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

library/bdm/stat/discrete.h

@@ -23 +23 @@
 
         //! Rectangular support
-        //! Support ponits are located inbetween ranges! For example:
+        //! Support points are located inbetween ranges! For example:
         //! For ranges=[0,1] and gridsizes=[1] the support point is 0.5
         class rectangular_support: public root {
@@ -89 +89 @@
                                 return ind;
                         } 
-                        //! set the first corner to actvec
+                        //! set the first vector to corner and store result in actvec
                         const vec& first_vec(){
                                 for ( int j = 0; j < dim; j++ ) {
@@ -146 +146 @@
         UIREGISTER(rectangular_support);
         
+        //! Discrete support with stored support points 
+        class discrete_support: public root{
+                protected:
+                        //! storage of support points
+                        Array<vec> Spoints;
+                        //! index in iterators
+                        int idx;
+                public:
+                        //! Default constructor
+                        discrete_support() : Spoints(0), idx(0){}
+                        //! Access function
+                        int points() const {return Spoints.length();}
+                        //! set the first vector to corner and store result in actvec
+                        const vec& first_vec(){bdm_assert_debug(Spoints.length()>0,"Empty support");idx=0; return Spoints(idx);}
+                        //! set next vector after calling first_vec()
+                        const vec& next_vec(){bdm_assert_debug(Spoints.length()>idx-1,"Out of support points"); return Spoints(++idx);}
+                        
+                        /*!
+                        \code
+                          class = "discrete_support";
+                          points = ( [1,2..], [2,2..], ...); // list of points
+                           === OR ===
+                          epdf = {class="epdf_offspring",...}; // epdf rfom which to sample
+                          npoints = 100;                     // number of samples
+                        \endcode
+                        */
+                        void from_setting (const Setting &set) {
+                                UI::get (Spoints, set, "points", UI::compulsory);
+                                if (points()<1){
+                                        int npoints;
+                                        shared_ptr<epdf> ep= UI::build<epdf>(set, "epdf", UI::compulsory);
+                                        if (!UI::get(npoints,set,"npoints",UI::optional)){npoints=100;}
+                                        
+                                        //sample
+                                        Spoints.set_size(npoints);
+                                        for(int i=0; i<points(); i++){Spoints(i)=ep->sample();}
+                                }
+                        }
+                        //! access function
+                        Array<vec> & _Spoints() {return Spoints;}
+        };
+        UIREGISTER(discrete_support);
+        
         class grid_fnc: public fnc{
                 protected: