ESyS-Particle  4.0.1
RandomBlock3D.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 __RANDOMBLOCK3D_H
00014 #define __RANDOMBLOCK3D_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/RandomAssembly3D.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 CRandomBlock3D : public ARandomAssembly3D
00041 {
00042  private:
00043   bool m_is_bonded;
00044 
00045  protected:
00046   virtual Vec3 getAPoint();
00047   virtual int getNParts() const{return m_bpart.size();};
00048   double m_maxConnDist;
00049 
00050  public:
00051   CRandomBlock3D(double,double,double,double,double,double,double,double,double,bool,bool is_bonded=true);
00052   virtual ~CRandomBlock3D();
00053 
00054   virtual void generate(int,unsigned int);
00055   virtual void insertParticle(const SimpleParticle);
00056   virtual void tagParticleClosestTo(const Vec3&,int); 
00057   virtual void tagEdgeY(int,int,double);
00058   virtual void tagEdgeZ(int,int,double);
00059  
00060   virtual void writeToGeoFile(const string&);
00061   virtual double calcPorosity();
00062   virtual vector<pair<double,double> > getSizeDistribution(int);
00063 };
00064 
00065 #endif //  __RANDOMBLOCK3D_H
00066