inet::math::PeriodicallyInterpolated2DFunction< R, X, Y > Class Template Reference

#include <PrimitiveFunctions.h>

Inheritance diagram for inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >:

Public Member Functions
	PeriodicallyInterpolated2DFunction (X startX, X endX, int sizeX, Y startY, Y endY, int sizeY, const IInterpolator< X, R > &interpolatorX, const IInterpolator< Y, R > &interpolatorY, const std::vector< R > &rs)
virtual R	getValue (const Point< X, Y > &p) const override
virtual void	partition (const Interval< X, Y > &i, const std::function< void(const Interval< X, Y > &, const IFunction< R, Domain< X, Y >> *)> callback) const override
Public Member Functions inherited from inet::math::FunctionBase< R, Domain< X, Y > >
virtual void	partition (const typename Domain< X, Y > ::I &i, const std::function< void(const typename Domain< X, Y > ::I &, const IFunction< R, Domain< X, Y > > *)> callback) const override
	Subdivides the provided domain and calls back f with the subdomains and the corresponding potentially simpler domain limited functions. More...
virtual Interval< R >	getRange () const override
	Returns the valid range of the function as an interval. More...
virtual Interval< R >	getRange (const typename Domain< X, Y > ::I &i) const override
	Returns the valid range of the function as an interval for the given domain. More...
virtual Domain< X, Y > ::I	getDomain () const override
	Returns the valid domain of the function as an interval. More...
virtual bool	isFinite () const override
	Returns true if the function value is finite in the whole domain. More...
virtual bool	isFinite (const typename Domain< X, Y > ::I &i) const override
	Returns true if the function value is finite in the given domain. More...
virtual bool	isNonZero () const override
	Returns true if the function value is non-zero in the whole domain. More...
virtual bool	isNonZero (const typename Domain< X, Y > ::I &i) const override
	Returns true if the function value is non-zero in the given domain. More...
virtual R	getMin () const override
	Returns the minimum value for the whole domain. More...
virtual R	getMin (const typename Domain< X, Y > ::I &i) const override
	Returns the minimum value for the given domain. More...
virtual R	getMax () const override
	Returns the maximum value for the whole domain. More...
virtual R	getMax (const typename Domain< X, Y > ::I &i) const override
	Returns the maximum value for the given domain. More...
virtual R	getMean () const override
	Returns the mean value for the whole domain. More...
virtual R	getMean (const typename Domain< X, Y > ::I &i) const override
	Returns the mean value for the given domain. More...
virtual R	getIntegral () const override
	Returns the integral value for the whole domain. More...
virtual R	getIntegral (const typename Domain< X, Y > ::I &i) const override
	Returns the integral value for the given domain. More...
virtual const Ptr< const IFunction< R, Domain< X, Y > > >	add (const Ptr< const IFunction< R, Domain< X, Y > >> &o) const override
	Adds the provided function to this function. More...
virtual const Ptr< const IFunction< R, Domain< X, Y > > >	subtract (const Ptr< const IFunction< R, Domain< X, Y > >> &o) const override
	Substracts the provided function from this function. More...
virtual const Ptr< const IFunction< R, Domain< X, Y > > >	multiply (const Ptr< const IFunction< double, Domain< X, Y > >> &o) const override
	Multiplies the provided function with this function. More...
virtual const Ptr< const IFunction< double, Domain< X, Y > > >	divide (const Ptr< const IFunction< R, Domain< X, Y > >> &o) const override
	Divides this function with the provided function. More...
virtual std::ostream &	printOn (std::ostream &os) const override
virtual void	print (std::ostream &os, int level=0) const override
	Prints this function in human readable form to the provided stream for the whole domain. More...
virtual void	print (std::ostream &os, const typename Domain< X, Y > ::I &i, int level=0) const override
	Prints this function in a human readable form to the provided stream for the given domain. More...
virtual void	printPartitioning (std::ostream &os, const typename Domain< X, Y > ::I &i, int level) const override
	Prints the partitioning of this function in a human readable form to the provided stream for the given domain. More...
virtual void	printPartition (std::ostream &os, const typename Domain< X, Y > ::I &i, int level=0) const override
	Prints a single partition of this function in a human readable form to the provided stream for the given domain. More...
