CellSearch.hpp

#ifndef __FIBER_BASEOP_CELLSEARCH_HPP
#define __FIBER_BASEOP_CELLSEARCH_HPP
 
 
#include "gridopDllApi.h"
 
#include <field/UniformCartesianArray.hpp>
#include <field/Field.hpp>
 
#include <aerie/BoundingBox.hpp>
 
namespace Fiber
{
 
struct  gridop_API      CellSearch
{
        RefPtr<UniformCartesianArray>   UniGridCoords; 
 
        typedef uint16_t  idx_t;
 
        typedef FixedArray<idx_t, 3>            CellIndex_t; 
 
        RefPtr<MemArray<4, idx_t>       >       FragmentLookup;
 
        RefPtr<MemArray<4, CellIndex_t> >       CellLookup;
 
        RefPtr<MemArray<4, index_t>     >       LinearCellLookup;
 
        typedef MultiArray<4, idx_t>            FragmentLookup_t;
        typedef MultiArray<4, CellIndex_t>      CellLookup_t;
        typedef MultiArray<4, index_t>          LinearCellLookup_t; 
 
 
        CellSearch(const Field&Coordinates   , 
                   const Eagle::BoundingBox&BBox, 
                   const MultiIndex<3>&UnigridSize,
                   index_t MaxCellsPerUnicell, const WeakPtr<FragmentSelector>&FS = NullPtr() );
 
 
        bool    addToUniCell(FragmentLookup_t&FL, CellLookup_t&CL,
                             const  MultiIndex<3>&UniCell, int FragmentIndex, const MultiIndex<3>&CellIndex); 
 
 
        struct  Location
        {
                int                     FragmentID;
                CellIndex_t             CellIndex;
                FixedArray<double,3>    VertexIndex;
 
                Location(int theFragmentID)
                : FragmentID( theFragmentID )
                {
                        CellIndex[0] = CellIndex[1] = CellIndex[2] = -1; 
                        VertexIndex[0] = VertexIndex[1] = VertexIndex[2] = 0.0;
                }
 
                Location(int theFragmentID, const CellIndex_t&theCellIndex,
                         const FixedArray<double,3>&theVertexIndex)
                : FragmentID( theFragmentID )
                , CellIndex( theCellIndex )
                , VertexIndex( theVertexIndex )
                {}
        };
 
        Location findLocation(const Eagle::PhysicalSpace::point&P, const Field&Coordinates) const; 
};
 
 
} // namespace Fiber
 
 
#endif // __FIBER_BASEOP_CELLSEARCH_HPP