LSST Applications g013ef56533+63812263fb,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+fde7a7a78c,g210f2d0738+db0c280453,g262e1987ae+abed931625,g29ae962dfc+058d1915d8,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+64337f1634,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5a60e81ecd+d9e514a434,g64539dfbff+db0c280453,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g7382096ae9+36d16ea71a,g74acd417e5+c70e70fbf6,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+1b779678e3,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+9583a964dd,g98a1a72a9c+028271c396,g98df359435+530b675b85,gb8cb2b794d+4e54f68785,gbf99507273+8c5ae1fdc5,gc2a301910b+db0c280453,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+c70e70fbf6,ge410e46f29+f459a6810c,ge41e95a9f2+db0c280453,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
Spline.h
Go to the documentation of this file.
1#if !defined(LSST_AFW_MATH_DETAIL_SPLINE)
2#define LSST_AFW_MATH_DETAIL_SPLINE 1
3#include <cmath>
4#include <vector>
5
6namespace lsst {
7namespace afw {
8namespace math {
9namespace detail {
10
11/*
12 * Splines
13 */
14class Spline {
15public:
16 virtual ~Spline() = default;
17
18 Spline(Spline const&) = default;
19 Spline(Spline&&) = default;
20 Spline& operator=(Spline const&) = default;
21 Spline& operator=(Spline&&) = default;
22
36 void derivative(std::vector<double> const& x, std::vector<double>& dydx) const;
37
46 std::vector<double> roots(double const value, double const x0, double const x1) const;
47
48protected:
49 Spline() = default;
53 void _allocateSpline(int const nknot);
54
55 std::vector<double> _knots; // positions of knots
56 std::vector<std::vector<double> > _coeffs; // and associated coefficients
57};
58
59class TautSpline : public Spline {
60public:
62
133 TautSpline(std::vector<double> const& x, std::vector<double> const& y, double const gamma = 0,
134 Symmetry type = Unknown);
135
136private:
144 void calculateTautSpline(std::vector<double> const& x, std::vector<double> const& y, double const gamma0);
155 void calculateTautSplineEvenOdd(std::vector<double> const& x, std::vector<double> const& y,
156 double const gamma0, bool even);
157};
158
159class SmoothedSpline : public Spline {
160public:
229 double s, double* chisq = nullptr, std::vector<double>* errs = nullptr);
230};
231} // namespace detail
232} // namespace math
233} // namespace afw
234} // namespace lsst
235#endif
SmoothedSpline(std::vector< double > const &x, std::vector< double > const &y, std::vector< double > const &dy, double s, double *chisq=nullptr, std::vector< double > *errs=nullptr)
Cubic spline data smoother.
Definition Spline.cc:534
std::vector< double > _knots
Definition Spline.h:55
Spline & operator=(Spline &&)=default
std::vector< double > roots(double const value, double const x0, double const x1) const
Find the roots of Spline - val = 0 in the range [x0, x1).
Definition Spline.cc:1226
Spline & operator=(Spline const &)=default
void derivative(std::vector< double > const &x, std::vector< double > &dydx) const
Find the derivative of a Spline.
Definition Spline.cc:57
Spline(Spline const &)=default
void _allocateSpline(int const nknot)
Allocate the storage a Spline needs.
Definition Spline.cc:21
void interpolate(std::vector< double > const &x, std::vector< double > &y) const
Interpolate a Spline.
Definition Spline.cc:29
std::vector< std::vector< double > > _coeffs
Definition Spline.h:56
TautSpline(std::vector< double > const &x, std::vector< double > const &y, double const gamma=0, Symmetry type=Unknown)
Construct cubic spline interpolant to given data.
Definition Spline.cc:86