| Alternative 1 | |
|---|---|
| Error | 99.4% |
| Cost | 13248.00 |
\[\frac{\frac{1}{s}}{\left(2 + e^{\frac{\left|x\right|}{s}}\right) + e^{\frac{\left|x\right|}{-s}}}
\]
(FPCore (x s) :precision binary32 (/ (exp (/ (- (fabs x)) s)) (* (* s (+ 1.0 (exp (/ (- (fabs x)) s)))) (+ 1.0 (exp (/ (- (fabs x)) s))))))
(FPCore (x s) :precision binary32 (let* ((t_0 (/ (fabs x) s))) (/ 1.0 (* s (+ (exp (- t_0)) (+ 2.0 (exp t_0)))))))
float code(float x, float s) {
return expf((-fabsf(x) / s)) / ((s * (1.0f + expf((-fabsf(x) / s)))) * (1.0f + expf((-fabsf(x) / s))));
}
float code(float x, float s) {
float t_0 = fabsf(x) / s;
return 1.0f / (s * (expf(-t_0) + (2.0f + expf(t_0))));
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = exp((-abs(x) / s)) / ((s * (1.0e0 + exp((-abs(x) / s)))) * (1.0e0 + exp((-abs(x) / s))))
end function
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: t_0
t_0 = abs(x) / s
code = 1.0e0 / (s * (exp(-t_0) + (2.0e0 + exp(t_0))))
end function
function code(x, s) return Float32(exp(Float32(Float32(-abs(x)) / s)) / Float32(Float32(s * Float32(Float32(1.0) + exp(Float32(Float32(-abs(x)) / s)))) * Float32(Float32(1.0) + exp(Float32(Float32(-abs(x)) / s))))) end
function code(x, s) t_0 = Float32(abs(x) / s) return Float32(Float32(1.0) / Float32(s * Float32(exp(Float32(-t_0)) + Float32(Float32(2.0) + exp(t_0))))) end
function tmp = code(x, s) tmp = exp((-abs(x) / s)) / ((s * (single(1.0) + exp((-abs(x) / s)))) * (single(1.0) + exp((-abs(x) / s)))); end
function tmp = code(x, s) t_0 = abs(x) / s; tmp = single(1.0) / (s * (exp(-t_0) + (single(2.0) + exp(t_0)))); end
\frac{e^{\frac{-\left|x\right|}{s}}}{\left(s \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\begin{array}{l}
t_0 := \frac{\left|x\right|}{s}\\
\frac{1}{s \cdot \left(e^{-t_0} + \left(2 + e^{t_0}\right)\right)}
\end{array}
Results
Initial program 99.5
Simplified99.4
[Start]99.5 | \[ \frac{e^{\frac{-\left|x\right|}{s}}}{\left(s \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\] |
|---|---|
associate-*l* [=>]99.5 | \[ \frac{e^{\frac{-\left|x\right|}{s}}}{\color{blue}{s \cdot \left(\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)\right)}}
\] |
associate-/r* [=>]99.5 | \[ \color{blue}{\frac{\frac{e^{\frac{-\left|x\right|}{s}}}{s}}{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}}
\] |
distribute-frac-neg [=>]99.5 | \[ \frac{\frac{e^{\color{blue}{-\frac{\left|x\right|}{s}}}}{s}}{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\] |
exp-neg [=>]99.5 | \[ \frac{\frac{\color{blue}{\frac{1}{e^{\frac{\left|x\right|}{s}}}}}{s}}{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\] |
associate-/l/ [=>]99.5 | \[ \frac{\color{blue}{\frac{1}{s \cdot e^{\frac{\left|x\right|}{s}}}}}{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\] |
associate-/r* [=>]99.4 | \[ \frac{\color{blue}{\frac{\frac{1}{s}}{e^{\frac{\left|x\right|}{s}}}}}{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}
\] |
associate-/r* [<=]99.4 | \[ \color{blue}{\frac{\frac{1}{s}}{e^{\frac{\left|x\right|}{s}} \cdot \left(\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)\right)}}
\] |
Taylor expanded in s around 0 99.6
Final simplification99.6
| Alternative 1 | |
|---|---|
| Error | 99.4% |
| Cost | 13248.00 |
| Alternative 2 | |
|---|---|
| Error | 99.5% |
| Cost | 6944.00 |
| Alternative 3 | |
|---|---|
| Error | 99.5% |
| Cost | 6912.00 |
| Alternative 4 | |
|---|---|
| Error | 99.5% |
| Cost | 6912.00 |
| Alternative 5 | |
|---|---|
| Error | 99.5% |
| Cost | 6880.00 |
| Alternative 6 | |
|---|---|
| Error | 96.2% |
| Cost | 6688.00 |
| Alternative 7 | |
|---|---|
| Error | 94.7% |
| Cost | 6656.00 |
| Alternative 8 | |
|---|---|
| Error | 95.1% |
| Cost | 3620.00 |
| Alternative 9 | |
|---|---|
| Error | 80.3% |
| Cost | 3556.00 |
| Alternative 10 | |
|---|---|
| Error | 95.0% |
| Cost | 3556.00 |
| Alternative 11 | |
|---|---|
| Error | 95.0% |
| Cost | 3556.00 |
| Alternative 12 | |
|---|---|
| Error | 66.1% |
| Cost | 352.00 |
| Alternative 13 | |
|---|---|
| Error | 63.2% |
| Cost | 297.00 |
| Alternative 14 | |
|---|---|
| Error | 27.4% |
| Cost | 96.00 |
herbie shell --seed 2023093
(FPCore (x s)
:name "Logistic distribution"
:precision binary32
:pre (and (<= 0.0 s) (<= s 1.0651631))
(/ (exp (/ (- (fabs x)) s)) (* (* s (+ 1.0 (exp (/ (- (fabs x)) s)))) (+ 1.0 (exp (/ (- (fabs x)) s))))))