potential.h

Go to the documentation of this file.

#ifndef POTENTIAL_H
#define POTENTIAL_H
#include "task.h"
#include "numeric.h"
#include "geometry.h"
#include "density.h"
#include "wall.h"
#define ID_POTENTIAL "Potential"
//---------------------------------------------------------------------------------------------------------------------

class Potential : public Task {
 private:
  std::string strGeometry;
 protected:
  int tstep;
  GridPosition NG;
  double dx;
  double Rd;
  double cd;
  grid<double> Pot;
  Geometry* pGeom;
 public:  
  Potential () : pGeom(NULL) {};
  virtual ~Potential () {};
#ifdef ONE_DIMENSIONAL

  double Vp;
#endif
#ifndef ONE_DIMENSIONAL

  double potd[20];
#endif

  virtual std::string GetClassName () const {return ID_POTENTIAL;}
  virtual grid<double>* GetGrid () {return &Pot;}
  virtual Position GetSpaceQuant () const {return dx;}
  virtual double& operator[] (GridPosition i) {return Pot[i];}
  virtual GridPosition GetNx () const {return NG;}
  virtual void Init ();
 protected:
  virtual PARAMETERMAP* MakeParamMap (PARAMETERMAP* pm);
};

class DustPotential : public Potential {
 private:
  std::string strDensity;
  double xq1;
  double xq2;
 protected:
  std::valarray<double> Nd_u;
  std::valarray<double> Nd_o;
  std::valarray<double> Hd;
  std::valarray<double> ni;
  Density* pRho;
  Wall* pWall;
  double n0;
 public:
  DustPotential () {}
  virtual ~DustPotential () {}
  virtual void Init ();
  virtual bool Execute ();
 protected:
  virtual PARAMETERMAP* MakeParamMap (PARAMETERMAP* pm = NULL);  
};
#endif