aslInterpolation.h
Go to the documentation of this file.
/*
* Advanced Simulation Library <http://asl.org.il>
*
* Copyright 2015 Avtech Scientific <http://avtechscientific.com>
*
*
* This file is part of Advanced Simulation Library (ASL).
*
* ASL is free software: you can redistribute it and/or modify it
* under the terms of the GNU Affero General Public License as
* published by the Free Software Foundation, version 3 of the License.
*
* ASL is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with ASL. If not, see <http://www.gnu.org/licenses/>.
*
*/
#ifndef ASLINTERPOLATION
#define ASLINTERPOLATION
#include "../aslUtilities.h"
#include <cmath>
namespace asl
{
class ElementBase;
typedef std::shared_ptr<ElementBase> Element;
template <class Func> class UniversalFunction
{
public:
inline double operator()(double x);
inline double operator()(Element x);
};
class LinearSpline
{
public:
template <typename T> inline T operator(T x,T r0)
{
return fabs(x)<r0 ? 1.- fabs(x)/r0 : 0;
}
};
class QuadraticSpline
{
public:
template <typename T> inline T operator(T x,T r0)
{
T a(fabs(x));
T b((1.- a/r0));
return a<r0 ? b*b : 0;
}
};
class QubicSpline
{
public:
template <typename T> inline T operator(T x,T r0)
{
T a(fabs(x));
T b((1.- a/r0));
return a<r0 ? b*b*b : 0;
}
};
class QubicSpline1
{
public:
template <typename T> inline T operator(T x,T r0)
{
T a(fabs(x));
T b(a/r0);
T b2(b*b);
return b < 1 ? 2. * b2*b - 3. * b2 + 1 : 0;
}
};
// --------------------------- Implementation ---------------------
template <class Func> inline double UniversalFunction<Func>::operator()(double x)
{
return Func(x);
}
template <class Func> inline Element UniversalFunction<Func>::operator()(Element x)
{
return Func(x);
}
} // asl
#endif
asl::LinearSpline
Linear spline function.
Definition: aslInterpolation.h:54
asl::LinearSpline::operator
T operator(T x, T r0)
Definition: aslInterpolation.h:56
asl::QuadraticSpline
Quadratic spline function.
Definition: aslInterpolation.h:70
asl::QuadraticSpline::operator
T operator(T x, T r0)
Definition: aslInterpolation.h:72
asl::QubicSpline1
Qubic spline function.
Definition: aslInterpolation.h:106
asl::QubicSpline1::operator
T operator(T x, T r0)
Definition: aslInterpolation.h:108
asl::QubicSpline
Qubic spline function.
Definition: aslInterpolation.h:88
asl::QubicSpline::operator
T operator(T x, T r0)
Definition: aslInterpolation.h:90
asl::UniversalFunction
Definition: aslInterpolation.h:40
asl::UniversalFunction::operator()
double operator()(double x)
Definition: aslInterpolation.h:119
acl::elementOperators::fabs
Element fabs(Element a)
asl
Advanced Simulation Library.
Definition: aslDataInc.h:31
asl::Element
std::shared_ptr< ElementBase > Element
Definition: aslInterpolation.h:36