\frac{\cosh x \cdot \frac{y}{x}}{z}
\begin{array}{l}
\mathbf{if}\;z \leq -40841048426.870995:\\
\;\;\;\;\cosh x \cdot \frac{y}{z \cdot x}\\
\mathbf{elif}\;z \leq 7.323266843150271 \cdot 10^{-31}:\\
\;\;\;\;\begin{array}{l}
t_0 := \sqrt{\cosh x}\\
\frac{t_0 \cdot \left(t_0 \cdot \frac{y}{x}\right)}{z}
\end{array}\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \frac{y \cdot \left(e^{x} + \frac{1}{e^{x}}\right)}{z \cdot x}\\
\end{array}
(FPCore (x y z) :precision binary64 (/ (* (cosh x) (/ y x)) z))
(FPCore (x y z)
:precision binary64
(if (<= z -40841048426.870995)
(* (cosh x) (/ y (* z x)))
(if (<= z 7.323266843150271e-31)
(let* ((t_0 (sqrt (cosh x)))) (/ (* t_0 (* t_0 (/ y x))) z))
(* 0.5 (/ (* y (+ (exp x) (/ 1.0 (exp x)))) (* z x))))))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 (z <= -40841048426.870995) {
tmp = cosh(x) * (y / (z * x));
} else if (z <= 7.323266843150271e-31) {
double t_0 = sqrt(cosh(x));
tmp = (t_0 * (t_0 * (y / x))) / z;
} else {
tmp = 0.5 * ((y * (exp(x) + (1.0 / exp(x)))) / (z * x));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.6 |
|---|---|
| Target | 0.4 |
| Herbie | 0.3 |
if z < -40841048426.8709946Initial program 12.1
Applied *-un-lft-identity_binary6412.1
Applied times-frac_binary6412.1
Simplified12.1
Simplified0.3
if -40841048426.8709946 < z < 7.3232668431502707e-31Initial program 0.3
Applied add-sqr-sqrt_binary640.3
Applied associate-*l*_binary640.3
Simplified0.3
if 7.3232668431502707e-31 < z Initial program 10.6
Taylor expanded in y around inf 0.4
Final simplification0.3
herbie shell --seed 2021275
(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.038530535935153e-39) (/ (/ (* (cosh x) y) x) z) (* (/ (/ y z) x) (cosh x))))
(/ (* (cosh x) (/ y x)) z))