OpenPaperOpt/SurfAcc.h

Go to the documentation of this file.

#ifndef SURFACC_H
#define SURFACC_H

#include "utilities.h"
#include "V2.h"
#include <fstream>
#include "Array2D.h"

namespace p3d{
        class SurfAcc {
        public:

                SurfAcc(){}
                ~SurfAcc(){}
                SurfAcc(double sizeX,double sizeY,int nSegX,int nSegY); //Nseg must be on form 2^n
                void ResetAccs(double value=0); //Can only be used with height operations!!
                void SetAccs(const SurfAcc* source);
                void SetAccs(const Array2D <float> *source);

                //==========================Area based operations===========================
                void CalcPollySupport(double totVal,V2 <double>* p1,V2 <double>* p2,V2 <double>* p3,V2 <double>* p4,double linSmpPerSq=5);

                void   AddSupportToAcc();
                double Calc3DegIntegral(double c1,double c2,double c3) const;
                inline double CalcPiecLinIntegral(double g1,double g2,double k) const;

                //========================Height based operations=======================
                void CalcPollyHeights(double h2,double h4,V2 <double>* p1,V2 <double>* p2,V2 <double>* p3,V2 <double>* p4,int useMin=1,double linSmpPerSq=5);
                double GetMinHeightSpan();


                void AddHeightsToAcc(int keepMax=1);

                void SaveSurfAccBinary(string filename);
                void ReadSurfAccFromBinary(string filename);
                //Sets base to max(support,base)
                //Accessors - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

                const double* GetAccumulator() {return mAccs;}
                double* GetAccumulatorModify() {return mAccs;} //You'd better know what you'r doing :)
                double  GetBoxSizeX() {return mSizeX/mNSegX;}
                double  GetBoxSizeY() {return mSizeY/mNSegY;}
                double  GetTotal() const {return mTotal;}
                double  GetMean() const;

        private:
                double  mSizeX,mSizeY;    //um, size of entire area
                double  mSWidX,mSWidY;    //um, size of segments
                double  mSWidXI,mSWidYI;  //Inverse of above
                int     mNSegX,mNSegY;    //Number of segments
                double  mSmalestSegInv;   //1/min(mSegX,mSegY);
                double  mTotal;           //Total accumulated value / (max/min) height

                unsigned int  mBinMaskX,mBinMaskY,mBinShiftY; //(x,y)->Idx


                AutoDelete<double>      mAccs;    //Accumulator array
                AutoDelete<double>      mSupport; //Current suport

                int                    mSupXOfs,mSupYOfs; //x,y-ofsets in array
                int                    mSupXSize,mSupYSize; //Number of elements in each direction (safe included)
                int                    mSupCapacity;
                int                    mSupValid; //Validity flag
                double                 mSupTotal; //Total volume
                double                 mSupScaler;
        };

}
#endif;