LSSTApplications  19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
LSSTDataManagementBasePackage
Distortion.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  */
28 
29 namespace lsst {
30 namespace afw {
31 namespace geom {
32 namespace ellipses {
33 
34 double Distortion::getAxisRatio() const {
35  double e = getE();
36  return std::sqrt((1.0 - e) / (1.0 + e));
37 }
38 
40  double eta = other.getE();
41  if (eta < 1E-8) {
42  _complex = other.getComplex() * (1.0 - eta * eta / 3.0);
43  } else {
44  double delta = std::tanh(eta);
45  _complex = other.getComplex() * delta / eta;
46  }
47  return *this;
48 }
49 
51  double g = other.getE();
52  _complex = other.getComplex() * 2.0 / (1 + g * g);
53  return *this;
54 }
55 
57  Jacobian result = Jacobian::Zero();
58  double eta = other.getE();
59  double alpha, beta;
60  if (eta < 1E-8) {
61  alpha = (1.0 - eta * eta / 3.0);
62  beta = -2.0 / 3.0;
63  } else {
64  double delta = std::tanh(eta);
65  alpha = delta / eta;
66  beta = (1.0 - delta * delta - alpha) / (eta * eta);
67  }
68  _complex = other.getComplex() * alpha;
69  result(0, 0) = alpha + other.getE1() * other.getE1() * beta;
70  result(1, 1) = alpha + other.getE2() * other.getE2() * beta;
71  result(1, 0) = result(0, 1) = other.getE1() * other.getE2() * beta;
72  return result;
73 }
74 
76  Jacobian result = Jacobian::Zero();
77  double g = other.getE();
78  double alpha = 2.0 / (1 + g * g);
79  double beta = -alpha * alpha;
80  _complex = other.getComplex() * alpha;
81  result(0, 0) = alpha + other.getE1() * other.getE1() * beta;
82  result(1, 1) = alpha + other.getE2() * other.getE2() * beta;
83  result(1, 0) = result(0, 1) = other.getE1() * other.getE2() * beta;
84  return result;
85 }
86 
88  if (getE() > 1.0) {
90  "Distortion magnitude cannot be greater than one.");
91  }
92 }
93 } // namespace ellipses
94 } // namespace geom
95 } // namespace afw
96 } // namespace lsst
lsst::afw::geom::ellipses::Distortion::getAxisRatio
double getAxisRatio() const
Definition: Distortion.cc:34
ellipses
Distortion.h
lsst::afw::geom::ellipses::Distortion::normalize
void normalize()
Definition: Distortion.cc:87
lsst::afw
Definition: imageAlgorithm.dox:1
lsst::afw::geom::ellipses::Distortion::dAssign
Jacobian dAssign(Distortion const &other)
Definition: Distortion.h:70
std::sqrt
T sqrt(T... args)
lsst::afw::geom::ellipses::Distortion::operator=
Distortion & operator=(Distortion const &other)
Definition: Distortion.h:59
lsst::afw::geom::ellipses::ConformalShear
A logarithmic complex ellipticity with magnitude .
Definition: ConformalShear.h:45
ConformalShear.h
lsst::afw::geom::ellipses::detail::EllipticityBase::Jacobian
Eigen::Matrix2d Jacobian
Definition: EllipticityBase.h:47
other
ItemVariant const * other
Definition: Schema.cc:56
lsst::afw::geom::ellipses::detail::EllipticityBase::_complex
std::complex< double > _complex
Definition: EllipticityBase.h:91
lsst::afw::geom::ellipses::ReducedShear
A complex ellipticity with magnitude .
Definition: ReducedShear.h:45
result
py::object result
Definition: _schema.cc:429
lsst::afw::geom::ellipses::detail::EllipticityBase::getE
double getE() const
Definition: EllipticityBase.h:75
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
lsst::geom
Definition: AffineTransform.h:36
lsst.pex::exceptions::InvalidParameterError
Reports invalid arguments.
Definition: Runtime.h:66
std::tanh
T tanh(T... args)
ReducedShear.h
lsst::afw::geom::ellipses::Distortion
A complex ellipticity with magnitude .
Definition: Distortion.h:44
BaseCore.h