LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
ConformalShear.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 */
24#include <cmath>
25
29
30namespace lsst {
31namespace afw {
32namespace geom {
33namespace ellipses {
34
36 double e = getE();
37 return std::exp(-e);
38}
39
41 double delta = other.getE();
42 if (delta < 1E-8) {
43 _complex = other.getComplex() * (1.0 + delta * delta / 3.0);
44 } else {
45 double eta = std::atanh(delta);
46 _complex = other.getComplex() * eta / delta;
47 }
48 return *this;
49}
50
52 double g = other.getE();
53 if (g < 1E-8) {
54 _complex = other.getComplex() * 2.0 * (1.0 + g * g / 3.0);
55 } else {
56 double eta = 2.0 * std::atanh(g);
57 _complex = other.getComplex() * eta / g;
58 }
59 return *this;
60}
61
63 Jacobian result = Jacobian::Zero();
64 double delta = other.getE();
65 double alpha, beta;
66 if (delta < 1E-8) {
67 alpha = 1.0 + delta * delta / 3.0;
68 beta = 2.0 / 3.0;
69 } else {
70 double eta = std::atanh(delta);
71 alpha = eta / delta;
72 beta = (1.0 / (1.0 - delta * delta) - alpha) / (delta * delta);
73 }
74 _complex = other.getComplex() * alpha;
75 result(0, 0) = alpha + other.getE1() * other.getE1() * beta;
76 result(1, 1) = alpha + other.getE2() * other.getE2() * beta;
77 result(1, 0) = result(0, 1) = other.getE1() * other.getE2() * beta;
78 return result;
79}
80
82 Jacobian result = Jacobian::Zero();
83 double g = other.getE();
84 double alpha, beta;
85 if (g < 1E-8) {
86 alpha = 2.0 * (1.0 + g * g / 3.0);
87 beta = 4.0 / 3.0;
88 } else {
89 double eta = 2.0 * std::atanh(g);
90 alpha = eta / g;
91 beta = 1.0 * (2.0 / (1.0 - g * g) - alpha) / (g * g);
92 }
93 _complex = other.getComplex() * alpha;
94 result(0, 0) = alpha + other.getE1() * other.getE1() * beta;
95 result(1, 1) = alpha + other.getE2() * other.getE2() * beta;
96 result(1, 0) = result(0, 1) = other.getE1() * other.getE2() * beta;
97 return result;
98}
99} // namespace ellipses
100} // namespace geom
101} // namespace afw
102} // namespace lsst
py::object result
Definition _schema.cc:429
T atanh(T... args)
A logarithmic complex ellipticity with magnitude .
Jacobian dAssign(ConformalShear const &other)
ConformalShear & operator=(ConformalShear const &other)
A complex ellipticity with magnitude .
Definition Distortion.h:44
A complex ellipticity with magnitude .
T exp(T... args)