#include <InsertionBoundary.h>

Inheritance diagram for InsertionBoundary:

Public Member Functions
virtual InsertionBoundary *	copy () const
	BaseBoundary copy method. More...

void	set (Vec3D posMin, Vec3D posMax, Vec3D velMin, Vec3D velMax, Mdouble radMin, Mdouble radMax, int maxFailed)

virtual void	checkBoundaryActionsBeforeTimeStep (ParticleHandler &pH, WallHandler &wH, RNG &random)

void	read (std::istream &is)
	reads wall More...

void	print (std::ostream &os) const
	outputs wall More...

Public Member Functions inherited from BaseBoundary
virtual	~BaseBoundary ()

virtual int	createPeriodicParticles (BaseParticle *P UNUSED, ParticleHandler &pH UNUSED)

virtual bool	checkBoundaryAfterParticleMoved (BaseParticle *P UNUSED, ParticleHandler &pH UNUSED)

virtual void	checkBoundaryActionsBeforeTimeStep (ParticleHandler &pH UNUSED, WallHandler &wH UNUSED, RNG &random UNUSED)

void	set_Index (int index)

void	set_Id (int id)

void	setHandler (BoundaryHandler *handler)

virtual void	read (std::istream &is UNUSED)=0
	reads boundary More...

virtual void	print (std::ostream &os UNUSED) const =0
	outputs boundary More...

virtual void	moveInHandler (int newPos)

Public Attributes
Vec3D	posMin_

Vec3D	posMax_

Vec3D	velMin_

Vec3D	velMax_

Mdouble	radMin_

Mdouble	radMax_

int	maxFailed_

Detailed Description

Definition at line 31 of file InsertionBoundary.h.

Member Function Documentation

virtual void InsertionBoundary::checkBoundaryActionsBeforeTimeStep	(	ParticleHandler &	pH,
		WallHandler &	wH,
		RNG &	random
	)

inlinevirtual

Definition at line 53 of file InsertionBoundary.h.

References WallHandler::begin(), BaseHandler< T >::begin(), BaseHandler< T >::copyAndAddObject(), WallHandler::end(), BaseHandler< T >::end(), BaseParticle::get_Position(), BaseParticle::get_Radius(), RNG::get_RN(), maxFailed_, posMax_, posMin_, radMax_, radMin_, BaseParticle::set_inertia(), BaseParticle::set_Mass(), BaseParticle::set_Position(), BaseParticle::set_Radius(), BaseParticle::set_Velocity(), sqr, velMax_, velMin_, Vec3D::X, Vec3D::Y, and Vec3D::Z.

     {
         int failed = 0;
         BaseParticle p0;
         bool isInsertable;
         Mdouble distance;
         Vec3D normal;
         //try max_failed times to find new insertable particle
         while (failed<=maxFailed_)
         {
             std::cout << "I" << std::endl;
             //Set random radius, position and velocity
             p0.set_Radius(random.get_RN(radMin_,radMax_));
             p0.set_Mass(1.0);                              
             p0.set_inertia(1.0);
             Vec3D position;
             position.X = random.get_RN(posMin_.X,posMax_.X);
             position.Y = random.get_RN(posMin_.Y,posMax_.Y);
             position.Z = random.get_RN(posMin_.Z,posMax_.Z);
             Vec3D velocity;
             velocity.X = random.get_RN(velMin_.X,velMax_.X);
             velocity.Y = random.get_RN(velMin_.Y,velMax_.Y);
             velocity.Z = random.get_RN(velMin_.Z,velMax_.Z);
             p0.set_Position(position);
             p0.set_Velocity(velocity);
             // p0.set_Position(Vec3D(random.get_RN(posMin_.X,posMax_.X),random.get_RN(posMin_.Y,posMax_.Y),random.get_RN(posMin_.Z,posMax_.Z)));
             // p0.set_Velocity(Vec3D(random.get_RN(velMin_.X,velMax_.X),random.get_RN(velMin_.Y,velMax_.Y),random.get_RN(velMin_.Z,velMax_.Z)));
             
             isInsertable=true;
             
 //            cout<<"Particle to insert "<<p0<<endl;
             //Check if it has no collision with walls
             for(std::vector<BaseWall*>::iterator it=wH.begin();isInsertable&&it!=wH.end();it++)
             {
                 (*it)->get_distance_and_normal(p0,distance,normal);
                 if(distance<p0.get_Radius())
                 {
  //                   cout<<"Collision with wall "<<**it<<endl;
                     isInsertable=false;
                 }
                 else
                 {
 //                    cout<<"No collision with wall "<<**it<<endl;
                 }
             }
             //Check if it has no collision with other particles
             for(std::vector<BaseParticle*>::iterator it=pH.begin();isInsertable&&it!=pH.end();it++)
             {
                 if(GetDistance2(p0.get_Position(),(*it)->get_Position())<sqr(p0.get_Radius()+(*it)->get_Radius()))
                 {
  //                   cout<<"Collision with particle "<<**it<<endl;
                     isInsertable=false;
                 }
                 else
                 {
 //                    cout<<"No collision with particle "<<**it<<endl;
                 }
             }
             
             if (isInsertable)
             {
                 //Still needs to set a lot of things like mass and HGRID stuff
  //               cout<<"Inserting particle"<<endl;
                 pH.copyAndAddObject(p0);
                 failed = 0; 
             }
             else
             {
 //                cout<<"Not inserting particle"<<endl;
                 failed++;
             }
         }      
     };

virtual InsertionBoundary* InsertionBoundary::copy ( ) const

inlinevirtual

BaseBoundary copy method.

It calls the copy constructor of this BaseBoundary, useful for polymorphism todo{Does this work correctly?}

Implements BaseBoundary.

Definition at line 34 of file InsertionBoundary.h.

     {
         #ifdef CONSTUCTOR_OUTPUT
             cerr << "virtual InsertionBoundary* copy() const finished" << endl;
         #endif      
         return new InsertionBoundary(*this);
     }

void InsertionBoundary::print ( std::ostream & os ) const

inline

outputs wall

Definition at line 134 of file InsertionBoundary.h.

References maxFailed_, posMax_, posMin_, radMax_, radMin_, velMax_, and velMin_.

                                      {
         os << "InsertionBoundary posMin " << posMin_
             << " posMax " << posMax_ 
             << " velMin " << velMin_
             << " velMax " << velMax_
             << " radMin " << radMin_
             << " radMax " << radMax_
             << " maxFailed " << maxFailed_;
     }

void InsertionBoundary::read ( std::istream & is )

inline

reads wall

Definition at line 128 of file InsertionBoundary.h.

References maxFailed_, posMax_, posMin_, radMax_, radMin_, velMax_, and velMin_.

                                { 
         std::string dummy;
         is >> dummy >> posMin_ >> dummy >> posMax_ >> dummy >> velMin_ >> dummy >> velMax_ >> dummy >> radMin_ >> dummy >> radMax_ >> dummy >> maxFailed_;
     }

void InsertionBoundary::set	(	Vec3D	posMin,
		Vec3D	posMax,
		Vec3D	velMin,
		Vec3D	velMax,
		Mdouble	radMin,
		Mdouble	radMax,
		int	maxFailed
	)

inline

Definition at line 42 of file InsertionBoundary.h.

References maxFailed_, posMax_, posMin_, radMax_, radMin_, velMax_, and velMin_.

     {
         posMin_=posMin;
         posMax_=posMax;
         velMin_=velMin;
         velMax_=velMax;
         radMin_=radMin;
         radMax_=radMax;
         maxFailed_=maxFailed;        
     }