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
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
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::vector< std::array< double, lsst::gauss2d::N_PARAMS_GAUSS2D > > GradParamFactors
std::vector< ParamBaseRef > ParamRefs
Definition param_defs.h:13
std::vector< std::array< size_t, lsst::gauss2d::N_PARAMS_GAUSS2D > > GradParamMap
std::vector< ExtraParamFactorValues > ExtraParamFactors
std::vector< std::array< size_t, lsst::gauss2d::N_EXTRA_MAP > > ExtraParamMap
std::vector< ParamBaseCRef > ParamCRefs
Definition param_defs.h:11
std::vector< std::shared_ptr< Component > > Components
std::map< ParamBaseCRef, size_t > ParameterMap
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 reserve(T... args)
T size(T... args)
Options for filtering Parameter instances.
T to_string(T... args)