LSST Applications g013ef56533+2edba5e46d,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+d594c53611,g210f2d0738+d003cec0be,g25ecf2a47a+204e5a9ad5,g262e1987ae+16552835ac,g29ae962dfc+867dbde878,g2cef7863aa+aef1011c0b,g30d7c61c20+990af31dd2,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+c8e8c50746,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5b326b94bb+01d4acd216,g64539dfbff+d003cec0be,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g74acd417e5+4b21cac47b,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+b19788b7f5,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+3bd86eedd0,g98a1a72a9c+67cf96cfb3,g98df359435+688a06938e,gbf99507273+8c5ae1fdc5,gc2a301910b+d003cec0be,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+4b21cac47b,ge410e46f29+f459a6810c,ge41e95a9f2+d003cec0be,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
transforms.cc
Go to the documentation of this file.
1#include <limits>
2
4
6
7namespace lsst::gauss2d::fit {
8static const double INF = std::numeric_limits<double>::infinity();
9
10double LogitLimitedTransform::derivative(double x) const {
11 double y = (x - _limits->get_min()) / _range;
12 if (y == 1)
13 return INF;
14 else if (y == 0)
15 return -INF;
16 return (1 / y + 1 / (1 - y)) * _factor / _range;
17}
18
19double LogitLimitedTransform::forward(double x) const {
20 double min = _limits->get_min();
21 double max = _limits->get_max();
22 if (x == min)
24 else if (x == max)
26 double y = (x - min) / _range;
27 if (!(y < 1) || !(y > 0)) return nan("");
28 return log(y / (1 - y)) * _factor;
29}
30
31double LogitLimitedTransform::reverse(double x) const {
32 double y = -x * _factor;
33 // math.log(np.finfo(np.float64) = 709.782712893384
34 // things will go badly well before then
35 if (y > 709.7827) return _limits->get_min();
36 y = 1 + exp(y);
37 return _range / y + _limits->get_min();
38}
39} // namespace lsst::gauss2d::fit
double forward(double x) const override
Return the transformed value of x.
Definition transforms.cc:19
double derivative(double x) const override
Return the derivative of this tranform at the value x.
Definition transforms.cc:10
double reverse(double x) const override
Return the original value of x given a transformed value.
Definition transforms.cc:31
T infinity(T... args)