
(FPCore (x s) :precision binary32 (let* ((t_0 (exp (/ (- (fabs x)) s))) (t_1 (+ 1.0 t_0))) (/ t_0 (* (* s t_1) t_1))))
float code(float x, float s) {
float t_0 = expf((-fabsf(x) / s));
float t_1 = 1.0f + t_0;
return t_0 / ((s * t_1) * t_1);
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: t_0
real(4) :: t_1
t_0 = exp((-abs(x) / s))
t_1 = 1.0e0 + t_0
code = t_0 / ((s * t_1) * t_1)
end function
function code(x, s) t_0 = exp(Float32(Float32(-abs(x)) / s)) t_1 = Float32(Float32(1.0) + t_0) return Float32(t_0 / Float32(Float32(s * t_1) * t_1)) end
function tmp = code(x, s) t_0 = exp((-abs(x) / s)); t_1 = single(1.0) + t_0; tmp = t_0 / ((s * t_1) * t_1); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\frac{-\left|x\right|}{s}}\\
t_1 := 1 + t\_0\\
\frac{t\_0}{\left(s \cdot t\_1\right) \cdot t\_1}
\end{array}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x s) :precision binary32 (let* ((t_0 (exp (/ (- (fabs x)) s))) (t_1 (+ 1.0 t_0))) (/ t_0 (* (* s t_1) t_1))))
float code(float x, float s) {
float t_0 = expf((-fabsf(x) / s));
float t_1 = 1.0f + t_0;
return t_0 / ((s * t_1) * t_1);
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: t_0
real(4) :: t_1
t_0 = exp((-abs(x) / s))
t_1 = 1.0e0 + t_0
code = t_0 / ((s * t_1) * t_1)
end function
function code(x, s) t_0 = exp(Float32(Float32(-abs(x)) / s)) t_1 = Float32(Float32(1.0) + t_0) return Float32(t_0 / Float32(Float32(s * t_1) * t_1)) end
function tmp = code(x, s) t_0 = exp((-abs(x) / s)); t_1 = single(1.0) + t_0; tmp = t_0 / ((s * t_1) * t_1); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\frac{-\left|x\right|}{s}}\\
t_1 := 1 + t\_0\\
\frac{t\_0}{\left(s \cdot t\_1\right) \cdot t\_1}
\end{array}
\end{array}
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (let* ((t_0 (/ x_m (- s)))) (/ (pow E t_0) (* (+ 1.0 (exp t_0)) (+ s (/ s (exp (/ x_m s))))))))
x_m = fabs(x);
float code(float x_m, float s) {
float t_0 = x_m / -s;
return powf(((float) M_E), t_0) / ((1.0f + expf(t_0)) * (s + (s / expf((x_m / s)))));
}
x_m = abs(x) function code(x_m, s) t_0 = Float32(x_m / Float32(-s)) return Float32((Float32(exp(1)) ^ t_0) / Float32(Float32(Float32(1.0) + exp(t_0)) * Float32(s + Float32(s / exp(Float32(x_m / s)))))) end
x_m = abs(x); function tmp = code(x_m, s) t_0 = x_m / -s; tmp = (single(2.71828182845904523536) ^ t_0) / ((single(1.0) + exp(t_0)) * (s + (s / exp((x_m / s))))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := \frac{x\_m}{-s}\\
\frac{{e}^{t\_0}}{\left(1 + e^{t\_0}\right) \cdot \left(s + \frac{s}{e^{\frac{x\_m}{s}}}\right)}
\end{array}
\end{array}
Initial program 99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
distribute-lft-in99.3%
*-rgt-identity99.3%
+-commutative99.3%
Simplified99.3%
add-sqr-sqrt99.1%
associate-/l*99.1%
remove-double-neg99.1%
add-sqr-sqrt99.1%
sqrt-unprod95.6%
sqr-neg95.6%
sqrt-unprod-0.0%
add-sqr-sqrt94.3%
frac-2neg94.3%
add-sqr-sqrt-0.0%
sqrt-unprod97.4%
sqr-neg97.4%
sqrt-unprod99.1%
add-sqr-sqrt99.1%
add-sqr-sqrt49.0%
fabs-sqr49.0%
add-sqr-sqrt97.6%
Applied egg-rr97.6%
associate-*r/97.6%
rem-square-sqrt97.8%
Simplified97.8%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr97.7%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified97.7%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr60.4%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified60.7%
*-un-lft-identity60.7%
exp-prod60.7%
Applied egg-rr60.7%
exp-1-e60.7%
Simplified60.7%
Final simplification60.7%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (let* ((t_0 (exp (/ x_m (- s))))) (/ t_0 (* (+ 1.0 t_0) (+ s (/ s (exp (/ x_m s))))))))
x_m = fabs(x);
float code(float x_m, float s) {
float t_0 = expf((x_m / -s));
return t_0 / ((1.0f + t_0) * (s + (s / expf((x_m / s)))));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
real(4) :: t_0
t_0 = exp((x_m / -s))
code = t_0 / ((1.0e0 + t_0) * (s + (s / exp((x_m / s)))))
end function
x_m = abs(x) function code(x_m, s) t_0 = exp(Float32(x_m / Float32(-s))) return Float32(t_0 / Float32(Float32(Float32(1.0) + t_0) * Float32(s + Float32(s / exp(Float32(x_m / s)))))) end
x_m = abs(x); function tmp = code(x_m, s) t_0 = exp((x_m / -s)); tmp = t_0 / ((single(1.0) + t_0) * (s + (s / exp((x_m / s))))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := e^{\frac{x\_m}{-s}}\\
\frac{t\_0}{\left(1 + t\_0\right) \cdot \left(s + \frac{s}{e^{\frac{x\_m}{s}}}\right)}
\end{array}
\end{array}
Initial program 99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
distribute-lft-in99.3%
*-rgt-identity99.3%
+-commutative99.3%
Simplified99.3%
add-sqr-sqrt99.1%
associate-/l*99.1%
remove-double-neg99.1%
add-sqr-sqrt99.1%
sqrt-unprod95.6%
sqr-neg95.6%
sqrt-unprod-0.0%
add-sqr-sqrt94.3%
frac-2neg94.3%
add-sqr-sqrt-0.0%
sqrt-unprod97.4%
sqr-neg97.4%
sqrt-unprod99.1%
add-sqr-sqrt99.1%
add-sqr-sqrt49.0%
fabs-sqr49.0%
add-sqr-sqrt97.6%
Applied egg-rr97.6%
associate-*r/97.6%
rem-square-sqrt97.8%
Simplified97.8%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr97.7%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified97.7%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr60.4%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified60.7%
Final simplification60.7%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (let* ((t_0 (exp (/ x_m (- s))))) (/ t_0 (* (+ 1.0 t_0) (+ s (/ s (+ 1.0 (/ x_m s))))))))
x_m = fabs(x);
float code(float x_m, float s) {
float t_0 = expf((x_m / -s));
return t_0 / ((1.0f + t_0) * (s + (s / (1.0f + (x_m / s)))));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
real(4) :: t_0
t_0 = exp((x_m / -s))
code = t_0 / ((1.0e0 + t_0) * (s + (s / (1.0e0 + (x_m / s)))))
end function
x_m = abs(x) function code(x_m, s) t_0 = exp(Float32(x_m / Float32(-s))) return Float32(t_0 / Float32(Float32(Float32(1.0) + t_0) * Float32(s + Float32(s / Float32(Float32(1.0) + Float32(x_m / s)))))) end
x_m = abs(x); function tmp = code(x_m, s) t_0 = exp((x_m / -s)); tmp = t_0 / ((single(1.0) + t_0) * (s + (s / (single(1.0) + (x_m / s))))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
t_0 := e^{\frac{x\_m}{-s}}\\
\frac{t\_0}{\left(1 + t\_0\right) \cdot \left(s + \frac{s}{1 + \frac{x\_m}{s}}\right)}
\end{array}
\end{array}
Initial program 99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
distribute-lft-in99.3%
*-rgt-identity99.3%
+-commutative99.3%
Simplified99.3%
add-sqr-sqrt99.1%
associate-/l*99.1%
remove-double-neg99.1%
add-sqr-sqrt99.1%
sqrt-unprod95.6%
sqr-neg95.6%
sqrt-unprod-0.0%
add-sqr-sqrt94.3%
frac-2neg94.3%
add-sqr-sqrt-0.0%
sqrt-unprod97.4%
sqr-neg97.4%
sqrt-unprod99.1%
add-sqr-sqrt99.1%
add-sqr-sqrt49.0%
fabs-sqr49.0%
add-sqr-sqrt97.6%
Applied egg-rr97.6%
associate-*r/97.6%
rem-square-sqrt97.8%
Simplified97.8%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr97.7%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified97.7%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr60.4%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified60.7%
Taylor expanded in x around 0 58.1%
Final simplification58.1%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ (pow E (/ x_m (- s))) (* (+ s (/ s (exp (/ x_m s)))) 2.0)))
x_m = fabs(x);
float code(float x_m, float s) {
return powf(((float) M_E), (x_m / -s)) / ((s + (s / expf((x_m / s)))) * 2.0f);
}
x_m = abs(x) function code(x_m, s) return Float32((Float32(exp(1)) ^ Float32(x_m / Float32(-s))) / Float32(Float32(s + Float32(s / exp(Float32(x_m / s)))) * Float32(2.0))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = (single(2.71828182845904523536) ^ (x_m / -s)) / ((s + (s / exp((x_m / s)))) * single(2.0)); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{{e}^{\left(\frac{x\_m}{-s}\right)}}{\left(s + \frac{s}{e^{\frac{x\_m}{s}}}\right) \cdot 2}
\end{array}
Initial program 99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
distribute-lft-in99.3%
*-rgt-identity99.3%
+-commutative99.3%
Simplified99.3%
add-sqr-sqrt99.1%
associate-/l*99.1%
remove-double-neg99.1%
add-sqr-sqrt99.1%
sqrt-unprod95.6%
sqr-neg95.6%
sqrt-unprod-0.0%
add-sqr-sqrt94.3%
frac-2neg94.3%
add-sqr-sqrt-0.0%
sqrt-unprod97.4%
sqr-neg97.4%
sqrt-unprod99.1%
add-sqr-sqrt99.1%
add-sqr-sqrt49.0%
fabs-sqr49.0%
add-sqr-sqrt97.6%
Applied egg-rr97.6%
associate-*r/97.6%
rem-square-sqrt97.8%
Simplified97.8%
Taylor expanded in s around inf 94.8%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr56.2%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified56.3%
*-un-lft-identity60.7%
exp-prod60.7%
Applied egg-rr56.3%
exp-1-e60.7%
Simplified56.3%
Final simplification56.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (* 0.5 (/ (exp (/ x_m (- s))) (+ s (/ s (exp (/ x_m s)))))))
x_m = fabs(x);
float code(float x_m, float s) {
return 0.5f * (expf((x_m / -s)) / (s + (s / expf((x_m / s)))));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = 0.5e0 * (exp((x_m / -s)) / (s + (s / exp((x_m / s)))))
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(0.5) * Float32(exp(Float32(x_m / Float32(-s))) / Float32(s + Float32(s / exp(Float32(x_m / s)))))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = single(0.5) * (exp((x_m / -s)) / (s + (s / exp((x_m / s))))); end
\begin{array}{l}
x_m = \left|x\right|
\\
0.5 \cdot \frac{e^{\frac{x\_m}{-s}}}{s + \frac{s}{e^{\frac{x\_m}{s}}}}
\end{array}
Initial program 99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
distribute-lft-in99.3%
*-rgt-identity99.3%
+-commutative99.3%
Simplified99.3%
add-sqr-sqrt99.1%
associate-/l*99.1%
remove-double-neg99.1%
add-sqr-sqrt99.1%
sqrt-unprod95.6%
sqr-neg95.6%
sqrt-unprod-0.0%
add-sqr-sqrt94.3%
frac-2neg94.3%
add-sqr-sqrt-0.0%
sqrt-unprod97.4%
sqr-neg97.4%
sqrt-unprod99.1%
add-sqr-sqrt99.1%
add-sqr-sqrt49.0%
fabs-sqr49.0%
add-sqr-sqrt97.6%
Applied egg-rr97.6%
associate-*r/97.6%
rem-square-sqrt97.8%
Simplified97.8%
Taylor expanded in s around inf 94.8%
distribute-frac-neg97.8%
add-sqr-sqrt49.1%
fabs-sqr49.1%
add-sqr-sqrt97.7%
rec-exp97.7%
Applied egg-rr56.2%
rec-exp97.7%
distribute-frac-neg97.7%
Simplified56.3%
Taylor expanded in x around inf 56.3%
associate-*r/56.3%
mul-1-neg56.3%
Simplified56.3%
Final simplification56.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ (* 0.5 (/ 1.0 (+ 1.0 (exp (/ x_m s))))) s))
x_m = fabs(x);
float code(float x_m, float s) {
return (0.5f * (1.0f / (1.0f + expf((x_m / s))))) / s;
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = (0.5e0 * (1.0e0 / (1.0e0 + exp((x_m / s))))) / s
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(Float32(0.5) * Float32(Float32(1.0) / Float32(Float32(1.0) + exp(Float32(x_m / s))))) / s) end
x_m = abs(x); function tmp = code(x_m, s) tmp = (single(0.5) * (single(1.0) / (single(1.0) + exp((x_m / s))))) / s; end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{0.5 \cdot \frac{1}{1 + e^{\frac{x\_m}{s}}}}{s}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
associate-*r/58.0%
+-commutative58.0%
Applied egg-rr58.0%
Final simplification58.0%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ (/ 0.5 s) (+ 1.0 (exp (/ x_m s)))))
x_m = fabs(x);
float code(float x_m, float s) {
return (0.5f / s) / (1.0f + expf((x_m / s)));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = (0.5e0 / s) / (1.0e0 + exp((x_m / s)))
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(Float32(0.5) / s) / Float32(Float32(1.0) + exp(Float32(x_m / s)))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = (single(0.5) / s) / (single(1.0) + exp((x_m / s))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{\frac{0.5}{s}}{1 + e^{\frac{x\_m}{s}}}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
Taylor expanded in x around inf 58.0%
associate-/r*58.0%
Simplified58.0%
Final simplification58.0%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (if (<= x_m 0.00019999999494757503) (/ 0.25 s) (* (/ 0.5 s) (/ 1.0 (/ x_m s)))))
x_m = fabs(x);
float code(float x_m, float s) {
float tmp;
if (x_m <= 0.00019999999494757503f) {
tmp = 0.25f / s;
} else {
tmp = (0.5f / s) * (1.0f / (x_m / s));
}
return tmp;
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
real(4) :: tmp
if (x_m <= 0.00019999999494757503e0) then
tmp = 0.25e0 / s
else
tmp = (0.5e0 / s) * (1.0e0 / (x_m / s))
end if
code = tmp
end function
x_m = abs(x) function code(x_m, s) tmp = Float32(0.0) if (x_m <= Float32(0.00019999999494757503)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(Float32(0.5) / s) * Float32(Float32(1.0) / Float32(x_m / s))); end return tmp end
x_m = abs(x); function tmp_2 = code(x_m, s) tmp = single(0.0); if (x_m <= single(0.00019999999494757503)) tmp = single(0.25) / s; else tmp = (single(0.5) / s) * (single(1.0) / (x_m / s)); end tmp_2 = tmp; end
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
\mathbf{if}\;x\_m \leq 0.00019999999494757503:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{s} \cdot \frac{1}{\frac{x\_m}{s}}\\
\end{array}
\end{array}
if x < 1.99999995e-4Initial program 99.1%
fabs-neg99.1%
distribute-frac-neg99.1%
distribute-frac-neg299.1%
fabs-neg99.1%
*-commutative99.1%
fabs-neg99.1%
+-commutative99.1%
fabs-neg99.1%
Simplified99.0%
Taylor expanded in s around inf 35.5%
if 1.99999995e-4 < x Initial program 100.0%
fabs-neg100.0%
distribute-frac-neg100.0%
distribute-frac-neg2100.0%
fabs-neg100.0%
*-commutative100.0%
fabs-neg100.0%
+-commutative100.0%
fabs-neg100.0%
Simplified100.0%
Applied egg-rr1.6%
Taylor expanded in x around 0 98.9%
Taylor expanded in x around 0 43.7%
Taylor expanded in x around inf 43.6%
Final simplification37.5%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (if (<= x_m 0.00019999999494757503) (/ 0.25 s) (* (/ 0.5 s) (/ s x_m))))
x_m = fabs(x);
float code(float x_m, float s) {
float tmp;
if (x_m <= 0.00019999999494757503f) {
tmp = 0.25f / s;
} else {
tmp = (0.5f / s) * (s / x_m);
}
return tmp;
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
real(4) :: tmp
if (x_m <= 0.00019999999494757503e0) then
tmp = 0.25e0 / s
else
tmp = (0.5e0 / s) * (s / x_m)
end if
code = tmp
end function
x_m = abs(x) function code(x_m, s) tmp = Float32(0.0) if (x_m <= Float32(0.00019999999494757503)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(Float32(0.5) / s) * Float32(s / x_m)); end return tmp end
x_m = abs(x); function tmp_2 = code(x_m, s) tmp = single(0.0); if (x_m <= single(0.00019999999494757503)) tmp = single(0.25) / s; else tmp = (single(0.5) / s) * (s / x_m); end tmp_2 = tmp; end
\begin{array}{l}
x_m = \left|x\right|
\\
\begin{array}{l}
\mathbf{if}\;x\_m \leq 0.00019999999494757503:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.5}{s} \cdot \frac{s}{x\_m}\\
\end{array}
\end{array}
if x < 1.99999995e-4Initial program 99.1%
fabs-neg99.1%
distribute-frac-neg99.1%
distribute-frac-neg299.1%
fabs-neg99.1%
*-commutative99.1%
fabs-neg99.1%
+-commutative99.1%
fabs-neg99.1%
Simplified99.0%
Taylor expanded in s around inf 35.5%
if 1.99999995e-4 < x Initial program 100.0%
fabs-neg100.0%
distribute-frac-neg100.0%
distribute-frac-neg2100.0%
fabs-neg100.0%
*-commutative100.0%
fabs-neg100.0%
+-commutative100.0%
fabs-neg100.0%
Simplified100.0%
Applied egg-rr1.6%
Taylor expanded in x around 0 98.9%
Taylor expanded in x around 0 43.7%
Taylor expanded in x around inf 36.3%
Final simplification35.7%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (* (/ 0.5 s) (/ 1.0 (+ (/ x_m s) 2.0))))
x_m = fabs(x);
float code(float x_m, float s) {
return (0.5f / s) * (1.0f / ((x_m / s) + 2.0f));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = (0.5e0 / s) * (1.0e0 / ((x_m / s) + 2.0e0))
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(Float32(0.5) / s) * Float32(Float32(1.0) / Float32(Float32(x_m / s) + Float32(2.0)))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = (single(0.5) / s) * (single(1.0) / ((x_m / s) + single(2.0))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{0.5}{s} \cdot \frac{1}{\frac{x\_m}{s} + 2}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
Taylor expanded in x around 0 47.3%
Final simplification47.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ 1.0 (/ (+ (/ x_m s) 2.0) (/ 0.5 s))))
x_m = fabs(x);
float code(float x_m, float s) {
return 1.0f / (((x_m / s) + 2.0f) / (0.5f / s));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = 1.0e0 / (((x_m / s) + 2.0e0) / (0.5e0 / s))
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(1.0) / Float32(Float32(Float32(x_m / s) + Float32(2.0)) / Float32(Float32(0.5) / s))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = single(1.0) / (((x_m / s) + single(2.0)) / (single(0.5) / s)); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{1}{\frac{\frac{x\_m}{s} + 2}{\frac{0.5}{s}}}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
Taylor expanded in x around 0 47.3%
associate-*l/47.3%
*-un-lft-identity47.3%
clear-num47.3%
+-commutative47.3%
Applied egg-rr47.3%
Final simplification47.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ 0.5 (* s (+ (/ x_m s) 2.0))))
x_m = fabs(x);
float code(float x_m, float s) {
return 0.5f / (s * ((x_m / s) + 2.0f));
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = 0.5e0 / (s * ((x_m / s) + 2.0e0))
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(0.5) / Float32(s * Float32(Float32(x_m / s) + Float32(2.0)))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = single(0.5) / (s * ((x_m / s) + single(2.0))); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{0.5}{s \cdot \left(\frac{x\_m}{s} + 2\right)}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
Taylor expanded in x around 0 47.3%
frac-times47.3%
metadata-eval47.3%
+-commutative47.3%
Applied egg-rr47.3%
Final simplification47.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ (/ 0.5 s) (+ (/ x_m s) 2.0)))
x_m = fabs(x);
float code(float x_m, float s) {
return (0.5f / s) / ((x_m / s) + 2.0f);
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = (0.5e0 / s) / ((x_m / s) + 2.0e0)
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(Float32(0.5) / s) / Float32(Float32(x_m / s) + Float32(2.0))) end
x_m = abs(x); function tmp = code(x_m, s) tmp = (single(0.5) / s) / ((x_m / s) + single(2.0)); end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{\frac{0.5}{s}}{\frac{x\_m}{s} + 2}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Applied egg-rr66.1%
Taylor expanded in x around 0 58.0%
Taylor expanded in x around 0 47.3%
associate-*l/47.3%
*-un-lft-identity47.3%
+-commutative47.3%
Applied egg-rr47.3%
Final simplification47.3%
x_m = (fabs.f32 x) (FPCore (x_m s) :precision binary32 (/ 0.25 s))
x_m = fabs(x);
float code(float x_m, float s) {
return 0.25f / s;
}
x_m = abs(x)
real(4) function code(x_m, s)
real(4), intent (in) :: x_m
real(4), intent (in) :: s
code = 0.25e0 / s
end function
x_m = abs(x) function code(x_m, s) return Float32(Float32(0.25) / s) end
x_m = abs(x); function tmp = code(x_m, s) tmp = single(0.25) / s; end
\begin{array}{l}
x_m = \left|x\right|
\\
\frac{0.25}{s}
\end{array}
Initial program 99.3%
fabs-neg99.3%
distribute-frac-neg99.3%
distribute-frac-neg299.3%
fabs-neg99.3%
*-commutative99.3%
fabs-neg99.3%
+-commutative99.3%
fabs-neg99.3%
Simplified99.3%
Taylor expanded in s around inf 28.1%
Final simplification28.1%
herbie shell --seed 2024053
(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))))))