123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191 |
-
-
-
-
-
-
- #ifndef CUBICSPLINEINTERPOLATOR_INC
- #define CUBICSPLINEINTERPOLATOR_INC
-
- #include "LemmaObject.h"
-
- namespace Lemma {
-
-
- struct SplineSet{
- VectorXr a;
- VectorXr b;
- VectorXr c;
- VectorXr d;
- VectorXr x;
-
- SplineSet( ) {
- }
-
- SplineSet(const int&n) {
- a = VectorXr::Zero(n+1);
- b = VectorXr::Zero(n);
- c = VectorXr::Zero(n+1);
- d = VectorXr::Zero(n);
- x = VectorXr::Zero(n+1);
- }
- };
-
-
- class CubicSplineInterpolator : public LemmaObject {
-
- friend std::ostream &operator<<(std::ostream &stream,
- const CubicSplineInterpolator& ob);
-
-
-
- public:
-
-
-
-
- explicit CubicSplineInterpolator ( const ctor_key& );
-
-
- CubicSplineInterpolator ( const YAML::Node& node, const ctor_key& );
-
-
- virtual ~CubicSplineInterpolator ();
-
-
-
- static std::shared_ptr<CubicSplineInterpolator> NewSP();
-
-
-
- virtual YAML::Node Serialize() const;
-
-
-
- static std::shared_ptr< CubicSplineInterpolator > DeSerialize(const YAML::Node& node);
-
-
-
-
-
-
-
- void SetKnots(const VectorXr& x, const VectorXr& y);
-
-
-
- void ResetKnotOrdinate( const VectorXr& y );
-
-
-
- VectorXr InterpolateOrderedSet(const VectorXr& x);
-
-
-
- Real Integrate(const Real& x0, const Real& x1, const int& n);
-
-
-
- Real Integrate(const Real& x0, const Real& x1);
-
-
-
- VectorXr GetKnotAbscissa();
-
-
-
- VectorXr GetKnotOrdinate();
-
-
-
- Real Interpolate(const Real& x, int& i);
-
-
-
- Real Interpolate(const Real& x);
-
-
-
-
-
-
- virtual inline std::string GetName() const {
- return CName;
- }
-
- protected:
-
-
-
-
-
- int Interval(const Real& x);
-
- private:
-
-
- CubicSplineInterpolator( const CubicSplineInterpolator& ) = delete;
-
-
- static constexpr auto CName = "CubicSplineInterpolator";
-
- SplineSet Spline;
-
- int ilo;
-
- int mflag;
-
-
-
- };
-
- }
-
- #endif
-
-
|