\frac{\cosh x \cdot \frac{y}{x}}{z}\begin{array}{l}
\mathbf{if}\;y \leq -9.905533536424995 \cdot 10^{+35}:\\
\;\;\;\;\frac{y \cdot \left(e^{x} + e^{-x}\right)}{z \cdot \left(x \cdot 2\right)}\\
\mathbf{elif}\;y \leq 741420780451993:\\
\;\;\;\;\frac{\left(\sqrt[3]{\cosh x} \cdot \sqrt[3]{\cosh x}\right) \cdot \left(\sqrt[3]{\cosh x} \cdot \frac{y}{x}\right)}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{y \cdot \frac{e^{x} + e^{-x}}{z}}{x \cdot 2}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (* (cosh x) (/ y x)) z))
(FPCore (x y z)
:precision binary64
(if (<= y -9.905533536424995e+35)
(/ (* y (+ (exp x) (exp (- x)))) (* z (* x 2.0)))
(if (<= y 741420780451993.0)
(/ (* (* (cbrt (cosh x)) (cbrt (cosh x))) (* (cbrt (cosh x)) (/ y x))) z)
(/ (* y (/ (+ (exp x) (exp (- x))) z)) (* x 2.0)))))double code(double x, double y, double z) {
return (cosh(x) * (y / x)) / z;
}
double code(double x, double y, double z) {
double tmp;
if (y <= -9.905533536424995e+35) {
tmp = (y * (exp(x) + exp(-x))) / (z * (x * 2.0));
} else if (y <= 741420780451993.0) {
tmp = ((cbrt(cosh(x)) * cbrt(cosh(x))) * (cbrt(cosh(x)) * (y / x))) / z;
} else {
tmp = (y * ((exp(x) + exp(-x)) / z)) / (x * 2.0);
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 8.1 |
|---|---|
| Target | 0.4 |
| Herbie | 0.4 |
if y < -9.90553353642499542e35Initial program 26.4
rmApplied cosh-def_binary64_1542526.4
Applied frac-times_binary64_1525226.4
Applied associate-/l/_binary64_151910.5
Simplified0.5
if -9.90553353642499542e35 < y < 741420780451993Initial program 0.3
rmApplied add-cube-cbrt_binary64_152740.4
Applied associate-*l*_binary64_151850.4
Simplified0.4
if 741420780451993 < y Initial program 23.3
rmApplied cosh-def_binary64_1542523.3
Applied frac-times_binary64_1525223.3
Applied associate-/l/_binary64_151910.4
Simplified0.4
rmApplied associate-/r*_binary64_151880.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020277
(FPCore (x y z)
:name "Linear.Quaternion:$ctan from linear-1.19.1.3"
:precision binary64
:herbie-target
(if (< y -4.618902267687042e-52) (* (/ (/ y z) x) (cosh x)) (if (< y 1.0385305359351529e-39) (/ (/ (* (cosh x) y) x) z) (* (/ (/ y z) x) (cosh x))))
(/ (* (cosh x) (/ y x)) z))