OpenPaperOpt/DistributionThreeDimension.h

Go to the documentation of this file.

#ifndef DISTRIBUTIONTHREEDIMENSION_H
#define DISTRIBUTIONTHREEDIMENSION_H

#include "utilities.h"
#include "Distribution.h"
#include "DistributionThreeDimensionBSP.h"
namespace Distributions
{
        class DistributionThreeDimension:public Distribution {
        public:
                DistributionThreeDimension();
                ~DistributionThreeDimension();
                
                DistributionThreeDimension(double minx,double maxx,double miny,double maxy,
                                                                                double maxsidex,double maxsidey,int minpoints,double maxSkew,
                                                                                int minPntsZ,double* data,int numData,int dataDist); 

                int  GetNum3dPatches(int FillBorders) {return mNumLief*(FillBorders?5:1);}
                void Buid3dPatchSet(int FillBorders,double* x,double* y,double* z,int count);
                
                //int GetDataSize();

                double  GenerateRandomNumber(double len ,double *rnd);

                
        
        protected:
                void   BuildScalerTree();
                double GetScaler(double y);

        private:

                double InvF(double X,double y,double* z=0); 
                double rec_InvF(double minx,double maxx,double miny,double maxy,double Xleft,
                                                double y,DistributionThreeDimensionBSP* tree,double* value,ZDist* bestzd,ZDist** zd);

                double rec_GetScaler(double miny,double maxy,double targety,void** tre);

                void rec_BuildScalerTree(DistributionThreeDimensionBSP* src,TreeNode* tre,double minx,double maxx,int* splitCnt);
                void rec_IntScalerTree(TreeNode* tre,double sum,void*** ary_tre,double** ary_dta);

                DistributionThreeDimensionBSP* RecBuild2D(double minx,double maxx,double miny,double maxy,
                                                                                                        fsamplesiter begin,fsamplesiter end,double maxsidex,
                                                                                                  double maxsidey,int minpoints,double maxSkew,int* liefCntr,int minPntsZ);

                double    mMinX,mMaxX;
                double    mMinY,mMaxY;
                int      mNumLief; //Number of lief nodes in hierarchy
                DistributionThreeDimensionBSP* mTree;    
                
                //Storage of 1/weights(y)
                AutoDelete<void*>  mYWTree; //Tree-structure, lobit set->ptr to mYWData
                AutoDelete<double> mYWData; //Data
                int                mYWSize; //Total datasize allocated for structure
        };
}
#endif