virtual void	printStructure (std::ostream &os, int level=0) const override
	Prints the internal data structure of this function in a human readable form to the provided stream. More...
Public Member Functions inherited from inet::math::IFunction< R, Domain< X, Y > >
virtual	~IFunction ()
virtual R	getValue (const typename Domain< X, Y > ::P &p) const=0
	Returns the value of the function at the given point. More...

Protected Member Functions
virtual R	getValueInterpolatedAlongX (X x, int indexX, int indexY) const
void	call (const Interval< X, Y > &i, const std::function< void(const Interval< X, Y > &, const IFunction< R, Domain< X, Y >> *)> callback) const

Protected Attributes
const X	startX
const X	endX
const X	stepX
int	sizeX
const Y	startY
const Y	endY
const Y	stepY
int	sizeY
const IInterpolator< X, R > &	interpolatorX
const IInterpolator< Y, R > &	interpolatorY
const std::vector< R >	rs

Constructor & Destructor Documentation

PeriodicallyInterpolated2DFunction()

template<typename R , typename X , typename Y >

inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::PeriodicallyInterpolated2DFunction ( X  startX, X  endX, int  sizeX, Y  startY, Y  endY, int  sizeY, const IInterpolator< X, R > &  interpolatorX, const IInterpolator< Y, R > &  interpolatorY, const std::vector< R > &  rs  )

inline

                                                                                                                                                                                                             :
        startX(startX), endX(endX), stepX((endX - startX) / (sizeX - 1)), sizeX(sizeX),
        startY(startY), endY(endY), stepY((endY - startY) / (sizeY - 1)), sizeY(sizeY),
        interpolatorX(interpolatorX), interpolatorY(interpolatorY), rs(rs)
    {
        ASSERT(rs.size() == sizeX * sizeY);
    }

Member Function Documentation

call()

template<typename R , typename X , typename Y >

void inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::call ( const Interval< X, Y > &  i, const std::function< void(const Interval< X, Y > &, const IFunction< R, Domain< X, Y >> *)>  callback  ) const

inlineprotected

                                                                                                                                          {
        if (!i.isEmpty()) {
            ConstantFunction<R, Domain<X, Y>> g(R(0));
            callback(i, &g);
        }
    }

getValue()

template<typename R , typename X , typename Y >

virtual R inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::getValue ( const Point< X, Y > &  p ) const

inlineoverridevirtual

                                                            {
        X x = std::get<0>(p);
        Y y = std::get<1>(p);
        int indexX = std::floor(toDouble(x - startX) / toDouble(stepX));
        int indexY = std::floor(toDouble(y - startY) / toDouble(stepY));
        if (indexX < 0 || indexX > sizeX - 2 || indexY < 0 || indexY > sizeY - 2)
            return R(0);
        else {
            R r1 = getValueInterpolatedAlongX(x, indexX, indexY);
            R r2 = getValueInterpolatedAlongX(x, indexX, indexY + 1);
            Y y1 = startY + stepY * indexY;
            Y y2 = y1 + stepY;
            return interpolatorY.getValue(y1, r1, y2, r2, y);
        }
    }

getValueInterpolatedAlongX()

template<typename R , typename X , typename Y >

virtual R inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::getValueInterpolatedAlongX ( X  x, int  indexX, int  indexY  ) const

inlineprotectedvirtual

                                                                            {
        R r1 = rs[sizeX * indexY + indexX];
        R r2 = rs[sizeX * indexY + indexX + 1];
        X x1 = startX + stepX * indexX;
        X x2 = x1 + stepX;
        return interpolatorX.getValue(x1, r1, x2, r2, x);
    }

partition()

template<typename R , typename X , typename Y >

virtual void inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::partition ( const Interval< X, Y > &  i, const std::function< void(const Interval< X, Y > &, const IFunction< R, Domain< X, Y >> *)>  callback  ) const

