particleranvd.cpp

Go to the documentation of this file.

#include <particleranvd.h>


particleranvd::particleranvd
(
  doublec radius_,
  doublec vmax_, 
  boxcollision const & box_
)
{
  radius = radius_;
  vmax = vmax_;
  box = box_;
}


void particleranvd::eval(particle & p) const
{
  assert(radius!=0.0);
  assert(vmax!=0.0);

  doublec t = 2.0*3.141592*r11();
  doublec rad = vmax*r11();
  p.vel[0] = rad*cos(t);
  p.vel[1] = rad*sin(t);

  p.radius = radius;

  randomposition(p);
}

particlev0::particlev0
(
  doublec radius_,
  doublec vmax_, 
  boxcollision const & box_
)
{
  radius = radius_;
  vmax = vmax_;
  box = box_;
}

void particlev0::eval(particle & p) const
{
  assert(radius!=0.0);
  assert(vmax!=0.0);

  doublec t = 2.0*3.141592*r11();
  doublec rad = vmax;
  p.vel[0] = rad*cos(t);
  p.vel[1] = rad*sin(t);

  p.radius = radius;

  randomposition(p);
}

particlev0theta0::particlev0theta0
(
  doublec radius_,
  doublec vmax_, 
  boxcollision const & box_
)
{
  radius = radius_;
  vmax = vmax_;
  box = box_;
}

void particlev0theta0::eval(particle & p) const
{
  assert(radius!=0.0);
  assert(vmax!=0.0);

  doublec t = theta0;
  doublec rad = vmax;
  p.vel[0] = rad*cos(t);
  p.vel[1] = rad*sin(t);

  p.radius = radius;

  randomposition(p);
}


particlev0spaced::particlev0spaced 
(
  doublec radius_,
  doublec vmax_, 
  boxcollision const & box_
)
{
  radius = radius_;
  vmax = vmax_;
  box = box_;
}


void particlev0spaced::eval(particle & p) const
{
  assert(radius!=0.0);
  assert(vmax!=0.0);


  double t = 2.0*3.141592*r11();
  doublec rad = vmax;
  p.vel[0] = rad*cos(t);
  p.vel[1] = rad*sin(t);

  p.radius = radius;


  // Even spaced algorithm

  static uint c=0;

  ++c;
   
  uint k;
  uint n;

  double sx(0.0);
  double sy(0.0);
  t=0.5;

  for ( n=c; n!=0; n/=4 )
  {
    k = n % 4;

    switch (k)
    {
      case 0: 
        break;
      case 1:
        sx += t; 
        break;
      case 2:
        sx += t;
        sy += t; 
        break;
      case 3:
        sy += t;
    }
    t *= .5;
  }

  p.pos[0] = box.x0 + (box.x1-box.x0)*(sx+.1);
  p.pos[1] = box.y0 + (box.y1-box.y0)*(sy+.1);
}