poisson.h

Go to the documentation of this file.

#ifndef POISSON_H
#define POISSON_H
#include "task.h"
#include "density.h"
#include "numeric.h"
#include "normalization.h"
#include "dimensionality.h"
#include "wall.h"
#include "potential.h"
#include <valarray>
//---------------------------------------------------------------------------------------------------------------------
#ifndef ONE_DIMENSIONAL

class Poisson {
 protected:
  double dx;
  double dy;
#ifdef THREE_DIMENSIONAL

  double dz;
#endif

  int nu1,nu2;
  int Rg;
  int nd;
  int ig;      
  double epsilon;
  std::valarray<double> potd;
  Geometry *geom;
  Wall *pWall;
  Potential *pPotential;
 public:

  Poisson(Geometry *geom_);
  
  void solve(grid<double> &u, grid<double> &f, double potd[20]);
 private:
  void solver(grid<double> &u, grid<double> &f, double potd[20], int time);

  void multigrid(grid<double> &u, grid<double> &f, double potd[20], int time);

  void gauss(grid<double> &u, grid<double> &f, double potd[20], int time);

  void itersteps(grid<double> &u, grid<double> &f, int i_start, int j_start, int i_stop, int j_stop, double potd[20], int time);
#ifdef THREE_DIMENSIONAL

  void itersteps(grid<double> &u, grid<double> &f, int i_start, int j_start, int k_start, int i_stop, int j_stop, int k_stop, double potd[20], int time);  
#endif

  void itersteps_fine(grid<double> &u, grid<double> &f, int i_start, int j_start, int i_stop, int j_stop, double potd[20], int time);
#ifdef THREE_DIMENSIONAL

  void itersteps_fine(grid<double> &u, grid<double> &f, int i_start, int j_start, int k_start, int i_stop, int j_stop, int k_stop, double potd[20], int time);
#endif

  void itersteps_coarse(grid<double> &u, grid<double> &f, double potd[20], int time);

  double distance(grid<double> &u, grid<double> &f);

  void boundary(grid<double> &u);

  void restriction(grid<double> &u, grid<double> &f, grid<double> &fn, int time);

  void prolongate(grid<double> &u, grid<double> &un, int t);
};
#endif
#endif