sersicmix.h

#ifndef LSST_GAUSS2D_FIT_SERSICMIX_H
#define LSST_GAUSS2D_FIT_SERSICMIX_H
 
#include <array>
#include <memory>
#include <stdexcept>
#include <vector>
 
#include "interpolation.h"
#include "lsst/gauss2d/object.h"
 
namespace lsst::gauss2d::fit {
 
static const unsigned short SERSICMIX_ORDER_DEFAULT = 4;

class IntegralSize : public Object {
public:
    explicit IntegralSize(double integral_, double sigma_);
 
    const double integral;
    const double sigma;
 
    std::string repr(bool name_keywords = false,
                     std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
    std::string str() const override;
};

class SersicMixInterpolator : public Object {
public:
    virtual std::vector<IntegralSize> get_integralsizes(double sersicindex) const = 0;
    virtual std::vector<IntegralSize> get_integralsizes_derivs(double sersicindex) const = 0;
 
    virtual InterpType get_interptype() const = 0;
    virtual unsigned short get_order() const = 0;
};
 
const std::shared_ptr<const SersicMixInterpolator> get_sersic_mix_interpolator_default(
        unsigned short order = SERSICMIX_ORDER_DEFAULT);

class SersicMixValues : public Object {
public:
    explicit SersicMixValues(double sersicindex, std::vector<IntegralSize> values);
 
    const double sersicindex;
    const std::vector<IntegralSize> values;
 
    std::string repr(bool name_keywords = false,
                     std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
    std::string str() const override;
};
 
inline bool operator<(const SersicMixValues &lhs, const SersicMixValues &rhs) {
    return lhs.sersicindex < rhs.sersicindex;
}
inline bool operator<(const SersicMixValues &lhs, double x) { return lhs.sersicindex < x; }
inline bool operator<(double x, const SersicMixValues &rhs) { return x < rhs.sersicindex; }
 
std::vector<SersicMixValues> get_sersic_mix_knots_copy(unsigned short order);
 
const std::vector<SersicMixValues> &get_sersic_mix_knots(unsigned short order);
 
const std::vector<SersicMixValues> &get_sersic_mix_knots_order4();
const std::vector<SersicMixValues> &get_sersic_mix_knots_order8();
 
}  // namespace lsst::gauss2d::fit
 
#endif