LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
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  */
26 
27 namespace lsst {
28 namespace afw {
29 namespace geom {
30 namespace ellipses {
31 
32 BaseCore::Registrar<Quadrupole> Quadrupole::registrar;
33 
34 std::string Quadrupole::getName() const { return "Quadrupole"; }
35 
37  if (_matrix(0, 1) != _matrix(1, 0))
39  "Quadrupole matrix must be symmetric.");
40  if (getIxx() < 0 || getIyy() < 0)
42  "Quadrupole matrix cannot have negative diagonal elements.");
43  if (getDeterminant() < 0)
45  "Quadrupole matrix cannot have negative determinant.");
46 }
47 
48 void Quadrupole::readParameters(double const* iter) {
49  setIxx(*iter++);
50  setIyy(*iter++);
51  setIxy(*iter++);
52 }
53 
54 void Quadrupole::writeParameters(double* iter) const {
55  *iter++ = getIxx();
56  *iter++ = getIyy();
57  *iter++ = getIxy();
58 }
59 
60 Quadrupole::Quadrupole(double ixx, double iyy, double ixy, bool normalize) {
61  setIxx(ixx);
62  setIyy(iyy);
63  setIxy(ixy);
64  if (normalize) this->normalize();
65 }
66 
67 Quadrupole::Quadrupole(BaseCore::ParameterVector const& vector, bool normalize) {
68  setIxx(vector[IXX]);
69  setIyy(vector[IYY]);
70  setIxy(vector[IXY]);
71  if (normalize) this->normalize();
72 }
73 
74 Quadrupole::Quadrupole(Matrix const& matrix, bool normalize) : _matrix(matrix) {
75  if (normalize) this->normalize();
76 }
77 
78 void Quadrupole::_assignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
79  ixx = getIxx();
80  iyy = getIyy();
81  ixy = getIxy();
82 }
83 
84 BaseCore::Jacobian Quadrupole::_dAssignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
85  ixx = getIxx();
86  iyy = getIyy();
87  ixy = getIxy();
88  return Jacobian::Identity();
89 }
90 
91 void Quadrupole::_assignToAxes(double& a, double& b, double& theta) const {
93 }
94 
95 BaseCore::Jacobian Quadrupole::_dAssignToAxes(double& a, double& b, double& theta) const {
96  return BaseCore::_dAssignQuadrupoleToAxes(getIxx(), getIyy(), getIxy(), a, b, theta);
97 }
98 
99 void Quadrupole::_assignFromQuadrupole(double ixx, double iyy, double ixy) {
100  setIxx(ixx);
101  setIyy(iyy);
102  setIxy(ixy);
103 }
104 
105 BaseCore::Jacobian Quadrupole::_dAssignFromQuadrupole(double ixx, double iyy, double ixy) {
106  setIxx(ixx);
107  setIyy(iyy);
108  setIxy(ixy);
109  return Jacobian::Identity();
110 }
111 
112 void Quadrupole::_assignFromAxes(double a, double b, double theta) {
113  BaseCore::_assignAxesToQuadrupole(a, b, theta, _matrix(0, 0), _matrix(1, 1), _matrix(0, 1));
114  _matrix(1, 0) = _matrix(0, 1);
115 }
116 
117 BaseCore::Jacobian Quadrupole::_dAssignFromAxes(double a, double b, double theta) {
118  Jacobian r = BaseCore::_dAssignAxesToQuadrupole(a, b, theta, _matrix(0, 0), _matrix(1, 1), _matrix(0, 1));
119  _matrix(1, 0) = _matrix(0, 1);
120  return r;
121 }
122 } // namespace ellipses
123 } // namespace geom
124 } // namespace afw
125 } // 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
Eigen::Vector3d ParameterVector
Parameter vector type.
Definition: BaseCore.h:61
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:232
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:219
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:194
Eigen::Matrix3d Jacobian
Parameter Jacobian matrix type.
Definition: BaseCore.h:64
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:184
void _assignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:91
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:36
Jacobian _dAssignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:84
Quadrupole(double ixx=1.0, double iyy=1.0, double ixy=0.0, bool normalize=false)
Construct from parameter values.
Definition: Quadrupole.cc:60
void _assignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:99
std::string getName() const override
Return a string that identifies this parametrization.
Definition: Quadrupole.cc:34
void writeParameters(double *iter) const override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:54
void _assignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:78
Eigen::Matrix< double, 2, 2, Eigen::DontAlign > Matrix
Matrix type for the matrix representation of Quadrupole parameters.
Definition: Quadrupole.h:52
Jacobian _dAssignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:95
Jacobian _dAssignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:117
void _assignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:112
void readParameters(double const *iter) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:48
Jacobian _dAssignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Quadrupole.cc:105
Reports invalid arguments.
Definition: Runtime.h:66
A base class for image defects.