lsst.afw g44befd6440+dd219f3bc1
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,
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
double x
table::Key< int > type
Definition: Detector.cc:163
int y
Definition: SpanSet.cc:48
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
A base class for image defects.