One-dimensional interpolated (e.g. More...

#include <PrimitiveFunctions.h>

Inheritance diagram for inet::math::PeriodicallyInterpolated1DFunction< R, X >:

Public Member Functions
	PeriodicallyInterpolated1DFunction (X start, X end, const IInterpolator< X, R > &interpolator, const std::vector< R > &rs)

virtual R	getValue (const Point< X > &p) const override

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

Public Member Functions inherited from inet::math::FunctionBase< R, Domain< X > >
virtual void	partition (const typename Domain< X > ::I &i, const std::function< void(const typename Domain< X > ::I &, const IFunction< R, Domain< X > > *)> 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 > ::I &i) const override
	Returns the valid range of the function as an interval for the given domain. More...

virtual Domain< X > ::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 > ::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 > ::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 > ::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 > ::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 > ::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 > ::I &i) const override
	Returns the integral value for the given domain. More...

virtual const Ptr< const IFunction< R, Domain< X > > >	add (const Ptr< const IFunction< R, Domain< X > >> &o) const override
	Adds the provided function to this function. More...

virtual const Ptr< const IFunction< R, Domain< X > > >	subtract (const Ptr< const IFunction< R, Domain< X > >> &o) const override
	Substracts the provided function from this function. More...

virtual const Ptr< const IFunction< R, Domain< X > > >	multiply (const Ptr< const IFunction< double, Domain< X > >> &o) const override
	Multiplies the provided function with this function. More...

virtual const Ptr< const IFunction< double, Domain< X > > >	divide (const Ptr< const IFunction< R, Domain< X > >> &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 > ::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 > ::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 > ::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 > >
virtual	~IFunction ()

virtual R	getValue (const typename Domain< X > ::P &p) const=0
	Returns the value of the function at the given point. More...

Protected Attributes
const X	start

const X	end

const X	step

const IInterpolator< X, R > &	interpolator

const std::vector< R >	rs

Detailed Description

template<typename R, typename X>
class inet::math::PeriodicallyInterpolated1DFunction< R, X >

One-dimensional interpolated (e.g.

constant, linear) function between intervals defined by points positioned periodically on the X axis.

Constructor & Destructor Documentation

◆ PeriodicallyInterpolated1DFunction()

template<typename R , typename X >

inet::math::PeriodicallyInterpolated1DFunction< R, X >::PeriodicallyInterpolated1DFunction	(	X	start,
		X	end,
		const IInterpolator< X, R > &	interpolator,
		const std::vector< R > &	rs
	)

inline

526 :

527 start(start), end(end), step((end - start) / (rs.size() - 1)), interpolator(interpolator), rs(rs) {}

Member Function Documentation

◆ getValue()

template<typename R , typename X >

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

inlineoverridevirtual

                                                          {
         X x = std::get<0>(p);
         int index = std::floor(toDouble(x - start) / toDouble(step));
         if (index < 0 || index > rs.size() - 2)
             return R(0);
         else {
             R r1 = rs[index];
             R r2 = rs[index + 1];
             X x1 = start + step * index;
             X x2 = x1 + step;
             return interpolator.getValue(x1, r1, x2, r2, x);
         }
     }

◆ partition()

template<typename R , typename X >

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

inlineoverridevirtual

                                                                                                                                                        {
         const auto& i1 = i.getIntersected(Interval<X>(getLowerBound<X>(), Point<X>(start), 0b0, 0b0, 0b0));
         if (!i1.isEmpty()) {
             ConstantFunction<R, Domain<X>> g(R(0));
             callback(i1, &g);
         }
         const auto& i2 = i.getIntersected(Interval<X>(Point<X>(start), Point<X>(end), 0b1, 0b0, 0b0));
         if (!i2.isEmpty()) {
             int startIndex = std::max(0, (int)std::floor(toDouble(std::get<0>(i.getLower()) - start) / toDouble(step)));
             int endIndex = std::min((int)rs.size() - 1, (int)std::ceil(toDouble(std::get<0>(i.getUpper()) - start) / toDouble(step)));
             for (int index = startIndex; index < endIndex; index++) {
                 Point<X> startPoint(start + step * index);
                 Point<X> endPoint(start + step * (index + 1));
                 const auto& i3 = i.getIntersected(Interval<X>(startPoint, endPoint, 0b1, 0b0, 0b0));
                 if (!i3.isEmpty()) {
                     if (dynamic_cast<const ConstantInterpolatorBase<X, R> *>(&interpolator)) {
                         R r = getValue((startPoint + endPoint) / 2);
                         ConstantFunction<R, Domain<X>> g(r);
                         callback(i3, &g);
                     }
                     else
                         throw cRuntimeError("TODO");
                 }
             }
         }
         const auto& i4 = i.getIntersected(Interval<X>(Point<X>(end), getUpperBound<X>(), 0b1, 0b0, 0b0));
         if (!i4.isEmpty()) {
             ConstantFunction<R, Domain<X>> g(R(0));
             callback(i4, &g);
         }
     }