pointsurface.h

Go to the documentation of this file.

#ifndef POINTSURFACE_H
#define POINTSURFACE_H

#include <cassert>

#include <point.h>
#include <gobj.h>

#include <random.h>


template< typename T=double, typename G=randomgenerator >
class pointsurface : public gobjContainerPrePost
{
  pointsurface()
    { assert(false); }
public:

  uint N;

  random11<T,G> rnd;

  pointsurface(uintc N_)
    : N(N_)
  {
    pre.push( new gobjglPushAttrib(GL_LIGHTING) );
    pre.push( new gobjglPushAttrib(GL_CURRENT_BIT) );
    pre.push( new gobjglDisable(GL_LIGHTING) );

    post.push( new gobjglPopAttrib() );
    post.push( new gobjglPopAttrib() );
  }

  template< typename FN >
  void addsurface2D( FN & fn )
  {
    point3<double> p;
    bool valid(false);

    nuke();

    push( new gobjglBegin(GL_POINTS) );

    for (uint i=0; i<N; )
    {
      point3<double> p;
      fn(valid,p.x,p.y,p.z,rnd(),rnd());
      if (valid)
      {
        push( new gobjglVertex3d(p) );
        ++i;
      }
    }

    push( new gobjglEnd() );
  }

  template< typename FN >
  void addsurface3D( FN & fn )
  {
    point3<double> p;
    bool valid(false);

    nuke();

    push( new gobjglBegin(GL_POINTS) );

    for (uint i=0; i<N; )
    {
      point3<double> p;
      fn(valid,p.x,p.y,p.z,rnd(),rnd(),rnd());
      if (valid)
      {
        push( new gobjglVertex3d(p) );
        ++i;
      }
    }

    push( new gobjglEnd() );
  }

};

// <TODO> Template the glVertex3d so that the type T determines which
//  OpenGL call is made. ie double -> glVertex3d, float -> glVertex3f.

class pointsurfaceSphere
{
public:

  void operator()
  (
    bool & accept,
    double & x,
    double & y,
    double & z,
    doublec u,
    doublec v,
    doublec w 
  );

};

class pointsurfaceParallelogram
{
public:

  point3<double> p0;

  point3<double> pu;
  point3<double> pv;

  pointsurfaceParallelogram
  (
    point3<double> const & p0_,
    point3<double> const & p1,
    point3<double> const & p2
  )
    : p0(p0_), pu(p1-p0), pv(p2-p0) {}

  void operator()
  (
    bool & accept,
    double & x,
    double & y,
    double & z,
    doublec u,
    doublec v
  ) const
  {
    accept = true;
    x = p0.x + pu.x*u + pv.x*v;
    y = p0.y + pu.y*u + pv.y*v;
    z = p0.z + pu.z*u + pv.z*v;
  }
};


#endif