\frac{\cosh x \cdot \frac{y}{x}}{z}\begin{array}{l}
\mathbf{if}\;y \le -2.1538486122792733 \cdot 10^{43}:\\
\;\;\;\;\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right) \cdot \left(\sqrt[3]{x} \cdot \left(\frac{1}{2} \cdot \frac{y}{z}\right)\right) + \frac{y}{x \cdot z}\\
\mathbf{elif}\;y \le 4.00822008275982678 \cdot 10^{-26}:\\
\;\;\;\;x \cdot \left(\frac{1}{2} \cdot \frac{y}{z}\right) + \frac{\frac{y}{x}}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{2} \cdot \left(e^{-1 \cdot x} + e^{x}\right)}{\frac{z \cdot x}{y}}\\
\end{array}double f(double x, double y, double z) {
double r611971 = x;
double r611972 = cosh(r611971);
double r611973 = y;
double r611974 = r611973 / r611971;
double r611975 = r611972 * r611974;
double r611976 = z;
double r611977 = r611975 / r611976;
return r611977;
}
double f(double x, double y, double z) {
double r611978 = y;
double r611979 = -2.1538486122792733e+43;
bool r611980 = r611978 <= r611979;
double r611981 = x;
double r611982 = cbrt(r611981);
double r611983 = r611982 * r611982;
double r611984 = 0.5;
double r611985 = z;
double r611986 = r611978 / r611985;
double r611987 = r611984 * r611986;
double r611988 = r611982 * r611987;
double r611989 = r611983 * r611988;
double r611990 = r611981 * r611985;
double r611991 = r611978 / r611990;
double r611992 = r611989 + r611991;
double r611993 = 4.008220082759827e-26;
bool r611994 = r611978 <= r611993;
double r611995 = r611981 * r611987;
double r611996 = r611978 / r611981;
double r611997 = r611996 / r611985;
double r611998 = r611995 + r611997;
double r611999 = -1.0;
double r612000 = r611999 * r611981;
double r612001 = exp(r612000);
double r612002 = exp(r611981);
double r612003 = r612001 + r612002;
double r612004 = r611984 * r612003;
double r612005 = r611985 * r611981;
double r612006 = r612005 / r611978;
double r612007 = r612004 / r612006;
double r612008 = r611994 ? r611998 : r612007;
double r612009 = r611980 ? r611992 : r612008;
return r612009;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.9 |
|---|---|
| Target | 0.5 |
| Herbie | 1.0 |
if y < -2.1538486122792733e+43Initial program 25.0
Taylor expanded around 0 1.6
rmApplied *-un-lft-identity1.6
Applied times-frac1.6
Applied associate-*r*1.6
Simplified1.6
rmApplied associate-*l*1.6
rmApplied add-cube-cbrt1.6
Applied associate-*l*1.6
if -2.1538486122792733e+43 < y < 4.008220082759827e-26Initial program 0.5
Taylor expanded around 0 11.1
rmApplied *-un-lft-identity11.1
Applied times-frac11.1
Applied associate-*r*11.1
Simplified11.1
rmApplied associate-*l*11.1
rmApplied associate-/r*1.1
if 4.008220082759827e-26 < y Initial program 19.7
Taylor expanded around inf 0.4
Simplified0.4
Final simplification1.0
herbie shell --seed 2020060
(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))