Average Error: 7.8 → 1.0
Time: 19.8s
Precision: 64
\[\frac{\cosh x \cdot \frac{y}{x}}{z}\]
↓
\[\begin{array}{l} \mathbf{if}\;z \le -712742474614758 \lor \neg \left(z \le 3.970994457643005255735446389224709150426 \cdot 10^{112}\right):\\ \;\;\;\;\frac{e^{x} + e^{-x}}{z \cdot \left(2 \cdot x\right)} \cdot y\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{\left(e^{x} + e^{-x}\right) \cdot y}{z}}{2 \cdot x}\\ \end{array}\]
\frac{\cosh x \cdot \frac{y}{x}}{z}
↓
\begin{array}{l}
\mathbf{if}\;z \le -712742474614758 \lor \neg \left(z \le 3.970994457643005255735446389224709150426 \cdot 10^{112}\right):\\
\;\;\;\;\frac{e^{x} + e^{-x}}{z \cdot \left(2 \cdot x\right)} \cdot y\\

\mathbf{else}:\\
\;\;\;\;\frac{\frac{\left(e^{x} + e^{-x}\right) \cdot y}{z}}{2 \cdot x}\\

\end{array}
double f(double x, double y, double z) {
        double r346534 = x;
        double r346535 = cosh(r346534);
        double r346536 = y;
        double r346537 = r346536 / r346534;
        double r346538 = r346535 * r346537;
        double r346539 = z;
        double r346540 = r346538 / r346539;
        return r346540;
}
↓
double f(double x, double y, double z) {
        double r346541 = z;
        double r346542 = -712742474614758.0;
        bool r346543 = r346541 <= r346542;
        double r346544 = 3.970994457643005e+112;
        bool r346545 = r346541 <= r346544;
        double r346546 = !r346545;
        bool r346547 = r346543 || r346546;
        double r346548 = x;
        double r346549 = exp(r346548);
        double r346550 = -r346548;
        double r346551 = exp(r346550);
        double r346552 = r346549 + r346551;
        double r346553 = 2.0;
        double r346554 = r346553 * r346548;
        double r346555 = r346541 * r346554;
        double r346556 = r346552 / r346555;
        double r346557 = y;
        double r346558 = r346556 * r346557;
        double r346559 = r346552 * r346557;
        double r346560 = r346559 / r346541;
        double r346561 = r346560 / r346554;
        double r346562 = r346547 ? r346558 : r346561;
        return r346562;
}

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.8
Target0.5
Herbie1.0
\[\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 < -712742474614758.0 or 3.970994457643005e+112 < z

    1. Initial program 13.3

      \[\frac{\cosh x \cdot \frac{y}{x}}{z}\]
    2. Using strategy rm

    3. Applied cosh-def13.3

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

    4. Applied frac-times13.3

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

    5. Applied associate-/l/0.4

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

    7. Applied associate-/l*0.8

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

    9. Applied associate-/r/0.4

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

    if -712742474614758.0 < z < 3.970994457643005e+112

    1. Initial program 1.8

      \[\frac{\cosh x \cdot \frac{y}{x}}{z}\]
    2. Using strategy rm

    3. Applied cosh-def1.8

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

    4. Applied frac-times1.8

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

    5. Applied associate-/l/14.0

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

    7. Applied associate-/r*1.6

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

  4. Final simplification1.0

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

Reproduce

herbie shell --seed 2019303 +o rules:numerics
(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.03853053593515302e-39) (/ (/ (* (cosh x) y) x) z) (* (/ (/ y z) x) (cosh x))))

  (/ (* (cosh x) (/ y x)) z))