Event.h

Go to the documentation of this file.

#ifndef  LS2_Event_H
#define  LS2_Event_H

#include "tarch/la/Vector.h"
#include "GlobalValue.h"
namespace algorithms{
namespace ls2{
class Event;
}
}
#define INF    100000000
#define dblINF 100000000.

class algorithms::ls2::Event {

public:
        Event(double time_i, int i_i, int j_i, tarch::la::Vector<DIM,int> v_i);
        Event(double time_i, int i_i, int j_i);
        Event(const Event& e);
        Event();
        ~Event();
        //friend class Box;
        bool operator<(const Event&) const;
        bool operator>(const Event&) const;
        void erase();
        //by Tan
        double getTime();
        void setTime(double time);
        int getCollisionPartnerI();
        void setCollisionPartnerI(int i);
        int getCollisionPartnerJ();
        void setCollisionPartnerJ(int j);
        tarch::la::Vector<DIM,int> getVector();
        void setVector(tarch::la::Vector<DIM,int> v);

private:

        double time;
        int i;
        int j;
        tarch::la::Vector<DIM,int> v;

        /* 0<=j<=N                binary collision between i and j
     j=N+DIM+1+x            transfer where x=-(k+1) for left wall
                            and x=k+1 for right wall
     j=INF                 both check after Event that did not altered motion of                           i and check after Event that altered motion of i, i.e                           rescaling of velocities. I currently don't see need t                           o separate the two

     j=-1                  check after collision

     Virtual identifiers as scalars...I think bad idea, but here's my work
     there will be easily fixed problems if DIM>=10
     -x<=v<=x               x=k+1, image in k direction
     v=xy                   x,y
      =-xy                  -x,y
      =-yx                  x,-y
      =yx                   -x,-y
     v=xyz                  x,y,z
      =-xyz                 -x,y,z
      =-yxz                 x,-y,z
      =-zxy                 x,y,-z
      =zyx                  -x,-y,z
      =xzy                  x,-y,-z
      =yzx                  -x,y,-z
      =-zyx                 -x,-y,-z
         */

};

#endif