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
linearintegralmodel.cc
Go to the documentation of this file.
2
6#include "lsst/gauss2d/fit/parameters.h"
8
9namespace lsst::gauss2d::fit {
10
11// not giving a nullptr default because users should explicitly use the null Channel
13 if (data_in != nullptr) {
14 size_t idx = 0;
15 for (const auto& datum : *data_in) {
16 if (datum.second == nullptr) {
17 throw std::runtime_error("LinearIntegralModel data[" + std::to_string(idx)
18 + "] can't be null");
19 }
20 if (_map.find(datum.first) != _map.end()) {
21 throw std::runtime_error("LinearIntegralModel data[" + std::to_string(idx)
22 + "] channel=" + datum.first.get().str() + " duplicated");
23 }
24 _data.emplace_back(datum.first, datum.second);
25 _map.insert(_data.back());
26 idx++;
27 }
28 } else {
30 _map.insert(_data.back());
31 }
32}
34
37 for (auto& datum : _data) rval.emplace_back(datum.first);
38 return rval;
39}
40
41double LinearIntegralModel::get_integral(const Channel& channel) const {
42 return _map.at(channel)->get_value();
43}
44
47 return {};
48}
49
51 return _map.at(channel);
52}
53
55 return _map.at(channel);
56}
57
58typename LinearIntegralModel::Data::iterator LinearIntegralModel::begin() noexcept { return _data.begin(); }
59typename LinearIntegralModel::Data::const_iterator LinearIntegralModel::cbegin() const noexcept {
60 return _data.begin();
61}
62
63typename LinearIntegralModel::Data::iterator LinearIntegralModel::end() noexcept { return _data.end(); }
64typename LinearIntegralModel::Data::const_iterator LinearIntegralModel::cend() const noexcept {
65 return _data.cend();
66}
67
69 for (const auto& p : _data) insert_param_channel(p.first, *p.second, params, filter);
70 return params;
71}
73 for (const auto& p : _data) insert_param_channel(p.first, *p.second, params, filter);
74 return params;
75}
76
77size_t LinearIntegralModel::size() const { return _data.size(); }
78
79std::string LinearIntegralModel::repr(bool name_keywords, std::string_view namespace_separator) const {
80 std::string s = type_name_str<LinearIntegralModel>(false, namespace_separator) + "("
81 + (name_keywords ? "data=" : "") + "{";
82 for (const auto& datum : _data) {
83 s += datum.first.get().repr(name_keywords, namespace_separator) + ": "
84 + datum.second->repr(name_keywords, namespace_separator) + ",";
85 }
86 return s + "})";
87}
88
91 for (const auto& datum : _data) {
92 s += datum.first.get().str() + ": " + datum.second->str() + ",";
93 }
94 return s + "})";
95}
96
97} // namespace lsst::gauss2d::fit
T at(T... args)
An observational channel, usually representing some range of wavelengths of light.
Definition channel.h:29
static const Channel & NONE()
Definition channel.cc:111
LinearIntegralModel(const Data *data_in)
Construct a LinearIntegralModel from input Data.
Data::const_iterator cbegin() const noexcept
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.
std::vector< std::reference_wrapper< const Channel > > get_channels() const override
Get the set of channels this instance is applicable for.
std::shared_ptr< IntegralParameterD > at(const Channel &channel)
Get the IntegralParameterD for the given Channel.
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.
double get_integral(const Channel &channel) const override
Get the value of the integral in a single Channel.
ParamRefs & get_parameters(ParamRefs &params, ParamFilter *filter=nullptr) const override
Add Parameter refs matching the filter to a vector, in order.
std::string str() const override
Return a brief, human-readable string representation of this.
std::vector< ChannelIntegralParameterD > Data
Data::const_iterator cend() const noexcept
size_t size() const
Return the size of Data (number of Channel/IntegralParameterD instances)
T emplace_back(T... args)
T make_shared(T... args)
std::vector< ParamBaseRef > ParamRefs
Definition param_defs.h:13
std::vector< ParamBaseCRef > ParamCRefs
Definition param_defs.h:11
void insert_param_channel(const lsst::gauss2d::fit::Channel &channel, g2f::ParamBase &param, t &params, ParamFilter *filter)
Same as insert_param, but only inserting if filter has a null or matching Channel.
std::string type_name_str(bool strip_namespace=false, std::string_view namespace_str=detail::NAMESPACE_SEPARATOR)
Get a string representation of an arbitrary C++ type, potentially modifying its namespace prefix.
Definition type_name.h:104
T size(T... args)
Options for filtering Parameter instances.
T to_string(T... args)