LSSTApplications  10.0-2-g4f67435,11.0.rc2+1,11.0.rc2+12,11.0.rc2+3,11.0.rc2+4,11.0.rc2+5,11.0.rc2+6,11.0.rc2+7,11.0.rc2+8
LSSTDataManagementBasePackage
EllipticityBase.h
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 
25 #ifndef LSST_AFW_GEOM_ELLIPSES_EllipticityBase_h_INCLUDED
26 #define LSST_AFW_GEOM_ELLIPSES_EllipticityBase_h_INCLUDED
27 
28 #include "Eigen/Core"
29 #include <complex>
30 
31 namespace lsst { namespace afw { namespace geom { namespace ellipses {
32 
33 namespace detail {
34 
43 public:
44 
45  typedef Eigen::Matrix2d Jacobian;
46 
47  enum ParameterEnum { E1=0, E2=1 };
48 
49  std::complex<double> & getComplex() { return _complex; }
50 
51  std::complex<double> const & getComplex() const { return _complex; }
52 
53  void setComplex(std::complex<double> const & v) { _complex = v; }
54 
55  double getE1() const { return _complex.real(); }
56  void setE1(double e1) {
57 #if __cplusplus < 201103L
58  _complex = std::complex<double>(e1, _complex.imag());
59 #else
60  _complex.real(e1);
61 #endif
62  }
63 
64  double getE2() const { return _complex.imag(); }
65  void setE2(double e2) {
66 #if __cplusplus < 201103L
67  _complex = std::complex<double>(_complex.real(), e2);
68 #else
69  _complex.imag(e2);
70 #endif
71  }
72 
73  double getE() const { return std::sqrt(std::norm(_complex)); }
74  void setE(double e) { _complex *= e / getE(); }
75 
76  double getTheta() const { return 0.5 * std::arg(_complex); }
77 
78 protected:
79 
80  explicit EllipticityBase(std::complex<double> const & complex) : _complex(complex) {}
81 
82  explicit EllipticityBase(double e1=0.0, double e2=0.0) : _complex(e1, e2) {}
83 
84  std::complex<double> _complex;
85 };
86 
87 } // namespace detail
88 
89 }}}} // namespace lsst::afw::geom::ellipses
90 
91 #endif // !LSST_AFW_GEOM_ELLIPSES_EllipticityBase_h_INCLUDED
EllipticityBase(std::complex< double > const &complex)
T norm(const T &x)
Definition: Integrate.h:191
std::complex< double > const & getComplex() const
EllipticityBase(double e1=0.0, double e2=0.0)
EllipticityBase is a base class for complex ellipticity types.
void setComplex(std::complex< double > const &v)