Average Error: 44.2 → 11.4
Time: 20.1s
Precision: 64
\[1.1102230246251565 \cdot 10^{-16} \lt a \lt 9007199254740992.0 \land 1.1102230246251565 \cdot 10^{-16} \lt b \lt 9007199254740992.0 \land 1.1102230246251565 \cdot 10^{-16} \lt c \lt 9007199254740992.0\]
\[\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}\]
↓
\[\begin{array}{l} \mathbf{if}\;b \le 0.43696855946829866:\\ \;\;\;\;\frac{\frac{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} - b \cdot \left(b \cdot b\right)}{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) + \left(b \cdot b + b \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b}\right)}}{a \cdot 3}\\ \mathbf{else}:\\ \;\;\;\;\frac{-1}{2} \cdot \frac{c}{b}\\ \end{array}\]
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}
↓
\begin{array}{l}
\mathbf{if}\;b \le 0.43696855946829866:\\
\;\;\;\;\frac{\frac{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} - b \cdot \left(b \cdot b\right)}{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) + \left(b \cdot b + b \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b}\right)}}{a \cdot 3}\\

\mathbf{else}:\\
\;\;\;\;\frac{-1}{2} \cdot \frac{c}{b}\\

\end{array}
double f(double a, double b, double c, double __attribute__((unused)) d) {
        double r2910635 = b;
        double r2910636 = -r2910635;
        double r2910637 = r2910635 * r2910635;
        double r2910638 = 3.0;
        double r2910639 = a;
        double r2910640 = r2910638 * r2910639;
        double r2910641 = c;
        double r2910642 = r2910640 * r2910641;
        double r2910643 = r2910637 - r2910642;
        double r2910644 = sqrt(r2910643);
        double r2910645 = r2910636 + r2910644;
        double r2910646 = r2910645 / r2910640;
        return r2910646;
}
↓
double f(double a, double b, double c, double __attribute__((unused)) d) {
        double r2910647 = b;
        double r2910648 = 0.43696855946829866;
        bool r2910649 = r2910647 <= r2910648;
        double r2910650 = a;
        double r2910651 = c;
        double r2910652 = r2910650 * r2910651;
        double r2910653 = -3.0;
        double r2910654 = r2910652 * r2910653;
        double r2910655 = r2910647 * r2910647;
        double r2910656 = r2910654 + r2910655;
        double r2910657 = sqrt(r2910656);
        double r2910658 = r2910656 * r2910657;
        double r2910659 = r2910647 * r2910655;
        double r2910660 = r2910658 - r2910659;
        double r2910661 = r2910647 * r2910657;
        double r2910662 = r2910655 + r2910661;
        double r2910663 = r2910656 + r2910662;
        double r2910664 = r2910660 / r2910663;
        double r2910665 = 3.0;
        double r2910666 = r2910650 * r2910665;
        double r2910667 = r2910664 / r2910666;
        double r2910668 = -0.5;
        double r2910669 = r2910651 / r2910647;
        double r2910670 = r2910668 * r2910669;
        double r2910671 = r2910649 ? r2910667 : r2910670;
        return r2910671;
}

Error

Bits error versus a

Bits error versus b

Bits error versus c

Bits error versus d

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes

  2. if b < 0.43696855946829866

    1. Initial program 24.8

      \[\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}\]

    2. Simplified24.8

      \[\leadsto \color{blue}{\frac{\sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} - b}{3 \cdot a}}\]
    3. Using strategy rm

    4. Applied flip3--24.9

      \[\leadsto \frac{\color{blue}{\frac{{\left(\sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}\right)}^{3} - {b}^{3}}{\sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} \cdot \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} + \left(b \cdot b + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} \cdot b\right)}}}{3 \cdot a}\]

    5. Simplified24.3

      \[\leadsto \frac{\frac{\color{blue}{\sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} \cdot \left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) - \left(b \cdot b\right) \cdot b}}{\sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} \cdot \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} + \left(b \cdot b + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} \cdot b\right)}}{3 \cdot a}\]

    6. Simplified24.3

      \[\leadsto \frac{\frac{\sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} \cdot \left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) - \left(b \cdot b\right) \cdot b}{\color{blue}{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) + \left(b \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} + b \cdot b\right)}}}{3 \cdot a}\]

    if 0.43696855946829866 < b

    1. Initial program 47.6

      \[\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}\]

    2. Simplified47.6

      \[\leadsto \color{blue}{\frac{\sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c} - b}{3 \cdot a}}\]

    3. Taylor expanded around inf 9.2

      \[\leadsto \color{blue}{\frac{-1}{2} \cdot \frac{c}{b}}\]
  3. Recombined 2 regimes into one program.

  4. Final simplification11.4

    \[\leadsto \begin{array}{l} \mathbf{if}\;b \le 0.43696855946829866:\\ \;\;\;\;\frac{\frac{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b} - b \cdot \left(b \cdot b\right)}{\left(\left(a \cdot c\right) \cdot -3 + b \cdot b\right) + \left(b \cdot b + b \cdot \sqrt{\left(a \cdot c\right) \cdot -3 + b \cdot b}\right)}}{a \cdot 3}\\ \mathbf{else}:\\ \;\;\;\;\frac{-1}{2} \cdot \frac{c}{b}\\ \end{array}\]

Reproduce

herbie shell --seed 2019135 
(FPCore (a b c d)
  :name "Cubic critical, medium range"
  :pre (and (< 1.1102230246251565e-16 a 9007199254740992.0) (< 1.1102230246251565e-16 b 9007199254740992.0) (< 1.1102230246251565e-16 c 9007199254740992.0))
  (/ (+ (- b) (sqrt (- (* b b) (* (* 3 a) c)))) (* 3 a)))