Point3DSet.cpp

Go to the documentation of this file.

#include "Point3DSet.h"
#include "PCSIMException.h"

#include <iostream>
using std::cerr;
using std::endl;

size_t Point3DSet::getIndex( Point3D<double> const& p )
{
  pointhash_t::iterator i = lookup.find( p );
  if ( i == lookup.end() ) {
     throw( PCSIM::NotFoundException( "Point3DSet::getIndex", "Point not found." ) );
  }
  return i->second;
}

size_t Point3DSet::getIndex(double x, double y, double z)
{
  return getIndex( Point3D<double>(x,y,z) );
}

shared_ptr< Point3DSet > Point3DSet::subspace( SubSpaceConstraint const& cnstr ) const
{
  return subset( *subspace_index( cnstr ) );
}

shared_ptr< vector<size_t> >  Point3DSet::subspace_index( SubSpaceConstraint const& cnstr ) const
{
  shared_ptr< vector<size_t> > idx( new vector<size_t> );
  return idx;
}

shared_ptr< Point3DSet > Point3DSet::subset( vector<size_t> const& indices ) const
{
   shared_ptr< Point3DSet > result( new Point3DSet( indices.size() ) );
   for( size_t i=0; i<indices.size(); i++ ) {
        size_t s = indices[i];
        if( s >= m_points.size() ) throw( PCSIM::Exception( "Point3DSet::subset", "Index out of range.") );
        result->setPoint( i, m_points[ s ] );
   }
   return result;
}

---