geometry.h

Go to the documentation of this file.

#ifndef GEOMETRY_H
#define GEOMETRY_H
#include "task.h"
#include "numeric.h"
#define ID_GEOMETRY "Geometry"

//---------------------------------------------------------------------------------------------------------------------
class Geometry : public Task {
 protected:
  GridPosition GridSize;
  Position GridSpacing;
 public:
  double Rd;
  int number_d;
  Geometry () : Task() {};
  virtual ~Geometry () {};
  virtual std::string GetClassName () const {return ID_GEOMETRY;}
  virtual void Init (); 
  virtual bool Execute () {return true;}; 
  GridPosition GetGridSize () const {return GridSize;}
  Position GetGridSpacing () const {return GridSpacing;}
  virtual Position MinPosition() const = 0;
  virtual Position MaxPosition() const = 0;


  virtual double GetDustRadius();
#ifndef ONE_DIMENSIONAL

  virtual int GetNumberDust();
  virtual GridPosition GetDustPosition(int);
#endif

  virtual Position GetLocation (const GridPosition&) const = 0;
  virtual double Distance(const GridPosition&, const GridPosition&) const = 0;
  double GetVolumeQuant () const {return GridSpacing.volume(); }
  virtual double GetBoxVolume () const = 0;
  
#ifdef ONE_DIMENSIONAL
  virtual bool outside(const Position&) const = 0 ;
  

  virtual bool insidedust(const Position&) const = 0; 
  virtual bool outsideRight(const Position&) const = 0;
#endif
#ifndef ONE_DIMENSIONAL
  int ig;
  Position pt,pd;
#ifdef TWO_DIMENSIONAL

  Position pc[2][2000];
#endif
#ifdef THREE_DIMENSIONAL

  Position pc[3][2000];
#endif

  GridPosition Dust[20];
  virtual bool outside(const Position&) const = 0;
  virtual bool outside_periodic(const Position&) const = 0;
  virtual bool outside_right(const Position&) const = 0; 
  virtual bool outside_left(const Position&) const = 0;
  virtual bool insidedust(const Position&, GridPosition) const = 0; 
#endif
  

  virtual void DivByVolume(grid<double> &G) = 0;
#ifndef THREE_DIMENSIONAL
  virtual void DivByVolume(grid<Position> &G) = 0;
  virtual void DivByVolume(grid<Velocity> &G) = 0;
#endif


protected:
  virtual PARAMETERMAP* MakeParamMap (PARAMETERMAP* pm = NULL);  
};

class ClosedGeometry : public Geometry {
 public:
  double Rd;
#ifndef ONE_DIMENSIONAL
  int x,y;
  int Rg;
  double dx2;
  double dy2;
  double Rq;
  double pos1x,pos1y,pos1z,pos2x,pos2y,pos2z,pos3x,pos3y,pos3z,pos4x,pos4y,pos4z,pos5x,pos5y,pos5z;
  double pos6x,pos6y,pos6z,pos7x,pos7y,pos7z,pos8x,pos8y,pos8z,pos9x,pos9y,pos9z,pos10x,pos10y,pos10z;
  double pos11x,pos11y,pos11z,pos12x,pos12y,pos12z,pos13x,pos13y,pos13z,pos14x,pos14y,pos14z,pos15x,pos15y,pos15z; 
  double pos16x,pos16y,pos16z,pos17x,pos17y,pos17z,pos18x,pos18y,pos18z,pos19x,pos19y,pos19z,pos20x,pos20y,pos20z;
#ifdef TWO_DIMENSIONAL
  int xp,yp;
#endif
#ifdef THREE_DIMENSIONAL
  int z;
  double dz2;
#endif
#endif

  int number_d;
  ClosedGeometry () : Geometry() {};
  virtual ~ClosedGeometry () {};
  virtual void Init();

  virtual Position MinPosition() const {return Position(0.);}
  virtual Position MaxPosition() const {return Position(GridSize);};


  virtual Position GetLocation (const GridPosition&) const;
  virtual double Distance(const GridPosition&, const GridPosition&) const;
  virtual double GetDustRadius();
#ifndef ONE_DIMENSIONAL
  virtual int GetNumberDust();
  virtual GridPosition GetDustPosition(int);
#endif
#ifdef ONE_DIMENSIONAL
  virtual double GetBoxVolume () const {return (GetGridSize()[0]-1)*GetVolumeQuant();};
  virtual bool outside(const Position&) const;  
  virtual bool insidedust(const Position&) const;  
  virtual bool outsideRight(const Position&) const;
#endif
#ifdef TWO_DIMENSIONAL
  virtual double GetBoxVolume () const {return (GetGridSize()[0]-1)*(GetGridSize()[1]-1)*GetVolumeQuant();};
#endif
#ifdef THREE_DIMENSIONAL
  virtual double GetBoxVolume () const {
    return (GetGridSize()[0]-1)*(GetGridSize()[1]-1)*(GetGridSize()[2]-1)*GetVolumeQuant();
  };
#endif
#ifndef ONE_DIMENSIONAL
  virtual bool outside(const Position&) const;
  virtual bool outside_periodic(const Position&) const; 
  virtual bool outside_right(const Position&) const;
  virtual bool outside_left(const Position&) const;
  virtual bool insidedust(const Position&, GridPosition) const;
#endif
  virtual void DivByVolume(grid<double> &G);
  virtual void DivByVolume(grid<Position> &G);
  virtual void DivByVolume(grid<Velocity> &G);
 protected:
  virtual PARAMETERMAP* MakeParamMap (PARAMETERMAP* pm = NULL);  
};
#endif