LSST Applications g044012fb7c+6976b5ec80,g04a91732dc+88a5fc122b,g07dc498a13+7e3c5f68a2,g114c6a66ad+09472d7a76,g1409bbee79+7e3c5f68a2,g1a7e361dbc+7e3c5f68a2,g1fd858c14a+3a43eabc0e,g35bb328faa+fcb1d3bbc8,g3bd4b5ce2c+2647bb081c,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g5477a8d5ce+b19c77c7ae,g58d0cdf3ff+4a2e102ff8,g60b5630c4e+09472d7a76,g623d845a50+09472d7a76,g6f0c2978f1+fcf1c0bcd6,g71fabbc107+09472d7a76,g75b6c65c88+d0b1dc44cc,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+349c7e81d4,g89139ef638+7e3c5f68a2,g9125e01d80+fcb1d3bbc8,g95236ca021+f7a31438ed,g989de1cb63+7e3c5f68a2,g9f33ca652e+f17d666fbc,gaaedd4e678+7e3c5f68a2,gabe3b4be73+1e0a283bba,gb1101e3267+f870f33517,gb58c049af0+f03b321e39,gc99c83e5f0+76d20ab76d,gcf25f946ba+349c7e81d4,gd0fa69b896+f3a65fa83c,gd6cbbdb0b4+c8606af20c,gde0f65d7ad+5bd27d919f,ge278dab8ac+932305ba37,gfba249425e+fcb1d3bbc8,w.2025.07
LSST Data Management Base Package
Loading...
Searching...
No Matches
Quadrupole.h
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008, 2009, 2010 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24
25#ifndef LSST_AFW_GEOM_ELLIPSES_Quadrupole_h_INCLUDED
26#define LSST_AFW_GEOM_ELLIPSES_Quadrupole_h_INCLUDED
27
28/*
29 * Definitions and inlines for Quadrupole.
30 *
31 * Note: do not include directly; use the main ellipse header file.
32 */
33
38
39namespace lsst {
40namespace afw {
41namespace geom {
42namespace ellipses {
43
47class Quadrupole : public BaseCore {
48public:
49 enum ParameterEnum { IXX = 0, IYY = 1, IXY = 2 };
50
52 using Matrix = Eigen::Matrix<double, 2, 2, Eigen::DontAlign>;
53
54 double const getIxx() const { return _matrix(0, 0); }
55 void setIxx(double ixx) { _matrix(0, 0) = ixx; }
56
57 double const getIyy() const { return _matrix(1, 1); }
58 void setIyy(double iyy) { _matrix(1, 1) = iyy; }
59
60 double const getIxy() const { return _matrix(1, 0); }
61 void setIxy(double ixy) { _matrix(0, 1) = _matrix(1, 0) = ixy; }
62
65
67 std::string getName() const override;
68
73 void normalize() override;
74
75 void readParameters(double const* iter) override;
76
77 void writeParameters(double* iter) const override;
78
80 Matrix const& getMatrix() const { return _matrix; }
81
83 double getDeterminant() const { return getIxx() * getIyy() - getIxy() * getIxy(); }
84
87 _matrix = other._matrix;
88 return *this;
89 }
90
91 // Delegate to copy-assignment for backwards compatibility
92 Quadrupole& operator=(Quadrupole&& other) { return *this = other; }
93
95 Quadrupole& operator=(BaseCore const& other) {
97 return *this;
98 }
99
101 explicit Quadrupole(double ixx = 1.0, double iyy = 1.0, double ixy = 0.0, bool normalize = false);
102
104 explicit Quadrupole(BaseCore::ParameterVector const& vector, bool normalize = false);
105
107 explicit Quadrupole(Matrix const& matrix, bool normalize = true);
108
110 Quadrupole(Quadrupole const& other) : _matrix(other._matrix) {}
111
112 // Delegate to copy-constructor for backwards compatibility
113 Quadrupole(Quadrupole&& other) : Quadrupole(other) {}
114
115 ~Quadrupole() override = default;
116
118 Quadrupole(BaseCore const& other) { *this = other; }
119
121 Quadrupole(BaseCore::Transformer const& transformer) { transformer.apply(*this); }
122
124 Quadrupole(BaseCore::Convolution const& convolution) { convolution.apply(*this); }
125
126protected:
128
129 void _assignToQuadrupole(double& ixx, double& iyy, double& ixy) const override;
130 void _assignFromQuadrupole(double ixx, double iyy, double ixy) override;
131
132 void _assignToAxes(double& a, double& b, double& theta) const override;
133 void _assignFromAxes(double a, double b, double theta) override;
134
135 Jacobian _dAssignToQuadrupole(double& ixx, double& iyy, double& ixy) const override;
136 Jacobian _dAssignFromQuadrupole(double ixx, double iyy, double ixy) override;
137
138 Jacobian _dAssignToAxes(double& a, double& b, double& theta) const override;
139 Jacobian _dAssignFromAxes(double a, double b, double theta) override;
140
141private:
142 static Registrar<Quadrupole> registrar;
143
144 Matrix _matrix;
145};
146} // namespace ellipses
147} // namespace geom
148} // namespace afw
149} // namespace lsst
150
151#endif // !LSST_AFW_GEOM_ELLIPSES_Quadrupole_h_INCLUDED
A temporary-only expression object for ellipse core convolution.
Definition Convolution.h:44
A temporary-only expression object for ellipse core transformations.
Definition Transformer.h:49
Eigen::Vector3d ParameterVector
Parameter vector type.
Definition BaseCore.h:63
BaseCore()=default
Return the size of the bounding box for the ellipse core.
BaseCore & operator=(BaseCore const &other)
Set the parameters of this ellipse core from another.
Definition BaseCore.cc:155
Eigen::Matrix3d Jacobian
Parameter Jacobian matrix type.
Definition BaseCore.h:64
Eigen::Matrix< double, 2, 2, Eigen::DontAlign > Matrix
Matrix type for the matrix representation of Quadrupole parameters.
Definition Quadrupole.h:52
Quadrupole & operator=(Quadrupole const &other)
Standard assignment.
Definition Quadrupole.h:86
void _assignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:90
double getDeterminant() const
Return the determinant of the matrix representation.
Definition Quadrupole.h:83
void normalize() override
Put the parameters into a "standard form", and throw InvalidParameterError if they cannot be normaliz...
Definition Quadrupole.cc:35
Matrix const & getMatrix() const
Return a 2x2 symmetric matrix of the parameters.
Definition Quadrupole.h:80
Quadrupole(Quadrupole const &other)
Copy constructor.
Definition Quadrupole.h:110
Quadrupole(BaseCore::Convolution const &convolution)
Converting copy constructor.
Definition Quadrupole.h:124
std::shared_ptr< Quadrupole > clone() const
Deep copy the ellipse core.
Definition Quadrupole.h:64
Jacobian _dAssignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:83
Quadrupole(double ixx=1.0, double iyy=1.0, double ixy=0.0, bool normalize=false)
Construct from parameter values.
Definition Quadrupole.cc:59
void _assignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:98
std::string getName() const override
Return a string that identifies this parametrization.
Definition Quadrupole.cc:33
Quadrupole & operator=(BaseCore const &other)
Converting assignment.
Definition Quadrupole.h:95
void writeParameters(double *iter) const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:53
void _assignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:77
Jacobian _dAssignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:94
Quadrupole(BaseCore const &other)
Converting copy constructor.
Definition Quadrupole.h:118
Quadrupole(BaseCore::Transformer const &transformer)
Converting copy constructor.
Definition Quadrupole.h:121
Jacobian _dAssignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
void _assignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
void readParameters(double const *iter) override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.cc:47
Jacobian _dAssignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Quadrupole & operator=(Quadrupole &&other)
Definition Quadrupole.h:92
std::shared_ptr< BaseCore > _clone() const override
Return the size of the bounding box for the ellipse core.
Definition Quadrupole.h:127
T make_shared(T... args)
BaseCore::Registrar< Separable< Ellipticity_, Radius_ > > Separable< Ellipticity_, Radius_ >::registrar
Definition Separable.cc:36
T static_pointer_cast(T... args)