math.hpp

Go to the documentation of this file.

#ifndef PLX_MATH_HPP
#define PLX_MATH_HPP

#include <gmtl/gmtl.h>

typedef gmtl::Rayd Ray;
typedef gmtl::Matrix44d Matrix;

using namespace gmtl;


template<typename T> inline T clamp(T x, T min, T max) { return (x<min) ? min : (x>max) ? max : x; }

// typedef float Scalar;

inline Vec3d spherical_coordinate_to_Vec3d(double azimuth, double polar)
{
  return Vec3d(
    cos(Math::deg2Rad(azimuth)) * sin(Math::deg2Rad(polar) )
  , sin(Math::deg2Rad(azimuth)) * sin(Math::deg2Rad(polar) )
  , cos(Math::deg2Rad(polar))
  );
}

template<typename IterT>
unsigned char* new_colormap(int steps, IterT keys, int nkeys)
{
  static const int CHANNELS = 3;

  unsigned char* cmap = new unsigned char[steps*CHANNELS];
  
  unsigned char* out = cmap;
  double window = 1.0/nkeys;

  for(int i = 0; i < steps ; ++i )
  {
    // f is the global step interpolation fraction [ = 0..1 ]
    double f = ( double(i) / double(steps-1) );
    
    // alpha mix factor between two keys [ = 0..1  ]

    double keyf = f * double(nkeys-1);
    double base;
    double alpha = modf( keyf, &base );

    int ia = int(base);    
    if (ia >= nkeys) ia = nkeys-1;
    int ib = int(base) + 1;
    if (ib >= nkeys) ib = nkeys-1;

    Vec3d col = keys[ia] * (1-alpha) + keys[ib] * alpha;
    
    for(int j = 0; j < CHANNELS ; ++j )
    {
      out[j] = static_cast<unsigned char>(  col[j]*255 );
    }
    out += CHANNELS;
    /*
    *out++ =unsigned char(  col[0]*255 );
    *out++ =unsigned char(  col[1]*255 );
    *out++ =unsigned char(  col[2]*255 );
    */
  }
  return cmap;
}

inline Point3d center(const AABoxd& aabox)
{
  return ( aabox.getMin() + aabox.getMax() ) * 0.5;
}

inline Vec3d size(const AABoxd& aabox)
{
  return ( aabox.getMax() - aabox.getMin() );
}

#endif // PLX_MATH_HPP