Loading [MathJax]/extensions/tex2jax.js
LSST Applications 28.0.0,g1653933729+a8ce1bb630,g1a997c3884+a8ce1bb630,g28da252d5a+5bd70b7e6d,g2bbee38e9b+638fca75ac,g2bc492864f+638fca75ac,g3156d2b45e+07302053f8,g347aa1857d+638fca75ac,g35bb328faa+a8ce1bb630,g3a166c0a6a+638fca75ac,g3e281a1b8c+7bbb0b2507,g4005a62e65+17cd334064,g414038480c+5b5cd4fff3,g41af890bb2+4ffae9de63,g4e1a3235cc+0f1912dca3,g6249c6f860+3c3976f90c,g80478fca09+46aba80bd6,g82479be7b0+77990446f6,g858d7b2824+78ba4d1ce1,g89c8672015+f667a5183b,g9125e01d80+a8ce1bb630,ga5288a1d22+2a6264e9ca,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc22bb204ba+78ba4d1ce1,gc28159a63d+638fca75ac,gcf0d15dbbd+32ddb6096f,gd6b7c0dfd1+3e339405e9,gda3e153d99+78ba4d1ce1,gda6a2b7d83+32ddb6096f,gdaeeff99f8+1711a396fd,gdd5a9049c5+b18c39e5e3,ge2409df99d+a5e4577cdc,ge33fd446bb+78ba4d1ce1,ge79ae78c31+638fca75ac,gf0baf85859+64e8883e75,gf5289d68f6+e1b046a8d7,gfa443fc69c+91d9ed1ecf,gfda6b12a05+8419469a56
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
gslsersicmixinterpolator.h
Go to the documentation of this file.
1#ifndef LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H
2#define LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H
3
4#ifdef LSST_GAUSS2D_FIT_HAS_GSL
5
6#include <memory>
7
8#include "gsl.h"
9#include "gslinterpolator.h"
10#include "sersicmix.h"
11
12#include <gsl/gsl_errno.h>
13#include <gsl/gsl_spline.h>
14
15namespace lsst::gauss2d::fit {
16
20class GSLSersicMixInterpolator : public SersicMixInterpolator {
21public:
22 explicit GSLSersicMixInterpolator(unsigned short order = SERSICMIX_ORDER_DEFAULT,
23 InterpType interp_type = INTERPTYPE_DEFAULT);
24 ~GSLSersicMixInterpolator();
25
26 bool correct_final_integral = true;
27 static constexpr InterpType INTERPTYPE_DEFAULT = InterpType::cspline;
29 const std::vector<SersicMixValues>& get_knots() const;
30
33 double get_final_correction() const;
34
35 std::vector<IntegralSize> get_integralsizes(double sersicindex) const override;
36 std::vector<IntegralSize> get_integralsizes_derivs(double sersicindex) const override;
37
38 InterpType get_interptype() const override;
39 unsigned short get_order() const override;
40
41 double get_sersicindex_min() const;
42 double get_sersicindex_max() const;
43
44 std::string repr(bool name_keywords = false,
45 std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
46 std::string str() const override;
47
48private:
49 // This is allows caching of the final_correction value between calls to
50 // get_integralsizes(_derivs); so get_final_correction can return it
51 mutable double _final_correction = 1.;
52 const InterpType _interp_type;
53 const unsigned short _order;
55 const std::vector<SersicMixValues>& _knots;
56 const double _sersicindex_min;
57 const double _sersicindex_max;
58};
59
60} // namespace lsst::gauss2d::fit
61
62#endif // LSST_GAUSS2D_FIT_HAS_GSL
63#endif // LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H
table::Key< int > order