LSST Applications g013ef56533+7c9321ec0f,g042eb84c57+c6cfa41bc3,g199a45376c+0ba108daf9,g1fd858c14a+fcad0d0313,g210f2d0738+c0f94c6586,g262e1987ae+a7e710680e,g29ae962dfc+fb55f2edb0,g2ac17093b6+61d6563b1e,g2b1d02342f+df6f932764,g2cef7863aa+aef1011c0b,g2f7ad74990+c0f94c6586,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+53cf87ae69,g47891489e3+4316d04fff,g511e8cfd20+baa56acf6c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+fd7ad03fde,g64539dfbff+c0f94c6586,g67b6fd64d1+4316d04fff,g67fd3c3899+c0f94c6586,g6985122a63+4316d04fff,g74acd417e5+ca833bee28,g786e29fd12+668abc6043,g81db2e9a8d+b2ec8e584f,g87389fa792+8856018cbb,g89139ef638+4316d04fff,g8d7436a09f+0a24083b20,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+11eb8fd5b8,gbf99507273+8c5ae1fdc5,gcdda8b9158+e4c84c9d5c,gce8aa8abaa+8c5ae1fdc5,gd7ef33dd92+4316d04fff,gdab6d2f7ff+ca833bee28,ge410e46f29+4316d04fff,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.40
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.