Average Error: 7.5 → 0.4
Time: 23.3s
Precision: 64
\[\frac{\cosh x \cdot \frac{y}{x}}{z}\]
↓
\[\begin{array}{l} \mathbf{if}\;z \le -7.499151431129345690915331400810285626293 \cdot 10^{-22} \lor \neg \left(z \le 3.123371243102039260157787505142381731931 \cdot 10^{-43}\right):\\ \;\;\;\;\frac{y}{\frac{z \cdot x}{\left(e^{x} + e^{-x}\right) \cdot \frac{1}{2}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{y}{\frac{z}{e^{x} + e^{-x}}}}{\frac{x}{\frac{1}{2}}}\\ \end{array}\]
\frac{\cosh x \cdot \frac{y}{x}}{z}
↓
\begin{array}{l}
\mathbf{if}\;z \le -7.499151431129345690915331400810285626293 \cdot 10^{-22} \lor \neg \left(z \le 3.123371243102039260157787505142381731931 \cdot 10^{-43}\right):\\
\;\;\;\;\frac{y}{\frac{z \cdot x}{\left(e^{x} + e^{-x}\right) \cdot \frac{1}{2}}}\\

\mathbf{else}:\\
\;\;\;\;\frac{\frac{y}{\frac{z}{e^{x} + e^{-x}}}}{\frac{x}{\frac{1}{2}}}\\

\end{array}
double f(double x, double y, double z) {
        double r396306 = x;
        double r396307 = cosh(r396306);
        double r396308 = y;
        double r396309 = r396308 / r396306;
        double r396310 = r396307 * r396309;
        double r396311 = z;
        double r396312 = r396310 / r396311;
        return r396312;
}
↓
double f(double x, double y, double z) {
        double r396313 = z;
        double r396314 = -7.499151431129346e-22;
        bool r396315 = r396313 <= r396314;
        double r396316 = 3.1233712431020393e-43;
        bool r396317 = r396313 <= r396316;
        double r396318 = !r396317;
        bool r396319 = r396315 || r396318;
        double r396320 = y;
        double r396321 = x;
        double r396322 = r396313 * r396321;
        double r396323 = exp(r396321);
        double r396324 = -r396321;
        double r396325 = exp(r396324);
        double r396326 = r396323 + r396325;
        double r396327 = 0.5;
        double r396328 = r396326 * r396327;
        double r396329 = r396322 / r396328;
        double r396330 = r396320 / r396329;
        double r396331 = r396313 / r396326;
        double r396332 = r396320 / r396331;
        double r396333 = r396321 / r396327;
        double r396334 = r396332 / r396333;
        double r396335 = r396319 ? r396330 : r396334;
        return r396335;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original7.5
Target0.4
Herbie0.4
\[\begin{array}{l} \mathbf{if}\;y \lt -4.618902267687041990497740832940559043667 \cdot 10^{-52}:\\ \;\;\;\;\frac{\frac{y}{z}}{x} \cdot \cosh x\\ \mathbf{elif}\;y \lt 1.038530535935153018369520384190862667426 \cdot 10^{-39}:\\ \;\;\;\;\frac{\frac{\cosh x \cdot y}{x}}{z}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{y}{z}}{x} \cdot \cosh x\\ \end{array}\]

Derivation

  1. Split input into 2 regimes

  2. if z < -7.499151431129346e-22 or 3.1233712431020393e-43 < z

    1. Initial program 10.7

      \[\frac{\cosh x \cdot \frac{y}{x}}{z}\]

    2. Taylor expanded around inf 0.4

      \[\leadsto \color{blue}{\frac{y \cdot \left(\frac{1}{2} \cdot e^{x} + \frac{1}{2} \cdot e^{-x}\right)}{x \cdot z}}\]

    3. Simplified0.4

      \[\leadsto \color{blue}{\frac{y}{\frac{z \cdot x}{\left(e^{x} + e^{-x}\right) \cdot \frac{1}{2}}}}\]

    if -7.499151431129346e-22 < z < 3.1233712431020393e-43

    1. Initial program 0.3

      \[\frac{\cosh x \cdot \frac{y}{x}}{z}\]

    2. Taylor expanded around inf 22.1

      \[\leadsto \color{blue}{\frac{y \cdot \left(\frac{1}{2} \cdot e^{x} + \frac{1}{2} \cdot e^{-x}\right)}{x \cdot z}}\]

    3. Simplified22.1

      \[\leadsto \color{blue}{\frac{y}{\frac{z \cdot x}{\left(e^{x} + e^{-x}\right) \cdot \frac{1}{2}}}}\]
    4. Using strategy rm

    5. Applied times-frac22.1

      \[\leadsto \frac{y}{\color{blue}{\frac{z}{e^{x} + e^{-x}} \cdot \frac{x}{\frac{1}{2}}}}\]

    6. Applied associate-/r*0.3

      \[\leadsto \color{blue}{\frac{\frac{y}{\frac{z}{e^{x} + e^{-x}}}}{\frac{x}{\frac{1}{2}}}}\]
  3. Recombined 2 regimes into one program.

  4. Final simplification0.4

    \[\leadsto \begin{array}{l} \mathbf{if}\;z \le -7.499151431129345690915331400810285626293 \cdot 10^{-22} \lor \neg \left(z \le 3.123371243102039260157787505142381731931 \cdot 10^{-43}\right):\\ \;\;\;\;\frac{y}{\frac{z \cdot x}{\left(e^{x} + e^{-x}\right) \cdot \frac{1}{2}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{y}{\frac{z}{e^{x} + e^{-x}}}}{\frac{x}{\frac{1}{2}}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019323 
(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))