particleranvd.h

Go to the documentation of this file.

#ifndef PARTICLERANVD_H
#define PARTICLERANVD_H

#include <cassert>
#include <cmath>
using namespace std;

#include <particle.h>
#include <boxcollision.h>
#include <random.h>


class particledistribution
{
protected:
  random11<> r11;

public:

  double radius;
  double vmax;
  double theta0;

  boxcollision box;

  particledistribution()
    : radius(0.0), vmax(0.0), theta0(0.0),
      box(0.0,0.0,0.0,0.0) {}

  void randomposition(particle & p) const
  {
    doublec pradius(p.radius);
    p.pos[0] = box.x0 + pradius + (box.x1-box.x0-2.0*pradius)*r11();
    p.pos[1] = box.y0 + pradius + (box.y1-box.y0-2.0*pradius)*r11();
  }
};


class particleranvd : public particledistribution 
{
public:

  particleranvd
  (
    doublec radius_,
    doublec vmax_, 
    boxcollision const & box_
  );
  
  void eval(particle & p) const;
};

class particlev0 : public particledistribution
{
public:

  particlev0
  (
    doublec radius_,
    doublec vmax_, 
    boxcollision const & box_
  );
  
  void eval(particle & p) const;
};

class particlev0theta0 : public particledistribution
{
public:

  particlev0theta0
  (
    doublec radius_,
    doublec vmax_, 
    boxcollision const & box_
  );
  
  void eval(particle & p) const;
};

class particlev0spaced : public particledistribution
{
public:

  particlev0spaced 
  (
    doublec radius_,
    doublec vmax_, 
    boxcollision const & box_
  );
  
  void eval(particle & p) const;
};

#endif