inlineoverridevirtual

                                                                                                                                                                {
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(getLowerBound<X>(), getLowerBound<Y>()), Point<X, Y>(X(startX), Y(startY)), 0b00, 0b00, 0b00)), callback);
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(X(startX), getLowerBound<Y>()), Point<X, Y>(X(endX), Y(startY)), 0b10, 0b00, 0b00)), callback);
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(X(endX), getLowerBound<Y>()), Point<X, Y>(getUpperBound<X>(), Y(startY)), 0b10, 0b00, 0b00)), callback);
 
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(getLowerBound<X>(), Y(startY)), Point<X, Y>(X(startX), Y(endY)), 0b01, 0b00, 0b00)), callback);
        const auto& i1 = i.getIntersected(Interval<X, Y>(Point<X, Y>(X(startX), Y(startY)), Point<X, Y>(X(endX), Y(endY)), 0b11, 0b00, 0b00));
        if (!i1.isEmpty()) {
            int startIndexX = std::max(0, (int)std::floor(toDouble(std::get<0>(i.getLower()) - startX) / toDouble(stepX)));
            int endIndexX = std::min(sizeX - 1, (int)std::ceil(toDouble(std::get<0>(i.getUpper()) - startX) / toDouble(stepX)));
            for (int indexX = startIndexX; indexX < endIndexX; indexX++) {
                int startIndexY = std::max(0, (int)std::floor(toDouble(std::get<1>(i.getLower()) - startY) / toDouble(stepY)));
                int endIndexY = std::min(sizeY - 1, (int)std::ceil(toDouble(std::get<1>(i.getUpper()) - startY) / toDouble(stepY)));
                for (int indexY = startIndexY; indexY < endIndexY; indexY++) {
                    Point<X, Y> startPoint(startX + stepX * indexX, startY + stepY * indexY);
                    Point<X, Y> endPoint(startX + stepX * (indexX + 1), startY + stepY * (indexY + 1));
                    const auto& i3 = i.getIntersected(Interval<X, Y>(startPoint, endPoint, 0b11, 0b00, 0b00));
                    if (!i3.isEmpty()) {
                        if (dynamic_cast<const ConstantInterpolatorBase<X, R> *>(&interpolatorX) && dynamic_cast<const ConstantInterpolatorBase<Y, R> *>(&interpolatorY)) {
                            R r = getValue((startPoint + endPoint) / 2);
                            ConstantFunction<R, Domain<X, Y>> g(r);
                            callback(i3, &g);
                        }
                        else
                            throw cRuntimeError("TODO");
                    }
                }
            }
        }
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(X(endX), Y(startY)), Point<X, Y>(getUpperBound<X>(), Y(endY)), 0b11, 0b00, 0b00)), callback);
 
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(getLowerBound<X>(), Y(endY)), Point<X, Y>(X(startX), getUpperBound<Y>()), 0b01, 0b00, 0b00)), callback);
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(X(startX), Y(endY)), Point<X, Y>(X(endX), getUpperBound<Y>()), 0b11, 0b00, 0b00)), callback);
        call(i.getIntersected(Interval<X, Y>(Point<X, Y>(X(endX), Y(endY)), Point<X, Y>(getUpperBound<X>(), getUpperBound<Y>()), 0b11, 0b00, 0b00)), callback);
    }

Member Data Documentation

endX

template<typename R , typename X , typename Y >

const X inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::endX

protected

endY

template<typename R , typename X , typename Y >

const Y inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::endY

protected

interpolatorX

template<typename R , typename X , typename Y >

const IInterpolator<X, R>& inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::interpolatorX

protected

interpolatorY

template<typename R , typename X , typename Y >

const IInterpolator<Y, R>& inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::interpolatorY

protected

rs

template<typename R , typename X , typename Y >

const std::vector<R> inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::rs

protected

sizeX

template<typename R , typename X , typename Y >

int inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::sizeX

protected

sizeY

template<typename R , typename X , typename Y >

int inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::sizeY

protected

startX

template<typename R , typename X , typename Y >

const X inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::startX

protected

startY

template<typename R , typename X , typename Y >

const Y inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::startY

protected

stepX

template<typename R , typename X , typename Y >

const X inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::stepX

protected

stepY

template<typename R , typename X , typename Y >

const Y inet::math::PeriodicallyInterpolated2DFunction< R, X, Y >::stepY

protected

---

The documentation for this class was generated from the following file:

• PrimitiveFunctions.h