UniformCartesianArray.hpp

#ifndef __FIBER__UNIFORMCARTESIANARRAY_HPP
#define __FIBER__UNIFORMCARTESIANARRAY_HPP

#include "DirectProductArray.hpp"
#include <eagle/PhysicalSpace.hpp>
#include <aerie/BoundingBox.hpp>

namespace Fiber
{

class   FIELD_API       UniformCartesianArray : public  LinearDirectProductMemArray<Eagle::point3>
{
        Eagle::point3 MakeDeltaPoint(const Eagle::tvector3&what)
        {
                return Eagle::point3(0.0,0.0,0.0) + what;
        } 

        Eagle::point3 DeltaPoint(const Eagle::tvector3&what, const MultiIndex<3>&divisor)
        {
        Eagle::point3 retval; 
                for(int i=0; i<3; i++)
                {
                        if (divisor[i]<2)
                                throw Exception( i, divisor[i]);

                        retval[i] = what[i] / (divisor[i]-1);
                }
                return retval;
        }

public: 
        struct  Exception : std::exception
        {
                int     dimension;
                index_t value;

                Exception(int Dimension, index_t Value)
                : dimension( Dimension)
                , value(Value)
                {}

                override const char*what() const throw()
                {
                        return "UniformCartesianArray: Specified Grid has less than two points, can't construct such.";
                }               
        };

        UniformCartesianArray(const Eagle::point3&start, const MultiIndex<3>&dims, 
                              const Eagle::point3&end,
                              const MemBase::Creator_t&C = MemCore::NullPtr() ) throw (const Exception&)
        : MemBase(C)
        , LinearDirectProductMemArray<Eagle::point3>(dims, start, DeltaPoint(end-start, dims), C )
        {}

        UniformCartesianArray(const Eagle::BoundingBox&BBox, const MultiIndex<3>&dims, const MemBase::Creator_t&C = MemCore::NullPtr() ) throw (const Exception&)
        : MemBase(C)
        , LinearDirectProductMemArray<Eagle::point3>(dims, BBox.mincoord(), DeltaPoint(BBox.maxcoord()-BBox.mincoord(), dims), C )
        {} 

        UniformCartesianArray(const MultiIndex<3>&dims, const Eagle::point3&start,
                              const Eagle::tvector3&delta, const MemBase::Creator_t&C = MemCore::NullPtr() )
        : MemBase(C)
        , LinearDirectProductMemArray<Eagle::point3>(dims, start, MakeDeltaPoint(delta), C )
        {} 

        UniformCartesianArray(const MultiIndex<3>&dims, const Eagle::point3&start,
                              const Eagle::point3&delta, const MemBase::Creator_t&C = MemCore::NullPtr() )
        : MemBase(C)
        , LinearDirectProductMemArray<Eagle::point3>(dims, start, delta, C)
        {}

        ~UniformCartesianArray(); 

        Eagle::tvector3 cell_diagonal() const
        {
                return Delta() - Eagle::point3(0,0,0);
        }
};


} // namespace Fiber

#endif // __FIBER__UNIFORMCARTESIANARRAY_HPP