LSST Applications g013ef56533+b8d55c8942,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+7a3b874d60,g210f2d0738+7416ca6900,g262e1987ae+1d557ba9a3,g29ae962dfc+519d34895e,g2cef7863aa+aef1011c0b,g30d7c61c20+36d16ea71a,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+cb326ad149,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5a60e81ecd+6240c63dbc,g64539dfbff+7416ca6900,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g74acd417e5+0bae3c876a,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+dee7680868,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+14b8164b5b,g98a1a72a9c+8389601a76,g98df359435+fff771c62d,gb8cb2b794d+6728931916,gbf99507273+8c5ae1fdc5,gc2a301910b+7416ca6900,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+0bae3c876a,ge410e46f29+f459a6810c,ge41e95a9f2+7416ca6900,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
fractionalintegralmodel.h
Go to the documentation of this file.
1#ifndef LSST_GAUSS2D_FIT_FRACTIONALINTEGRALMODEL_H
2#define LSST_GAUSS2D_FIT_FRACTIONALINTEGRALMODEL_H
3
4#include <map>
5#include <memory>
6#include <optional>
7#include <vector>
8
9#include "component.h"
10#include "param_defs.h"
11#include "param_filter.h"
12#include "parameters.h"
13#include "integralmodel.h"
14
15namespace lsst::gauss2d::fit {
31public:
35
38
40
43
44 typename Data::iterator begin() noexcept;
45 typename Data::const_iterator cbegin() const noexcept;
46
47 typename Data::iterator end() noexcept;
48 typename Data::const_iterator cend() const noexcept;
49
56 static std::shared_ptr<FractionalIntegralModel> find_model(const IntegralModel& model) {
57 const auto found = _registry_rev.find(model);
58 return (found == _registry_rev.end()) ? nullptr : (*found).second.lock();
59 }
60
62 const IntegralModel& get_parent_model() const;
63 double get_integral(const Channel& channel) const override;
65 const Channel& channel) const override;
66 double get_integral_remainder(const Channel& channel) const;
67
69 ParamRefs& get_parameters(ParamRefs& params, ParamFilter* filter = nullptr) const override;
70 ParamCRefs& get_parameters_const(ParamCRefs& params, ParamFilter* filter = nullptr) const override;
71
72 bool is_final() const;
73
86 static std::shared_ptr<FractionalIntegralModel> make(std::optional<const Data> data,
88 bool is_final = false);
90 std::optional<const Data> data, std::shared_ptr<const IntegralModel> model,
91 bool is_final = false);
92
93 size_t size() const;
94
95 std::string repr(bool name_keywords = false,
96 std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR) const override;
97 std::string str() const override;
98
99private:
100 Data _data = {};
101 // This could be unordered, but std::hash<std::string> won't take const strings
102 // (see also linearintegralmodel.h)
104
105 // TODO: See if all raw pointers can be changed to reference_wrappers or weak_ptrs
111 _registry = {};
114 _registry_rev = {};
115
116 struct Shared_enabler;
117
118 bool _is_final;
119
120 // not giving a nullptr default data_in because the map needs to match the model's channels
121 FractionalIntegralModel(std::optional<const Data> data, std::shared_ptr<const IntegralModel> model,
122 bool is_final);
123};
124
125} // namespace lsst::gauss2d::fit
126
127#endif
static constexpr std::string_view CC_NAMESPACE_SEPARATOR
The C++ namespace separator.
Definition object.h:45
An observational channel, usually representing some range of wavelengths of light.
Definition channel.h:29
A list of Observation instances that can be modelled.
Definition data.h:32
static std::shared_ptr< FractionalIntegralModel > find_model(const IntegralModel &model)
Find the FractionalIntegralModel that depends on a given IntegralModel, if any.
ParamCRefs & get_parameters_const(ParamCRefs &params, ParamFilter *filter=nullptr) const override
Same as get_parameters(), but for const refs.
std::string repr(bool name_keywords=false, std::string_view namespace_separator=Object::CC_NAMESPACE_SEPARATOR) const override
Return a full, callable string representation of this.
Data::const_iterator cend() const noexcept
ProperFractionParameterD & get_parameter_frac(const Channel &channel) const
FractionalIntegralModel & operator=(const FractionalIntegralModel &)=delete
std::pair< std::reference_wrapper< const Channel >, std::shared_ptr< ProperFractionParameterD > > ChannelIntegralParameterD
std::string str() const override
Return a brief, human-readable string representation of this.
double get_integral_remainder(const Channel &channel) const
std::vector< std::pair< ParamBaseCRef, ExtraParamFactorValues > > get_integral_derivative_factors(const Channel &channel) const override
Return the partial derivative of the model w.r.t.
ParamRefs & get_parameters(ParamRefs &params, ParamFilter *filter=nullptr) const override
Add Parameter refs matching the filter to a vector, in order.
std::vector< std::reference_wrapper< const Channel > > get_channels() const override
Get the set of channels this instance is applicable for.
double get_integral(const Channel &channel) const override
Get the value of the integral in a single Channel.
std::vector< ChannelIntegralParameterD > Data
static const std::shared_ptr< const FractionalIntegralModel > make_const(std::optional< const Data > data, std::shared_ptr< const IntegralModel > model, bool is_final=false)
Data::const_iterator cbegin() const noexcept
std::shared_ptr< ProperFractionParameterD > at(const Channel &channel)
FractionalIntegralModel(const FractionalIntegralModel &)=delete
static std::shared_ptr< FractionalIntegralModel > make(std::optional< const Data > data, std::shared_ptr< const IntegralModel > model, bool is_final=false)
Construct a FractionalIntegralModel and add to registry.
A Parametric model for the integral of a 2D distribution.
std::vector< ParamBaseRef > ParamRefs
Definition param_defs.h:13
std::vector< ParamBaseCRef > ParamCRefs
Definition param_defs.h:11
STL namespace.
Options for filtering Parameter instances.