LSSTApplications  20.0.0
LSSTDataManagementBasePackage
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
lsst::afw::geom::ellipses::Axes::getName
std::string getName() const override
Return a string that identifies this parametrization.
Definition: Axes.cc:34
lsst::afw::geom::ellipses::Axes::readParameters
void readParameters(double const *iter) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:49
std::string
STL class.
ellipses
lsst::afw::geom::ellipses::Axes::getB
double const getB() const
Definition: Axes.h:54
lsst::afw::geom::ellipses::Axes::normalize
void normalize() override
Put the parameters into a "standard form", if possible, and throw InvalidParameterError if they canno...
Definition: Axes.cc:36
lsst::afw::geom::ellipses::Axes::_assignFromQuadrupole
void _assignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:82
Axes.h
lsst::afw::geom::ellipses::Axes::_dAssignToAxes
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
lsst::afw::geom::ellipses::Axes::writeParameters
void writeParameters(double *iter) const override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:55
lsst::afw::geom::ellipses::BaseCore::_dAssignAxesToQuadrupole
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
lsst::afw::geom::ellipses::Axes::_assignToQuadrupole
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
lsst::afw
Definition: imageAlgorithm.dox:1
Quadrupole.h
lsst::afw::geom::ellipses::Axes::getTheta
double const getTheta() const
Definition: Axes.h:57
lsst::afw::geom::ellipses::Axes::THETA
@ THETA
Definition: Axes.h:49
lsst::afw::geom::ellipses::Axes::A
@ A
Definition: Axes.h:49
lsst::afw::geom::ellipses::Axes::getA
double const getA() const
Definition: Axes.h:51
lsst::afw::geom::ellipses::Axes::_dAssignToQuadrupole
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
lsst::afw::geom::ellipses::BaseCore::_dAssignQuadrupoleToAxes
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
lsst::afw::geom::ellipses::Axes::_dAssignFromAxes
Jacobian _dAssignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:98
std::ceil
T ceil(T... args)
b
table::Key< int > b
Definition: TransmissionCurve.cc:467
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
lsst::afw::geom::ellipses::Axes::setA
void setA(double a)
Definition: Axes.h:52
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
std::swap
T swap(T... args)
lsst::afw::geom::ellipses::Axes::setTheta
void setTheta(double theta)
Definition: Axes.h:58
lsst::geom
Definition: geomOperators.dox:4
lsst::afw::geom::ellipses::BaseCore::Jacobian
Eigen::Matrix3d Jacobian
Parameter Jacobian matrix type.
Definition: BaseCore.h:64
lsst::pex::exceptions::InvalidParameterError
Reports invalid arguments.
Definition: Runtime.h:66
lsst::afw::geom::ellipses::BaseCore::_assignAxesToQuadrupole
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
a
table::Key< int > a
Definition: TransmissionCurve.cc:466
lsst::afw::geom::ellipses::Axes::setB
void setB(double b)
Definition: Axes.h:55
M_PI
#define M_PI
Definition: ListMatch.cc:31
lsst::afw::geom::ellipses::Axes::_dAssignFromQuadrupole
Jacobian _dAssignFromQuadrupole(double ixx, double iyy, double ixy) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:87
lsst::afw::geom::ellipses::Axes::_assignFromAxes
void _assignFromAxes(double a, double b, double theta) override
Return the size of the bounding box for the ellipse core.
Definition: Axes.cc:92
astshim.fitsChanContinued.iter
def iter(self)
Definition: fitsChanContinued.py:88
lsst::afw::geom::ellipses::Axes::_assignToAxes
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
lsst::afw::geom::ellipses::BaseCore::_assignQuadrupoleToAxes
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
lsst::afw::geom::ellipses::Axes::B
@ B
Definition: Axes.h:49