| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 6976 |
\[\frac{1}{\left(1 + \left(1 + \mathsf{expm1}\left(\frac{x}{s}\right)\right)\right) \cdot \left(s \cdot \left(1 + e^{-\frac{x}{s}}\right)\right)}
\]
(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 (/ (/ 1.0 s) (+ (exp (/ (fabs x) s)) (+ (exp (/ (fabs x) (- s))) 2.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) {
return (1.0f / s) / (expf((fabsf(x) / s)) + (expf((fabsf(x) / -s)) + 2.0f));
}
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
code = (1.0e0 / s) / (exp((abs(x) / s)) + (exp((abs(x) / -s)) + 2.0e0))
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) return Float32(Float32(Float32(1.0) / s) / Float32(exp(Float32(abs(x) / s)) + Float32(exp(Float32(abs(x) / Float32(-s))) + Float32(2.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) tmp = (single(1.0) / s) / (exp((abs(x) / s)) + (exp((abs(x) / -s)) + single(2.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)}
\frac{\frac{1}{s}}{e^{\frac{\left|x\right|}{s}} + \left(e^{\frac{\left|x\right|}{-s}} + 2\right)}
Results
Initial program 0.2
Simplified0.2
[Start]0.2 | \[ \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)}
\] |
|---|---|
*-lft-identity [<=]0.2 | \[ \color{blue}{1 \cdot \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-*r/ [=>]0.2 | \[ \color{blue}{\frac{1 \cdot 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* [=>]0.2 | \[ \frac{1 \cdot 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)}}
\] |
times-frac [=>]0.2 | \[ \color{blue}{\frac{1}{s} \cdot \frac{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/ [=>]0.2 | \[ \color{blue}{\frac{\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-/l* [=>]0.2 | \[ \color{blue}{\frac{\frac{1}{s}}{\frac{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}{e^{\frac{-\left|x\right|}{s}}}}}
\] |
distribute-frac-neg [=>]0.2 | \[ \frac{\frac{1}{s}}{\frac{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}{e^{\color{blue}{-\frac{\left|x\right|}{s}}}}}
\] |
exp-neg [=>]0.2 | \[ \frac{\frac{1}{s}}{\frac{\left(1 + e^{\frac{-\left|x\right|}{s}}\right) \cdot \left(1 + e^{\frac{-\left|x\right|}{s}}\right)}{\color{blue}{\frac{1}{e^{\frac{\left|x\right|}{s}}}}}}
\] |
Final simplification0.2
| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 6976 |
| Alternative 2 | |
|---|---|
| Error | 0.2 |
| Cost | 6912 |
| Alternative 3 | |
|---|---|
| Error | 0.2 |
| Cost | 6912 |
| Alternative 4 | |
|---|---|
| Error | 0.2 |
| Cost | 6880 |
| Alternative 5 | |
|---|---|
| Error | 1.2 |
| Cost | 6688 |
| Alternative 6 | |
|---|---|
| Error | 1.7 |
| Cost | 6656 |
| Alternative 7 | |
|---|---|
| Error | 1.5 |
| Cost | 3588 |
| Alternative 8 | |
|---|---|
| Error | 3.1 |
| Cost | 3556 |
| Alternative 9 | |
|---|---|
| Error | 1.6 |
| Cost | 3556 |
| Alternative 10 | |
|---|---|
| Error | 3.5 |
| Cost | 3492 |
| Alternative 11 | |
|---|---|
| Error | 5.4 |
| Cost | 553 |
| Alternative 12 | |
|---|---|
| Error | 5.5 |
| Cost | 480 |
| Alternative 13 | |
|---|---|
| Error | 6.1 |
| Cost | 425 |
| Alternative 14 | |
|---|---|
| Error | 7.3 |
| Cost | 416 |
| Alternative 15 | |
|---|---|
| Error | 11.7 |
| Cost | 360 |
| Alternative 16 | |
|---|---|
| Error | 11.7 |
| Cost | 297 |
| Alternative 17 | |
|---|---|
| Error | 23.3 |
| Cost | 96 |
herbie shell --seed 2023066
(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))))))