gslsersicmixinterpolator.h

#ifndef LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H
#define LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H
 
#ifdef LSST_GAUSS2D_FIT_HAS_GSL
 
#include <memory>
 
#include "gsl.h"
#include "gslinterpolator.h"
#include "sersicmix.h"
 
#include <gsl/gsl_errno.h>
#include <gsl/gsl_spline.h>
 
namespace lsst::gauss2d::fit {

class GSLSersicMixInterpolator : public SersicMixInterpolator {
public:
    explicit GSLSersicMixInterpolator(unsigned short order = SERSICMIX_ORDER_DEFAULT,
                                      InterpType interp_type = INTERPTYPE_DEFAULT);
    ~GSLSersicMixInterpolator();
 
    bool correct_final_integral = true;
    static constexpr InterpType INTERPTYPE_DEFAULT = InterpType::cspline;
    const std::vector<SersicMixValues>& get_knots() const;

    double get_final_correction() const;
 
    std::vector<IntegralSize> get_integralsizes(double sersicindex) const override;
    std::vector<IntegralSize> get_integralsizes_derivs(double sersicindex) const override;
 
    InterpType get_interptype() const override;
    unsigned short get_order() const override;
 
    double get_sersicindex_min() const;
    double get_sersicindex_max() const;
 
    std::string repr(bool name_keywords = false,
                     std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
    std::string str() const override;
 
private:
    // This is allows caching of the final_correction value between calls to
    // get_integralsizes(_derivs); so get_final_correction can return it
    mutable double _final_correction = 1.;
    const InterpType _interp_type;
    const unsigned short _order;
    std::vector<std::pair<std::unique_ptr<GSLInterpolator>, std::unique_ptr<GSLInterpolator>>> _interps;
    const std::vector<SersicMixValues>& _knots;
    const double _sersicindex_min;
    const double _sersicindex_max;
};
 
}  // namespace lsst::gauss2d::fit
 
#endif  // LSST_GAUSS2D_FIT_HAS_GSL
#endif  // LSST_GAUSS2D_FIT_GSLSERSICMIXINTERPOLATOR_H