\frac{\cosh x \cdot \frac{y}{x}}{z}\begin{array}{l}
\mathbf{if}\;z \leq -0.06013142179743951:\\
\;\;\;\;\frac{\left(e^{x} + e^{-x}\right) \cdot y}{z \cdot \left(x \cdot 2\right)}\\
\mathbf{elif}\;z \leq 3.797577207409332 \cdot 10^{+78}:\\
\;\;\;\;\frac{1}{x \cdot \frac{z}{y \cdot \cosh x}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(e^{x} + e^{-x}\right) \cdot y}{z \cdot \left(x \cdot 2\right)}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (* (cosh x) (/ y x)) z))
(FPCore (x y z)
:precision binary64
(if (<= z -0.06013142179743951)
(/ (* (+ (exp x) (exp (- x))) y) (* z (* x 2.0)))
(if (<= z 3.797577207409332e+78)
(/ 1.0 (* x (/ z (* y (cosh x)))))
(/ (* (+ (exp x) (exp (- x))) y) (* 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 (z <= -0.06013142179743951) {
tmp = ((exp(x) + exp(-x)) * y) / (z * (x * 2.0));
} else if (z <= 3.797577207409332e+78) {
tmp = 1.0 / (x * (z / (y * cosh(x))));
} else {
tmp = ((exp(x) + exp(-x)) * y) / (z * (x * 2.0));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.7 |
|---|---|
| Target | 0.5 |
| Herbie | 0.6 |
if z < -0.060131421797439508 or 3.7975772074093323e78 < z Initial program 12.6
rmApplied cosh-def_binary6412.6
Applied frac-times_binary6412.6
Applied associate-/l/_binary640.3
Simplified0.3
if -0.060131421797439508 < z < 3.7975772074093323e78Initial program 1.1
rmApplied cosh-def_binary641.1
Applied frac-times_binary641.1
Applied associate-/l/_binary6416.3
Simplified16.3
rmApplied clear-num_binary6416.4
Simplified1.0
Final simplification0.6
herbie shell --seed 2020273
(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))