density.h

Go to the documentation of this file.

#ifndef DENSITY_H
#define DENSITY_H
#include "task.h"
#include "numeric.h"
#include "movement.h"
#include "species.h"
#include <vector>
#define ID_DENSITY "Density"
//#define ID_CURDENSITY "CurDensity"
//---------------------------------------------------------------------------------------------------------------------
class Geometry;

class Density : public Task {
 private:
  std::string strGeometry, strWall;

 protected:
  Geometry* pGeom;
  GridPosition NG;
  double dx;
  double dC;
  double dV;
  grid<double> den;
  std::vector<WSpeciesBase> species;
#ifdef ONE_DIMENSIONAL

  int Rg;
#endif
 public:  
  Density () {}
  virtual ~Density () {}
  std::string GetClassName () const {return ID_DENSITY;};
  
  void Init ();
  
  bool Execute ();
  void AddSpecies (SpeciesBase* pS);
  virtual grid<double>* GetGrid () {return &den;};
  virtual const GridPosition& GetNx () const {return NG;}; 
  double& operator[] (GridPosition i) {return den[i];}; 
 protected:
  PARAMETERMAP* MakeParamMap (PARAMETERMAP* pm = NULL);  
};
#endif