
(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 21 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}
(FPCore (x s) :precision binary32 (let* ((t_0 (exp (/ (- 0.0 (fabs x)) s)))) (/ t_0 (* (* s (+ 1.0 (/ 1.0 (exp (/ (fabs x) s))))) (+ t_0 1.0)))))
float code(float x, float s) {
float t_0 = expf(((0.0f - fabsf(x)) / s));
return t_0 / ((s * (1.0f + (1.0f / expf((fabsf(x) / s))))) * (t_0 + 1.0f));
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: t_0
t_0 = exp(((0.0e0 - abs(x)) / s))
code = t_0 / ((s * (1.0e0 + (1.0e0 / exp((abs(x) / s))))) * (t_0 + 1.0e0))
end function
function code(x, s) t_0 = exp(Float32(Float32(Float32(0.0) - abs(x)) / s)) return Float32(t_0 / Float32(Float32(s * Float32(Float32(1.0) + Float32(Float32(1.0) / exp(Float32(abs(x) / s))))) * Float32(t_0 + Float32(1.0)))) end
function tmp = code(x, s) t_0 = exp(((single(0.0) - abs(x)) / s)); tmp = t_0 / ((s * (single(1.0) + (single(1.0) / exp((abs(x) / s))))) * (t_0 + single(1.0))); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\frac{0 - \left|x\right|}{s}}\\
\frac{t\_0}{\left(s \cdot \left(1 + \frac{1}{e^{\frac{\left|x\right|}{s}}}\right)\right) \cdot \left(t\_0 + 1\right)}
\end{array}
\end{array}
Initial program 99.8%
distribute-frac-negN/A
exp-negN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (x s) :precision binary32 (let* ((t_0 (exp (/ (- 0.0 (fabs x)) s)))) (/ (/ t_0 s) (pow (+ t_0 1.0) 2.0))))
float code(float x, float s) {
float t_0 = expf(((0.0f - fabsf(x)) / s));
return (t_0 / s) / powf((t_0 + 1.0f), 2.0f);
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: t_0
t_0 = exp(((0.0e0 - abs(x)) / s))
code = (t_0 / s) / ((t_0 + 1.0e0) ** 2.0e0)
end function
function code(x, s) t_0 = exp(Float32(Float32(Float32(0.0) - abs(x)) / s)) return Float32(Float32(t_0 / s) / (Float32(t_0 + Float32(1.0)) ^ Float32(2.0))) end
function tmp = code(x, s) t_0 = exp(((single(0.0) - abs(x)) / s)); tmp = (t_0 / s) / ((t_0 + single(1.0)) ^ single(2.0)); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\frac{0 - \left|x\right|}{s}}\\
\frac{\frac{t\_0}{s}}{{\left(t\_0 + 1\right)}^{2}}
\end{array}
\end{array}
Initial program 99.8%
Taylor expanded in x around 0
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
neg-mul-1N/A
distribute-frac-neg2N/A
neg-mul-1N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32N/A
neg-mul-1N/A
neg-lowering-neg.f32N/A
pow-lowering-pow.f32N/A
Simplified99.8%
Final simplification99.8%
(FPCore (x s) :precision binary32 (/ (pow (sqrt (/ 1.0 (+ 2.0 (* 2.0 (cosh (/ (fabs x) s)))))) 2.0) s))
float code(float x, float s) {
return powf(sqrtf((1.0f / (2.0f + (2.0f * coshf((fabsf(x) / s)))))), 2.0f) / s;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = (sqrt((1.0e0 / (2.0e0 + (2.0e0 * cosh((abs(x) / s)))))) ** 2.0e0) / s
end function
function code(x, s) return Float32((sqrt(Float32(Float32(1.0) / Float32(Float32(2.0) + Float32(Float32(2.0) * cosh(Float32(abs(x) / s)))))) ^ Float32(2.0)) / s) end
function tmp = code(x, s) tmp = (sqrt((single(1.0) / (single(2.0) + (single(2.0) * cosh((abs(x) / s)))))) ^ single(2.0)) / s; end
\begin{array}{l}
\\
\frac{{\left(\sqrt{\frac{1}{2 + 2 \cdot \cosh \left(\frac{\left|x\right|}{s}\right)}}\right)}^{2}}{s}
\end{array}
Initial program 99.8%
Simplified99.5%
inv-powN/A
sqr-powN/A
pow2N/A
pow-lowering-pow.f32N/A
Applied egg-rr99.5%
unpow1N/A
sqr-powN/A
pow-prod-downN/A
unpow2N/A
pow-lowering-pow.f32N/A
pow-powN/A
metadata-evalN/A
inv-powN/A
/-lowering-/.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32N/A
Applied egg-rr99.5%
unpow1/2N/A
sqrt-lowering-sqrt.f32N/A
/-lowering-/.f32N/A
+-lowering-+.f32N/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.5%
Applied egg-rr99.5%
(FPCore (x s) :precision binary32 (/ (pow (pow (+ 2.0 (* 2.0 (cosh (/ (fabs x) s)))) -0.5) 2.0) s))
float code(float x, float s) {
return powf(powf((2.0f + (2.0f * coshf((fabsf(x) / s)))), -0.5f), 2.0f) / s;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = (((2.0e0 + (2.0e0 * cosh((abs(x) / s)))) ** (-0.5e0)) ** 2.0e0) / s
end function
function code(x, s) return Float32(((Float32(Float32(2.0) + Float32(Float32(2.0) * cosh(Float32(abs(x) / s)))) ^ Float32(-0.5)) ^ Float32(2.0)) / s) end
function tmp = code(x, s) tmp = (((single(2.0) + (single(2.0) * cosh((abs(x) / s)))) ^ single(-0.5)) ^ single(2.0)) / s; end
\begin{array}{l}
\\
\frac{{\left({\left(2 + 2 \cdot \cosh \left(\frac{\left|x\right|}{s}\right)\right)}^{-0.5}\right)}^{2}}{s}
\end{array}
Initial program 99.8%
Simplified99.5%
inv-powN/A
sqr-powN/A
pow2N/A
pow-lowering-pow.f32N/A
Applied egg-rr99.5%
(FPCore (x s) :precision binary32 (/ 1.0 (* s (+ 2.0 (* 2.0 (cosh (/ (fabs x) s)))))))
float code(float x, float s) {
return 1.0f / (s * (2.0f + (2.0f * coshf((fabsf(x) / s)))));
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = 1.0e0 / (s * (2.0e0 + (2.0e0 * cosh((abs(x) / s)))))
end function
function code(x, s) return Float32(Float32(1.0) / Float32(s * Float32(Float32(2.0) + Float32(Float32(2.0) * cosh(Float32(abs(x) / s)))))) end
function tmp = code(x, s) tmp = single(1.0) / (s * (single(2.0) + (single(2.0) * cosh((abs(x) / s))))); end
\begin{array}{l}
\\
\frac{1}{s \cdot \left(2 + 2 \cdot \cosh \left(\frac{\left|x\right|}{s}\right)\right)}
\end{array}
Initial program 99.8%
Simplified99.5%
associate-/l/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.5%
Applied egg-rr99.5%
(FPCore (x s) :precision binary32 (/ (exp (* (fabs x) (/ -1.0 s))) (* s 4.0)))
float code(float x, float s) {
return expf((fabsf(x) * (-1.0f / s))) / (s * 4.0f);
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = exp((abs(x) * ((-1.0e0) / s))) / (s * 4.0e0)
end function
function code(x, s) return Float32(exp(Float32(abs(x) * Float32(Float32(-1.0) / s))) / Float32(s * Float32(4.0))) end
function tmp = code(x, s) tmp = exp((abs(x) * (single(-1.0) / s))) / (s * single(4.0)); end
\begin{array}{l}
\\
\frac{e^{\left|x\right| \cdot \frac{-1}{s}}}{s \cdot 4}
\end{array}
Initial program 99.8%
Taylor expanded in s around inf
*-commutativeN/A
*-lowering-*.f3295.7%
Simplified95.7%
distribute-frac-negN/A
clear-numN/A
associate-/r/N/A
distribute-lft-neg-inN/A
*-lowering-*.f32N/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3295.7%
Applied egg-rr95.7%
Final simplification95.7%
(FPCore (x s) :precision binary32 (/ (exp (/ (- 0.0 (fabs x)) s)) (* s 4.0)))
float code(float x, float s) {
return expf(((0.0f - fabsf(x)) / s)) / (s * 4.0f);
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = exp(((0.0e0 - abs(x)) / s)) / (s * 4.0e0)
end function
function code(x, s) return Float32(exp(Float32(Float32(Float32(0.0) - abs(x)) / s)) / Float32(s * Float32(4.0))) end
function tmp = code(x, s) tmp = exp(((single(0.0) - abs(x)) / s)) / (s * single(4.0)); end
\begin{array}{l}
\\
\frac{e^{\frac{0 - \left|x\right|}{s}}}{s \cdot 4}
\end{array}
Initial program 99.8%
Taylor expanded in s around inf
*-commutativeN/A
*-lowering-*.f3295.7%
Simplified95.7%
Final simplification95.7%
(FPCore (x s) :precision binary32 (/ (/ 0.25 (exp (/ (fabs x) s))) s))
float code(float x, float s) {
return (0.25f / expf((fabsf(x) / s))) / s;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = (0.25e0 / exp((abs(x) / s))) / s
end function
function code(x, s) return Float32(Float32(Float32(0.25) / exp(Float32(abs(x) / s))) / s) end
function tmp = code(x, s) tmp = (single(0.25) / exp((abs(x) / s))) / s; end
\begin{array}{l}
\\
\frac{\frac{0.25}{e^{\frac{\left|x\right|}{s}}}}{s}
\end{array}
Initial program 99.8%
Taylor expanded in s around inf
*-commutativeN/A
*-lowering-*.f3295.7%
Simplified95.7%
Taylor expanded in x around 0
associate-*r/N/A
/-lowering-/.f32N/A
mul-1-negN/A
rec-expN/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3295.6%
Simplified95.6%
(FPCore (x s)
:precision binary32
(if (<= x 1.999999936531045e-21)
(/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s))))
(/
(/ 1.0 s)
(+
2.0
(*
2.0
(+
1.0
(*
(* x x)
(+
(/ (* (* x x) 0.041666666666666664) (* s (* s (* s s))))
(/ 0.5 (* s s))))))))))
float code(float x, float s) {
float tmp;
if (x <= 1.999999936531045e-21f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else {
tmp = (1.0f / s) / (2.0f + (2.0f * (1.0f + ((x * x) * ((((x * x) * 0.041666666666666664f) / (s * (s * (s * s)))) + (0.5f / (s * s)))))));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 1.999999936531045e-21) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else
tmp = (1.0e0 / s) / (2.0e0 + (2.0e0 * (1.0e0 + ((x * x) * ((((x * x) * 0.041666666666666664e0) / (s * (s * (s * s)))) + (0.5e0 / (s * s)))))))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(1.999999936531045e-21)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); else tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(2.0) + Float32(Float32(2.0) * Float32(Float32(1.0) + Float32(Float32(x * x) * Float32(Float32(Float32(Float32(x * x) * Float32(0.041666666666666664)) / Float32(s * Float32(s * Float32(s * s)))) + Float32(Float32(0.5) / Float32(s * s)))))))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(1.999999936531045e-21)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); else tmp = (single(1.0) / s) / (single(2.0) + (single(2.0) * (single(1.0) + ((x * x) * ((((x * x) * single(0.041666666666666664)) / (s * (s * (s * s)))) + (single(0.5) / (s * s))))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.999999936531045 \cdot 10^{-21}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{s}}{2 + 2 \cdot \left(1 + \left(x \cdot x\right) \cdot \left(\frac{\left(x \cdot x\right) \cdot 0.041666666666666664}{s \cdot \left(s \cdot \left(s \cdot s\right)\right)} + \frac{0.5}{s \cdot s}\right)\right)}\\
\end{array}
\end{array}
if x < 1.9999999e-21Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.4%
Simplified74.4%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.6%
Applied egg-rr76.6%
if 1.9999999e-21 < x Initial program 99.7%
Simplified99.1%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.1%
Applied egg-rr99.1%
Taylor expanded in s around -inf
Simplified33.9%
Taylor expanded in x around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified93.9%
Final simplification84.3%
(FPCore (x s)
:precision binary32
(if (<= x 5.000000136226006e-28)
(/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s))))
(if (<= x 4.999999980020986e-12)
(/ (/ 1.0 (+ 4.0 (* x (/ x (* s s))))) s)
(/
(/ 1.0 s)
(+
2.0
(*
2.0
(*
(* (* x x) (* x x))
(/ 0.041666666666666664 (* s (* s (* s s)))))))))))
float code(float x, float s) {
float tmp;
if (x <= 5.000000136226006e-28f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else if (x <= 4.999999980020986e-12f) {
tmp = (1.0f / (4.0f + (x * (x / (s * s))))) / s;
} else {
tmp = (1.0f / s) / (2.0f + (2.0f * (((x * x) * (x * x)) * (0.041666666666666664f / (s * (s * (s * s)))))));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 5.000000136226006e-28) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else if (x <= 4.999999980020986e-12) then
tmp = (1.0e0 / (4.0e0 + (x * (x / (s * s))))) / s
else
tmp = (1.0e0 / s) / (2.0e0 + (2.0e0 * (((x * x) * (x * x)) * (0.041666666666666664e0 / (s * (s * (s * s)))))))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(5.000000136226006e-28)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); elseif (x <= Float32(4.999999980020986e-12)) tmp = Float32(Float32(Float32(1.0) / Float32(Float32(4.0) + Float32(x * Float32(x / Float32(s * s))))) / s); else tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(2.0) + Float32(Float32(2.0) * Float32(Float32(Float32(x * x) * Float32(x * x)) * Float32(Float32(0.041666666666666664) / Float32(s * Float32(s * Float32(s * s)))))))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(5.000000136226006e-28)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); elseif (x <= single(4.999999980020986e-12)) tmp = (single(1.0) / (single(4.0) + (x * (x / (s * s))))) / s; else tmp = (single(1.0) / s) / (single(2.0) + (single(2.0) * (((x * x) * (x * x)) * (single(0.041666666666666664) / (s * (s * (s * s))))))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.000000136226006 \cdot 10^{-28}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{elif}\;x \leq 4.999999980020986 \cdot 10^{-12}:\\
\;\;\;\;\frac{\frac{1}{4 + x \cdot \frac{x}{s \cdot s}}}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{s}}{2 + 2 \cdot \left(\left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right) \cdot \frac{0.041666666666666664}{s \cdot \left(s \cdot \left(s \cdot s\right)\right)}\right)}\\
\end{array}
\end{array}
if x < 5.00000014e-28Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.5%
Simplified74.5%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.9%
Applied egg-rr76.9%
if 5.00000014e-28 < x < 4.99999998e-12Initial program 99.9%
Simplified100.0%
associate-+r+N/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
unpow2N/A
sqr-absN/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3288.8%
Simplified88.8%
if 4.99999998e-12 < x Initial program 99.7%
Simplified98.9%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.0%
Applied egg-rr99.0%
Taylor expanded in s around -inf
Simplified33.0%
Taylor expanded in x around inf
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
*-lowering-*.f32N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f32N/A
cube-multN/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3291.5%
Simplified91.5%
Final simplification83.7%
(FPCore (x s)
:precision binary32
(if (<= x 5.000000136226006e-28)
(/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s))))
(if (<= x 1.4999999397961872e-13)
(/ (/ 1.0 (+ 4.0 (* x (/ x (* s s))))) s)
(/ (/ (* s (* (* s s) -4.0)) (* x x)) (* x x)))))
float code(float x, float s) {
float tmp;
if (x <= 5.000000136226006e-28f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else if (x <= 1.4999999397961872e-13f) {
tmp = (1.0f / (4.0f + (x * (x / (s * s))))) / s;
} else {
tmp = ((s * ((s * s) * -4.0f)) / (x * x)) / (x * x);
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 5.000000136226006e-28) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else if (x <= 1.4999999397961872e-13) then
tmp = (1.0e0 / (4.0e0 + (x * (x / (s * s))))) / s
else
tmp = ((s * ((s * s) * (-4.0e0))) / (x * x)) / (x * x)
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(5.000000136226006e-28)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); elseif (x <= Float32(1.4999999397961872e-13)) tmp = Float32(Float32(Float32(1.0) / Float32(Float32(4.0) + Float32(x * Float32(x / Float32(s * s))))) / s); else tmp = Float32(Float32(Float32(s * Float32(Float32(s * s) * Float32(-4.0))) / Float32(x * x)) / Float32(x * x)); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(5.000000136226006e-28)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); elseif (x <= single(1.4999999397961872e-13)) tmp = (single(1.0) / (single(4.0) + (x * (x / (s * s))))) / s; else tmp = ((s * ((s * s) * single(-4.0))) / (x * x)) / (x * x); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.000000136226006 \cdot 10^{-28}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{elif}\;x \leq 1.4999999397961872 \cdot 10^{-13}:\\
\;\;\;\;\frac{\frac{1}{4 + x \cdot \frac{x}{s \cdot s}}}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{s \cdot \left(\left(s \cdot s\right) \cdot -4\right)}{x \cdot x}}{x \cdot x}\\
\end{array}
\end{array}
if x < 5.00000014e-28Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.5%
Simplified74.5%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.9%
Applied egg-rr76.9%
if 5.00000014e-28 < x < 1.49999994e-13Initial program 99.9%
Simplified100.0%
associate-+r+N/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
unpow2N/A
sqr-absN/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3295.6%
Simplified95.6%
if 1.49999994e-13 < x Initial program 99.7%
Simplified99.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.0%
Applied egg-rr99.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3278.1%
Simplified78.1%
Taylor expanded in x around inf
/-lowering-/.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
cube-multN/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3256.5%
Simplified56.5%
Taylor expanded in s around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3287.0%
Simplified87.0%
Final simplification82.6%
(FPCore (x s)
:precision binary32
(if (<= x 5.000000136226006e-28)
(/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s))))
(if (<= x 4.999999980020986e-12)
(/ (/ 1.0 (+ 4.0 (* x (/ x (* s s))))) s)
(/ (* s (* (* s s) -4.0)) (* (* x x) (* x x))))))
float code(float x, float s) {
float tmp;
if (x <= 5.000000136226006e-28f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else if (x <= 4.999999980020986e-12f) {
tmp = (1.0f / (4.0f + (x * (x / (s * s))))) / s;
} else {
tmp = (s * ((s * s) * -4.0f)) / ((x * x) * (x * x));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 5.000000136226006e-28) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else if (x <= 4.999999980020986e-12) then
tmp = (1.0e0 / (4.0e0 + (x * (x / (s * s))))) / s
else
tmp = (s * ((s * s) * (-4.0e0))) / ((x * x) * (x * x))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(5.000000136226006e-28)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); elseif (x <= Float32(4.999999980020986e-12)) tmp = Float32(Float32(Float32(1.0) / Float32(Float32(4.0) + Float32(x * Float32(x / Float32(s * s))))) / s); else tmp = Float32(Float32(s * Float32(Float32(s * s) * Float32(-4.0))) / Float32(Float32(x * x) * Float32(x * x))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(5.000000136226006e-28)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); elseif (x <= single(4.999999980020986e-12)) tmp = (single(1.0) / (single(4.0) + (x * (x / (s * s))))) / s; else tmp = (s * ((s * s) * single(-4.0))) / ((x * x) * (x * x)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.000000136226006 \cdot 10^{-28}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{elif}\;x \leq 4.999999980020986 \cdot 10^{-12}:\\
\;\;\;\;\frac{\frac{1}{4 + x \cdot \frac{x}{s \cdot s}}}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{s \cdot \left(\left(s \cdot s\right) \cdot -4\right)}{\left(x \cdot x\right) \cdot \left(x \cdot x\right)}\\
\end{array}
\end{array}
if x < 5.00000014e-28Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.5%
Simplified74.5%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.9%
Applied egg-rr76.9%
if 5.00000014e-28 < x < 4.99999998e-12Initial program 99.9%
Simplified100.0%
associate-+r+N/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
unpow2N/A
sqr-absN/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3288.8%
Simplified88.8%
if 4.99999998e-12 < x Initial program 99.7%
Simplified98.9%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.0%
Applied egg-rr99.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3279.4%
Simplified79.4%
Taylor expanded in x around inf
/-lowering-/.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
cube-multN/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3258.1%
Simplified58.1%
Taylor expanded in s around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
cube-multN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3287.6%
Simplified87.6%
Final simplification82.2%
(FPCore (x s) :precision binary32 (if (<= x 5.000000136226006e-28) (/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s)))) (/ (/ 1.0 (+ 4.0 (* x (/ x (* s s))))) s)))
float code(float x, float s) {
float tmp;
if (x <= 5.000000136226006e-28f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else {
tmp = (1.0f / (4.0f + (x * (x / (s * s))))) / s;
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 5.000000136226006e-28) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else
tmp = (1.0e0 / (4.0e0 + (x * (x / (s * s))))) / s
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(5.000000136226006e-28)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); else tmp = Float32(Float32(Float32(1.0) / Float32(Float32(4.0) + Float32(x * Float32(x / Float32(s * s))))) / s); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(5.000000136226006e-28)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); else tmp = (single(1.0) / (single(4.0) + (x * (x / (s * s))))) / s; end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.000000136226006 \cdot 10^{-28}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{4 + x \cdot \frac{x}{s \cdot s}}}{s}\\
\end{array}
\end{array}
if x < 5.00000014e-28Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.5%
Simplified74.5%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.9%
Applied egg-rr76.9%
if 5.00000014e-28 < x Initial program 99.7%
Simplified99.2%
associate-+r+N/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.2%
Applied egg-rr99.2%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
unpow2N/A
sqr-absN/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3281.4%
Simplified81.4%
Final simplification79.1%
(FPCore (x s) :precision binary32 (if (<= x 9.999999998199587e-24) (/ (/ 1.0 s) (+ 4.0 (* (/ x s) (/ x s)))) (/ 1.0 (* s (+ 4.0 (/ (* x x) (* s s)))))))
float code(float x, float s) {
float tmp;
if (x <= 9.999999998199587e-24f) {
tmp = (1.0f / s) / (4.0f + ((x / s) * (x / s)));
} else {
tmp = 1.0f / (s * (4.0f + ((x * x) / (s * s))));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 9.999999998199587e-24) then
tmp = (1.0e0 / s) / (4.0e0 + ((x / s) * (x / s)))
else
tmp = 1.0e0 / (s * (4.0e0 + ((x * x) / (s * s))))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(9.999999998199587e-24)) tmp = Float32(Float32(Float32(1.0) / s) / Float32(Float32(4.0) + Float32(Float32(x / s) * Float32(x / s)))); else tmp = Float32(Float32(1.0) / Float32(s * Float32(Float32(4.0) + Float32(Float32(x * x) / Float32(s * s))))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(9.999999998199587e-24)) tmp = (single(1.0) / s) / (single(4.0) + ((x / s) * (x / s))); else tmp = single(1.0) / (s * (single(4.0) + ((x * x) / (s * s)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.999999998199587 \cdot 10^{-24}:\\
\;\;\;\;\frac{\frac{1}{s}}{4 + \frac{x}{s} \cdot \frac{x}{s}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{s \cdot \left(4 + \frac{x \cdot x}{s \cdot s}\right)}\\
\end{array}
\end{array}
if x < 1e-23Initial program 99.8%
Simplified99.8%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.8%
Applied egg-rr99.8%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3274.0%
Simplified74.0%
times-fracN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3276.4%
Applied egg-rr76.4%
if 1e-23 < x Initial program 99.7%
Simplified99.1%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.1%
Applied egg-rr99.1%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3279.9%
Simplified79.9%
associate-/l/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3279.9%
Applied egg-rr79.9%
Final simplification78.0%
(FPCore (x s) :precision binary32 (if (<= x 9.999999998199587e-24) (/ 0.25 s) (/ 1.0 (* s (+ 4.0 (/ (* x x) (* s s)))))))
float code(float x, float s) {
float tmp;
if (x <= 9.999999998199587e-24f) {
tmp = 0.25f / s;
} else {
tmp = 1.0f / (s * (4.0f + ((x * x) / (s * s))));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 9.999999998199587e-24) then
tmp = 0.25e0 / s
else
tmp = 1.0e0 / (s * (4.0e0 + ((x * x) / (s * s))))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(9.999999998199587e-24)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(1.0) / Float32(s * Float32(Float32(4.0) + Float32(Float32(x * x) / Float32(s * s))))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(9.999999998199587e-24)) tmp = single(0.25) / s; else tmp = single(1.0) / (s * (single(4.0) + ((x * x) / (s * s)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.999999998199587 \cdot 10^{-24}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{s \cdot \left(4 + \frac{x \cdot x}{s \cdot s}\right)}\\
\end{array}
\end{array}
if x < 1e-23Initial program 99.8%
Taylor expanded in s around inf
/-lowering-/.f3238.1%
Simplified38.1%
if 1e-23 < x Initial program 99.7%
Simplified99.1%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.1%
Applied egg-rr99.1%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3279.9%
Simplified79.9%
associate-/l/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3279.9%
Applied egg-rr79.9%
Final simplification56.9%
(FPCore (x s) :precision binary32 (if (<= x 1.199999961918627e-19) (/ 0.25 s) (/ (/ 1.0 s) (* x (/ x (* s s))))))
float code(float x, float s) {
float tmp;
if (x <= 1.199999961918627e-19f) {
tmp = 0.25f / s;
} else {
tmp = (1.0f / s) / (x * (x / (s * s)));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 1.199999961918627e-19) then
tmp = 0.25e0 / s
else
tmp = (1.0e0 / s) / (x * (x / (s * s)))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(1.199999961918627e-19)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(Float32(1.0) / s) / Float32(x * Float32(x / Float32(s * s)))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(1.199999961918627e-19)) tmp = single(0.25) / s; else tmp = (single(1.0) / s) / (x * (x / (s * s))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.199999961918627 \cdot 10^{-19}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{s}}{x \cdot \frac{x}{s \cdot s}}\\
\end{array}
\end{array}
if x < 1.19999996e-19Initial program 99.8%
Taylor expanded in s around inf
/-lowering-/.f3239.7%
Simplified39.7%
if 1.19999996e-19 < x Initial program 99.7%
Simplified99.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f3299.1%
Applied egg-rr99.1%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3279.7%
Simplified79.7%
Taylor expanded in x around inf
unpow2N/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3275.1%
Simplified75.1%
(FPCore (x s) :precision binary32 (if (<= x 4.0000000467443897e-7) (/ 0.25 s) (/ 1.0 (/ x (/ s x)))))
float code(float x, float s) {
float tmp;
if (x <= 4.0000000467443897e-7f) {
tmp = 0.25f / s;
} else {
tmp = 1.0f / (x / (s / x));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 4.0000000467443897e-7) then
tmp = 0.25e0 / s
else
tmp = 1.0e0 / (x / (s / x))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(4.0000000467443897e-7)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(1.0) / Float32(x / Float32(s / x))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(4.0000000467443897e-7)) tmp = single(0.25) / s; else tmp = single(1.0) / (x / (s / x)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.0000000467443897 \cdot 10^{-7}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{x}{\frac{s}{x}}}\\
\end{array}
\end{array}
if x < 4.00000005e-7Initial program 99.6%
Taylor expanded in s around inf
/-lowering-/.f3238.3%
Simplified38.3%
if 4.00000005e-7 < x Initial program 100.0%
Simplified100.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3281.1%
Simplified81.1%
Taylor expanded in s around 0
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3264.7%
Simplified64.7%
associate-/r*N/A
clear-numN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3268.5%
Applied egg-rr68.5%
(FPCore (x s) :precision binary32 (if (<= x 4.0000000467443897e-7) (/ 0.25 s) (* (/ s x) (/ 1.0 x))))
float code(float x, float s) {
float tmp;
if (x <= 4.0000000467443897e-7f) {
tmp = 0.25f / s;
} else {
tmp = (s / x) * (1.0f / x);
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 4.0000000467443897e-7) then
tmp = 0.25e0 / s
else
tmp = (s / x) * (1.0e0 / x)
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(4.0000000467443897e-7)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(Float32(s / x) * Float32(Float32(1.0) / x)); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(4.0000000467443897e-7)) tmp = single(0.25) / s; else tmp = (s / x) * (single(1.0) / x); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.0000000467443897 \cdot 10^{-7}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{s}{x} \cdot \frac{1}{x}\\
\end{array}
\end{array}
if x < 4.00000005e-7Initial program 99.6%
Taylor expanded in s around inf
/-lowering-/.f3238.3%
Simplified38.3%
if 4.00000005e-7 < x Initial program 100.0%
Simplified100.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3281.1%
Simplified81.1%
Taylor expanded in s around 0
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3264.7%
Simplified64.7%
associate-/r*N/A
div-invN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3264.7%
Applied egg-rr64.7%
(FPCore (x s) :precision binary32 (if (<= x 4.0000000467443897e-7) (/ 0.25 s) (* s (/ 1.0 (* x x)))))
float code(float x, float s) {
float tmp;
if (x <= 4.0000000467443897e-7f) {
tmp = 0.25f / s;
} else {
tmp = s * (1.0f / (x * x));
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 4.0000000467443897e-7) then
tmp = 0.25e0 / s
else
tmp = s * (1.0e0 / (x * x))
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(4.0000000467443897e-7)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(s * Float32(Float32(1.0) / Float32(x * x))); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(4.0000000467443897e-7)) tmp = single(0.25) / s; else tmp = s * (single(1.0) / (x * x)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.0000000467443897 \cdot 10^{-7}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;s \cdot \frac{1}{x \cdot x}\\
\end{array}
\end{array}
if x < 4.00000005e-7Initial program 99.6%
Taylor expanded in s around inf
/-lowering-/.f3238.3%
Simplified38.3%
if 4.00000005e-7 < x Initial program 100.0%
Simplified100.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3281.1%
Simplified81.1%
Taylor expanded in s around 0
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3264.7%
Simplified64.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f3264.7%
Applied egg-rr64.7%
Final simplification47.3%
(FPCore (x s) :precision binary32 (if (<= x 4.0000000467443897e-7) (/ 0.25 s) (/ s (* x x))))
float code(float x, float s) {
float tmp;
if (x <= 4.0000000467443897e-7f) {
tmp = 0.25f / s;
} else {
tmp = s / (x * x);
}
return tmp;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
real(4) :: tmp
if (x <= 4.0000000467443897e-7) then
tmp = 0.25e0 / s
else
tmp = s / (x * x)
end if
code = tmp
end function
function code(x, s) tmp = Float32(0.0) if (x <= Float32(4.0000000467443897e-7)) tmp = Float32(Float32(0.25) / s); else tmp = Float32(s / Float32(x * x)); end return tmp end
function tmp_2 = code(x, s) tmp = single(0.0); if (x <= single(4.0000000467443897e-7)) tmp = single(0.25) / s; else tmp = s / (x * x); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 4.0000000467443897 \cdot 10^{-7}:\\
\;\;\;\;\frac{0.25}{s}\\
\mathbf{else}:\\
\;\;\;\;\frac{s}{x \cdot x}\\
\end{array}
\end{array}
if x < 4.00000005e-7Initial program 99.6%
Taylor expanded in s around inf
/-lowering-/.f3238.3%
Simplified38.3%
if 4.00000005e-7 < x Initial program 100.0%
Simplified100.0%
associate-/l/N/A
associate-/r*N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f32N/A
distribute-frac-negN/A
+-commutativeN/A
distribute-frac-negN/A
cosh-undefN/A
*-lowering-*.f32N/A
cosh-lowering-cosh.f32N/A
/-lowering-/.f32N/A
fabs-lowering-fabs.f32100.0%
Applied egg-rr100.0%
Taylor expanded in s around inf
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
unpow2N/A
sqr-absN/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3281.1%
Simplified81.1%
Taylor expanded in s around 0
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3264.7%
Simplified64.7%
(FPCore (x s) :precision binary32 (/ 0.25 s))
float code(float x, float s) {
return 0.25f / s;
}
real(4) function code(x, s)
real(4), intent (in) :: x
real(4), intent (in) :: s
code = 0.25e0 / s
end function
function code(x, s) return Float32(Float32(0.25) / s) end
function tmp = code(x, s) tmp = single(0.25) / s; end
\begin{array}{l}
\\
\frac{0.25}{s}
\end{array}
Initial program 99.8%
Taylor expanded in s around inf
/-lowering-/.f3226.9%
Simplified26.9%
herbie shell --seed 2024191
(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))))))