OpenPaperOpt/Staticsheet.h

Go to the documentation of this file.

#ifndef STATICSHEET_H
#define STATICSHEET_H

#include "utilities.h"
#include "StructureObject.h"
#include "Array2D.h"
#include "V2.h"
#include "SurfAcc.h"
#include "FiberStatic.h"
#include "Distribution.h"
#include "DistributionVonMises.h"
#include "DistributionComponents.h"
#include "DistributionOneDimension.h"
#include "DistributionThreeDimension.h"
#include "DistributionFromFiberSamples.h"

namespace p3d{
        
        class Fiber;
        class Staticsheet: public StructureObject{
        
        public:
                Staticsheet(double xsize,double ysize,double seglen);
                ~Staticsheet();

                int GetDataSize();
                void SaveStaticsheet(string filename);
                int ReadFiberFromFile(string filename);

                int     Process(WavePacket & wp, int up_down,int internal) const {return 0;};
                
                int AddFibers(int num,double length,double radius,double formFC);
                int AddFibers(Distributions::Distribution *pDistlen, Distributions::Distribution *pDistwid, Distributions::Distribution * formFC, int num);
                
                int AddFibers(Distributions::DistributionOneDimension *pDistlen,Distributions::DistributionThreeDimension * pDsitWidFactor, int num);
                
                int AddFibers(Distributions::DistributionFromFiberSamples *pDistSource, int num);
                
                void ShuffleFibers();

                void SetRandomOrientation_vonMises(double kapa,double theta);

                void ApproximateLocalCurl();

                void PlaceFibers(bool killAngularInformation=true,int numPlaceTests=5);

                void MakePaper3DStructure(Array2D<float> * surf); 

                void CalcHeightData(double* mat,int numCells);


                //Accessors / estimators
                double  GetSizeX() const {return mSizeX;} //um
                double  GetSizeY() const {return mSizeY;} //um
                double  GetTotalWeight()   const {return mTotalWeight;}
                //add new set function
                void    SetTotalWeight(double totalWeight) {mTotalWeight=totalWeight;}
                double  GetSurfaceWeight() const {return mTotalWeight/(mSizeX*mSizeY)*1e12;} //gram/m^2, verry bad approximation!!!!
                void   EstimateLocalSurfWeight(int res,double* data); //Res must be on form 2^n. Assumes little bending of each segment
                void   EstimateSurfaces(double probesize,int* resp,float** datau,float** datal,double* meandiff); 
                                                      //probesize in um. Selects res on form 2^n
                                                      //alocates matrixes, if set, user responcibility to clear
                void SaveEstimateSurfaces(string filename,int res,float* datau,float* datal,double meandiff);
                int    GetSegType() const {return mSegType;}
                int    GetFlags()   const {return mInitFlg;}

                void    SetSegType(int segType) {mSegType=segType;}
                void    SetFlags(int initFlg) {mInitFlg=initFlg;}

                double GetTop()     const {if (!(mInitFlg&SIMVOL_FLAG_POS3D)) throw StdExcept("Staticsheet:: :No 3d positions");return mSStr_top;}
                double GetMedTop()  const {if (!(mInitFlg&SIMVOL_FLAG_POS3D)) throw StdExcept("Staticsheet:: :No 3d positions");return mSStr_mtop;}
                double GetBot()     const {if (!(mInitFlg&SIMVOL_FLAG_POS3D)) throw StdExcept("Staticsheet:: :No 3d positions");return mSStr_bot;}
                double GetMedBot()  const {if (!(mInitFlg&SIMVOL_FLAG_POS3D)) throw StdExcept("Staticsheet:: :No 3d positions");return mSStr_mbot;}

                void SetTop(double SStr_top)    {mSStr_top=SStr_top;}
                void SetMedTop(double SStr_mtop)        {mSStr_mtop=SStr_mtop;}
                void SetBot(double SStr_bot)    {mSStr_bot=SStr_bot;}
                void SetMedBot(double SStr_mbot)        {mSStr_mbot=SStr_mbot;}

                vector<Fiber*>::const_iterator FiberBegin() const {return mFibers.begin();}
                vector<Fiber*>::const_iterator FiberEnd()   const {return mFibers.end();}
        
        private:

                //Size of surface
                double mSizeX;   //um
                double mSizeY;   //um

                //Surface structure (generated by 3d-generator, avaliable when SIMVOL_FLAG_POS3D set)
                double mSStr_top;
                double mSStr_mtop;
                double mSStr_bot;
                double mSStr_mbot;

                //Segment length for fibers
                double mSegLen;  //um

                //Data validity of fibers
                int    mInitFlg; //SIMVOL_FLAG_xxx
                int    mSegType; //SEGT_xxx

                //Orgin
                string mOrgin;

                vector<Fiber*> mFibers;
                double         mTotalWeight;
        };

}

#endif