NPlane.cpp

Go to the documentation of this file.

#include "Mint/NPlane.h"

NPlane::NPlane(const HyperPointSet& a) :
  _origin ( a.at(0)          ),
  _v      ( a.getDimension() )
{
  
  WELCOME_LOG << "Hello from the NPlane() Constructor";
  
  for(unsigned i = 1; i < a.size(); i++){
    _v.push_back( a.at(i) - a.at(0) );
  }

  if ( _v.getDimension() < (int)_v.size() ) ERROR_LOG << "You are providing too many directions in this dimension - your points must be L.D.";
  if ( 0                == (int)_v.size() ) ERROR_LOG << "You have not provided any directions";


}

NPlane::NPlane(const HyperPoint& origin, const HyperPointSet& v) :
  _origin ( origin ),
  _v      ( v      )
{
  WELCOME_LOG << "Hello from the NPlane() Constructor";

  if ( _v.getDimension() < (int)_v.size() ) ERROR_LOG << "You are providing too many directions in this dimension - your points must be L.D.";
  if ( 0                == (int)_v.size() ) ERROR_LOG << "You have not provided any directions";

}


HyperPoint NPlane::getParametricPoint(const HyperPoint& t) const{
  if ( t.size() != (int)_v.size() ) ERROR_LOG << "You have not provided the correct number of parameters";

  HyperPoint point = _origin;

  for(unsigned int i = 0; i < _v.size(); i++){
    point = point + _v.at(i)*t.at(i); //std::vector times a double
  }

  return point;

}

void NPlane::print(std::ostream& os, int endline) const{
  os << "NPlane: v = ";
  _origin.print(os, 0);

  for (unsigned int i = 0; i < _v.size(); i++){
    os << " + ";
    _v.at(i).print(os, 0);
    os << "t_" << i;
  }
  if (endline) os << std::endl;
}

NPlane::~NPlane(){

}