ESyS-Particle  4.0.1
RandomBlock.h
00001 
00002 //                                                         //
00003 // Copyright (c) 2003-2011 by The University of Queensland //
00004 // Earth Systems Science Computational Centre (ESSCC)      //
00005 // http://www.uq.edu.au/esscc                              //
00006 //                                                         //
00007 // Primary Business: Brisbane, Queensland, Australia       //
00008 // Licensed under the Open Software License version 3.0    //
00009 // http://www.opensource.org/licenses/osl-3.0.php          //
00010 //                                                         //
00012 
00013 #ifndef __RANDOMBLOCK_H
00014 #define __RANDOMBLOCK_H
00015 
00016 //-- project includes --
00017 #include "Geometry/SimpleParticle.h"
00018 #include "Geometry/SimpleNTable.h"
00019 #include "Geometry/Sphere2d.h"
00020 #include "Geometry/Line.h"
00021 #include "Geometry/RandomAssembly2D.h"
00022 
00023 //-- STL includes --
00024 #include <vector>
00025 #include <string>
00026 #include <utility>
00027 
00028 using std::vector; 
00029 using std::string;
00030 using std::pair;
00031 
00040 class CRandomBlock2D : public ARandomAssembly2D
00041 {
00042  protected:
00043   virtual Vec3 getAPoint();
00044   virtual int getNParts() const{return m_bpart.size();};
00045   double m_maxConnDist;
00046 
00047  public:
00048   CRandomBlock2D(double,double,double,double,double,double,double,bool circ_x=false);
00049   virtual ~CRandomBlock2D();
00050 
00051   virtual void generate(int,unsigned int);
00052   virtual void insertParticle(const SimpleParticle);
00053   virtual void tagParticleClosestTo(const Vec3&,int); 
00054   virtual void tagEdgeY(int,int,double);
00055   virtual void tagEdgeZ(int,int,double){}; // do nothing
00056  
00057   virtual void writeToGeoFile(const string&);
00058   virtual double calcPorosity();
00059   virtual vector<pair<double,double> > getSizeDistribution(int);
00060 };
00061 
00062 #endif // __RANDOMBLOCK_H