TreeCellMapper.hpp

#ifndef __TREECELLMAPPER_HPP
#define __TREECELLMAPPER_HPP
 
#include <aerie/BoundingBox.hpp>
#include <aerie/KDTree.hpp>
#include <eagle/PhysicalSpace.hpp>
#include <eagle/QuadraticMatrix.hpp> 
 
#include <field/UniformCartesianArray.hpp>
 
#include <fish/fiber/vector/Interpolate.hpp>
#include <fish/fiber/vector/LinearIpol.hpp>
 
#include <memcore/RefPtr.hpp>
 
using namespace Fiber;
using namespace Eagle;
using namespace std;
 
using namespace Eagle::PhysicalSpace;
 
class TreeCellMapper : public ReferenceBase<TreeCellMapper>
{
RefPtr<MemArray<3, point> > m_crds_arr;
MultiArray<3, point>&       m_crds;
 
RefPtr<Eagle::KDTree<3, MultiIndex<3> > > m_tree;
        
double                      m_eps; 
unsigned                    m_spread;
double                      m_max_len;
tvector                     m_range_scale;
MultiIndex<3>               m_ix, m_iy, m_iz, m_max;
tvector                     m_max_vec;
 
 
public:
 
        TreeCellMapper( RefPtr<MemArray<3, point> > crds, tvector range_scale = tvector(1.,1.,1.), double epsilon = 1./50, unsigned spread = 13 );
        
 
static inline double maxi(const double a, const double b)
{
        return ( a < b ) ? b:a;
}
 
static inline void maxVec( double& max_len, tvector& max_vec, const tvector& v )
{
tvector tmp(0,0,0); 
double len  = norm(v);
 
        for(unsigned i = 0; i < tvector::Dims; i++)
        {
                tmp[i] = abs( v[i] ); 
                max_vec[i] =  tmp[i] > max_vec[i] ? tmp[i] : max_vec[i];
        } 
 
        if( len > max_len  )
                max_len = len;
}
        
        unsigned localCoordinates( point&uvw, tvector& cell_size, const point&p );
 
        void setQueryScale( tvector& scale )
        {
                m_range_scale = scale;
        }
 
}; 
 
 
 
#endif // __UNIGRIDMAPPER_HPP