LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
Quadrupole.cc
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  */
25 
26 namespace lsst {
27 namespace afw {
28 namespace geom {
29 namespace ellipses {
30 
31 BaseCore::Registrar<Quadrupole> Quadrupole::registrar;
32 
33 std::string Quadrupole::getName() const { return "Quadrupole"; }
34 
36  if (_matrix(0, 1) != _matrix(1, 0))
38  "Quadrupole matrix must be symmetric.");
39  if (getIxx() < 0 || getIyy() < 0)
41  "Quadrupole matrix cannot have negative diagonal elements.");
42  if (getDeterminant() < 0)
44  "Quadrupole matrix cannot have negative determinant.");
45 }
46 
47 void Quadrupole::readParameters(double const* iter) {
48  setIxx(*iter++);
49  setIyy(*iter++);
50  setIxy(*iter++);
51 }
52 
53 void Quadrupole::writeParameters(double* iter) const {
54  *iter++ = getIxx();
55  *iter++ = getIyy();
56  *iter++ = getIxy();
57 }
58 
59 Quadrupole::Quadrupole(double ixx, double iyy, double ixy, bool normalize) {
60  setIxx(ixx);
61  setIyy(iyy);
62  setIxy(ixy);
63  if (normalize) this->normalize();
64 }
65 
66 Quadrupole::Quadrupole(BaseCore::ParameterVector const& vector, bool normalize) {
67  setIxx(vector[IXX]);
68  setIyy(vector[IYY]);
69  setIxy(vector[IXY]);
70  if (normalize) this->normalize();
71 }
72 
73 Quadrupole::Quadrupole(Matrix const& matrix, bool normalize) : _matrix(matrix) {
74  if (normalize) this->normalize();
75 }
76 
77 void Quadrupole::_assignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
78  ixx = getIxx();
79  iyy = getIyy();
80  ixy = getIxy();
81 }
82 
83 BaseCore::Jacobian Quadrupole::_dAssignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
84  ixx = getIxx();
85  iyy = getIyy();
86  ixy = getIxy();
87  return Jacobian::Identity();
88 }
89 
90 void Quadrupole::_assignToAxes(double& a, double& b, double& theta) const {
92 }
93 
94 BaseCore::Jacobian Quadrupole::_dAssignToAxes(double& a, double& b, double& theta) const {
95  return BaseCore::_dAssignQuadrupoleToAxes(getIxx(), getIyy(), getIxy(), a, b, theta);
96 }
97 
98 void Quadrupole::_assignFromQuadrupole(double ixx, double iyy, double ixy) {
99  setIxx(ixx);
100  setIyy(iyy);
101  setIxy(ixy);
102 }
103 
104 BaseCore::Jacobian Quadrupole::_dAssignFromQuadrupole(double ixx, double iyy, double ixy) {
105  setIxx(ixx);
106  setIyy(iyy);
107  setIxy(ixy);
108  return Jacobian::Identity();
109 }
110 
111 void Quadrupole::_assignFromAxes(double a, double b, double theta) {
112  BaseCore::_assignAxesToQuadrupole(a, b, theta, _matrix(0, 0), _matrix(1, 1), _matrix(0, 1));
113  _matrix(1, 0) = _matrix(0, 1);
114 }
115 
116 BaseCore::Jacobian Quadrupole::_dAssignFromAxes(double a, double b, double theta) {
117  Jacobian r = BaseCore::_dAssignAxesToQuadrupole(a, b, theta, _matrix(0, 0), _matrix(1, 1), _matrix(0, 1));
118  _matrix(1, 0) = _matrix(0, 1);
119  return r;
120 }
121 } // namespace ellipses
122 } // namespace geom
123 } // namespace afw
124 } // namespace lsst
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
table::Key< int > b
table::Key< int > a
static Jacobian _dAssignAxesToQuadrupole(double a, double b, double theta, double &ixx, double &iyy, double &ixy)
Return the size of the bounding box for the ellipse core.
Definition: BaseCore.cc:231
Eigen::Vector3d ParameterVector
Parameter vector type.
Definition: BaseCore.h:63
static void _assignAxesToQuadrupole(double a, double b, double theta, double &ixx, double &iyy, double &ixy)
Return the size of the bounding box for the ellipse core.
Definition: BaseCore.cc:218
static Jacobian _dAssignQuadrupoleToAxes(double ixx, double iyy, double ixy, double &a, double &b, double &theta)
Return the size of the bounding box for the ellipse core.
Definition: BaseCore.cc:193
static void _assignQuadrupoleToAxes(double ixx, double iyy, double ixy, double &a, double &b, double &theta)
Return the size of the bounding box for the ellipse core.
Definition: BaseCore.cc:183
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
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
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
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
Jacobian _dAssignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:116
void _assignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:111
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.
Definition: Quadrupole.cc:104
Reports invalid arguments.
Definition: Runtime.h:66
A base class for image defects.