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
Axes.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<Axes> Axes::registrar;
33 
34 std::string Axes::getName() const { return "Axes"; }
35 
37  if (_vector[A] < 0 || _vector[B] < 0)
39  "Major and minor axes cannot be negative.");
40  if (_vector[A] < _vector[B]) {
41  std::swap(_vector[A], _vector[B]);
42  _vector[THETA] += M_PI_2;
43  }
44  if (_vector[THETA] > M_PI_2 || _vector[THETA] <= -M_PI_2) {
45  _vector[THETA] -= M_PI * std::ceil(_vector[THETA] / M_PI - 0.5);
46  }
47 }
48 
49 void Axes::readParameters(double const* iter) {
50  setA(*iter++);
51  setB(*iter++);
52  setTheta(*iter++);
53 }
54 
55 void Axes::writeParameters(double* iter) const {
56  *iter++ = getA();
57  *iter++ = getB();
58  *iter++ = getTheta();
59 }
60 
61 void Axes::_assignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
62  BaseCore::_assignAxesToQuadrupole(_vector[A], _vector[B], _vector[THETA], ixx, iyy, ixy);
63 }
64 
65 BaseCore::Jacobian Axes::_dAssignToQuadrupole(double& ixx, double& iyy, double& ixy) const {
66  return BaseCore::_dAssignAxesToQuadrupole(_vector[A], _vector[B], _vector[THETA], ixx, iyy, ixy);
67 }
68 
69 void Axes::_assignToAxes(double& a, double& b, double& theta) const {
70  a = _vector[A];
71  b = _vector[B];
72  theta = _vector[THETA];
73 }
74 
75 BaseCore::Jacobian Axes::_dAssignToAxes(double& a, double& b, double& theta) const {
76  a = _vector[A];
77  b = _vector[B];
78  theta = _vector[THETA];
79  return Jacobian::Identity();
80 }
81 
82 void Axes::_assignFromQuadrupole(double ixx, double iyy, double ixy) {
83  BaseCore::_assignQuadrupoleToAxes(ixx, iyy, ixy, _vector[A], _vector[B], _vector[THETA]);
84  normalize();
85 }
86 
87 BaseCore::Jacobian Axes::_dAssignFromQuadrupole(double ixx, double iyy, double ixy) {
88  return BaseCore::_dAssignQuadrupoleToAxes(ixx, iyy, ixy, _vector[A], _vector[B], _vector[THETA]);
89  normalize();
90 }
91 
92 void Axes::_assignFromAxes(double a, double b, double theta) {
93  _vector[A] = a;
94  _vector[B] = b;
95  _vector[THETA] = theta;
96 }
97 
98 BaseCore::Jacobian Axes::_dAssignFromAxes(double a, double b, double theta) {
99  _vector[A] = a;
100  _vector[B] = b;
101  _vector[THETA] = theta;
102  return Jacobian::Identity();
103 }
104 } // namespace ellipses
105 } // namespace geom
106 } // namespace afw
107 } // namespace lsst
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
#define M_PI
Definition: ListMatch.cc:31
table::Key< int > b
table::Key< int > a
T ceil(T... args)
void _assignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:92
std::string getName() const override
Return a string that identifies this parametrization.
Definition: Axes.cc:34
Jacobian _dAssignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:75
void setA(double a)
Definition: Axes.h:52
double const getTheta() const
Definition: Axes.h:57
Jacobian _dAssignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:87
void setB(double b)
Definition: Axes.h:55
void writeParameters(double *iter) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:55
double const getA() const
Definition: Axes.h:51
void _assignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:82
void normalize() override
Put the parameters into a "standard form", if possible, and throw InvalidParameterError if they canno...
Definition: Axes.cc:36
void setTheta(double theta)
Definition: Axes.h:58
double const getB() const
Definition: Axes.h:54
void readParameters(double const *iter) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:49
void _assignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:61
Jacobian _dAssignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:98
Jacobian _dAssignToQuadrupole(double &ixx, double &iyy, double &ixy) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:65
void _assignToAxes(double &a, double &b, double &theta) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:69
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
Reports invalid arguments.
Definition: Runtime.h:66
A base class for image defects.
T swap(T... args)