DataInterpolation.h

#ifndef gridripper_amr1d_interpolation_DataInterpolation_h
#define gridripper_amr1d_interpolation_DataInterpolation_h


#include <gridripper/tvalarray.h>
#include <gridripper/config.h> /* GReal_t */
#include <gridripper/lang/IllegalArgumentException.h>


namespace gridripper { namespace amr1d {

class PDE;

namespace interpolation {

class DataInterpolation
{
    private:
        tvalarray<int> outKxArray;
        tvalarray<GReal_t> outData;

    public:

        DataInterpolation(const tvalarray<GReal_t>& data, const tvalarray<int>& kxarr, int nc, GReal_t deltaX,
                             int mul, int order, const PDE *pde, bool even)
                throw(IllegalArgumentException&);

        ~DataInterpolation();

        const tvalarray<GReal_t>& getData() const;

        const tvalarray<int>& getKxArray() const;

    private:
        static tvalarray<GReal_t> extend(const tvalarray<GReal_t>& data, int nc, GReal_t dx, int mul,
                                        int order, const PDE *pde, bool even);

        static tvalarray<int> extend(const tvalarray<int>& kxarr, int mul, int order,
                                        const PDE *pde, bool even);

        static tvalarray<GReal_t> restrict(const tvalarray<GReal_t>& data, int nc, int mul, int order,
                                        const PDE *pde, bool even);

        static tvalarray<int> restrict(const tvalarray<int>& kxarr, int mul, int order,
                                        const PDE *pde, bool even);

        static int getLeftMargin(const PDE *pde, bool even, int order);

        static int getRightMargin(const PDE *pde, bool even, int order);

        static int getMarginSize(int order) throw(IllegalArgumentException&);
};


} } } // namespace gridripper::amr1d::interpolation


#endif /* gridripper_amr1d_interpolation_DataInterpolation_h */