35 return (1.0 - e) / (1.0 + e);
39 double delta = other.
getE();
43 double g = (1.0 -
std::sqrt(1.0 - delta * delta)) / delta;
50 double eta = other.
getE();
62 double delta = other.
getE();
63 double s =
std::sqrt(1.0 - delta * delta);
66 alpha = 0.5 + 0.125 * delta * delta;
69 alpha = (1.0 - s) / (delta * delta);
70 beta = (2.0 * alpha - 1.0) / (delta * delta * s);
73 result(0, 0) = alpha + other.
getE1() * other.
getE1() * beta;
74 result(1, 1) = alpha + other.
getE2() * other.
getE2() * beta;
75 result(1, 0) = result(0, 1) = other.
getE1() * other.
getE2() * beta;
81 double eta = other.
getE();
84 alpha = 0.5 - eta * eta / 24.0;
89 beta = (0.5 * (1.0 - g * g) - alpha) / (eta * eta);
92 result(0, 0) = alpha + other.
getE1() * other.
getE1() * beta;
93 result(1, 1) = alpha + other.
getE2() * other.
getE2() * beta;
94 result(1, 0) = result(0, 1) = other.
getE1() * other.
getE2() * beta;
A complex ellipticity with magnitude .
A complex ellipticity with magnitude .
ReducedShear & operator=(ReducedShear const &other)
Jacobian dAssign(ReducedShear const &other)
double getAxisRatio() const
std::complex< double > _complex
std::complex< double > & getComplex()
A base class for image defects.