Loading [MathJax]/extensions/tex2jax.js
LSST Applications 28.0.0,g1653933729+a8ce1bb630,g1a997c3884+a8ce1bb630,g28da252d5a+5bd70b7e6d,g2bbee38e9b+638fca75ac,g2bc492864f+638fca75ac,g3156d2b45e+07302053f8,g347aa1857d+638fca75ac,g35bb328faa+a8ce1bb630,g3a166c0a6a+638fca75ac,g3e281a1b8c+7bbb0b2507,g4005a62e65+17cd334064,g414038480c+5b5cd4fff3,g41af890bb2+4ffae9de63,g4e1a3235cc+0f1912dca3,g6249c6f860+3c3976f90c,g80478fca09+46aba80bd6,g82479be7b0+77990446f6,g858d7b2824+78ba4d1ce1,g89c8672015+f667a5183b,g9125e01d80+a8ce1bb630,ga5288a1d22+2a6264e9ca,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc22bb204ba+78ba4d1ce1,gc28159a63d+638fca75ac,gcf0d15dbbd+32ddb6096f,gd6b7c0dfd1+3e339405e9,gda3e153d99+78ba4d1ce1,gda6a2b7d83+32ddb6096f,gdaeeff99f8+1711a396fd,gdd5a9049c5+b18c39e5e3,ge2409df99d+a5e4577cdc,ge33fd446bb+78ba4d1ce1,ge79ae78c31+638fca75ac,gf0baf85859+64e8883e75,gf5289d68f6+e1b046a8d7,gfa443fc69c+91d9ed1ecf,gfda6b12a05+8419469a56
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
psfmodel.cc
Go to the documentation of this file.
1#include <memory>
2#include <optional>
3
5
9
10namespace lsst::gauss2d::fit {
11
13 size_t i = 0;
14 _components.reserve(components.size());
15 for (auto& component : components) {
16 if (component == nullptr)
17 throw std::invalid_argument("PsfModel components[" + std::to_string(i) + "] can't be null");
18 auto channels = component->get_integralmodel().get_channels();
19 if ((channels.size() != 1) || ((*channels.begin()).get() != Channel::NONE())) {
20 throw std::invalid_argument("PsfModel components[" + std::to_string(i)
21 + "].get_integralmodel().get_channels()="
22 + str_iter_ref<true>(channels) + " must only contain None");
23 }
24 _components.push_back(std::move(component));
25 i++;
26 }
27}
28
30
31void PsfModel::add_extra_param_map(const Channel& channel, ExtraParamMap& map_extra,
32 const GradParamMap& map_grad, ParameterMap& offsets) const {
33 for (auto& component : _components) component->add_extra_param_map(channel, map_extra, map_grad, offsets);
34}
35
36void PsfModel::add_extra_param_factors(const Channel& channel, ExtraParamFactors& factors) const {
37 for (auto& component : _components) component->add_extra_param_factors(channel, factors);
38}
39
40void PsfModel::add_grad_param_map(const Channel& channel, GradParamMap& map, ParameterMap& offsets) const {
41 for (auto& component : _components) component->add_grad_param_map(channel, map, offsets);
42}
43
44void PsfModel::add_grad_param_factors(const Channel& channel, GradParamFactors& factors) const {
45 for (auto& component : _components) component->add_grad_param_factors(channel, factors);
46}
47
48Components PsfModel::get_components() const { return _components; }
49
52 in.reserve(_components.size());
53 for (auto& component : _components) {
54 in.push_back(component->get_gaussians(channel)->get_data());
55 }
56 return std::make_unique<lsst::gauss2d::Gaussians>(in);
57}
58
59size_t PsfModel::get_n_gaussians(const Channel& channel) const {
60 size_t n_g = 0;
61 for (auto& component : _components) n_g += component->get_n_gaussians(channel);
62 return n_g;
63}
64
66 for (auto& component : _components) component->get_parameters(params, filter);
67 return params;
68}
70 for (auto& component : _components) component->get_parameters_const(params, filter);
71 return params;
72}
73
75 size_t index) const {
76 for (auto& component : _components) component->set_extra_param_factors(channel, factors, index);
77}
78
79void PsfModel::set_grad_param_factors(const Channel& channel, GradParamFactors& factors, size_t index) const {
80 for (auto& component : _components) component->set_grad_param_factors(channel, factors, index);
81}
82
83std::string PsfModel::repr(bool name_keywords, std::string_view namespace_separator) const {
84 std::string str = type_name_str<PsfModel>(false, namespace_separator) + "("
85 + (name_keywords ? "components=[" : "[");
86 for (const auto& s : _components) str += s->repr(name_keywords, namespace_separator) + ",";
87 return str + "])";
88}
89
91 std::string str = type_name_str<PsfModel>(true) + "(components=[";
92 for (const auto& s : _components) str += s->str() + ",";
93 return str + "])";
94}
95
96} // namespace lsst::gauss2d::fit
table::Key< table::Array< int > > components
An observational channel, usually representing some range of wavelengths of light.
Definition channel.h:29
static const Channel & NONE()
Definition channel.cc:111
void set_grad_param_factors(const Channel &channel, GradParamFactors &factor, size_t index) const override
Set Parameter gradient factors in an existing map.
Definition psfmodel.cc:79
void add_grad_param_map(const Channel &channel, GradParamMap &map, ParameterMap &offsets) const override
Add Parameter gradient indices to an existing map.
Definition psfmodel.cc:40
void add_extra_param_map(const Channel &channel, ExtraParamMap &map_extra, const GradParamMap &map_grad, ParameterMap &offsets) const override
Add extra Parameter indices to a map.
Definition psfmodel.cc:31
void add_extra_param_factors(const Channel &channel, ExtraParamFactors &factors) const override
Add extra Parameter gradient factors to an existing vector.
Definition psfmodel.cc:36
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.
Definition psfmodel.cc:83
void add_grad_param_factors(const Channel &channel, GradParamFactors &factor) const override
Add Parameter gradient factors to an existing map.
Definition psfmodel.cc:44
void set_extra_param_factors(const Channel &channel, ExtraParamFactors &factors, size_t index) const override
Set extra Parameter gradient factors in an existing map.
Definition psfmodel.cc:74
std::unique_ptr< const lsst::gauss2d::Gaussians > get_gaussians(const Channel &channel=Channel::NONE()) const override
Return the vector of Gaussian sub-components controlled by this model.
Definition psfmodel.cc:50
PsfModel(Components &components)
Definition psfmodel.cc:12
std::string str() const override
Return a brief, human-readable string representation of this.
Definition psfmodel.cc:90
size_t get_n_gaussians(const Channel &channel=Channel::NONE()) const override
Return the number of Gaussian sub-components controlled by this model.
Definition psfmodel.cc:59
ParamCRefs & get_parameters_const(ParamCRefs &params, ParamFilter *filter=nullptr) const override
Same as get_parameters(), but for const refs.
Definition psfmodel.cc:69
Components get_components() const override
Definition psfmodel.cc:48
ParamRefs & get_parameters(ParamRefs &params, ParamFilter *filter=nullptr) const override
Add Parameter refs matching the filter to a vector, in order.
Definition psfmodel.cc:65
T move(T... args)
std::string str_iter_ref(const T &container)
Definition object.h:148
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 push_back(T... args)
T reserve(T... args)
T size(T... args)
Options for filtering Parameter instances.
T to_string(T... args)