The quadratic formula (r1)

Average Error: 33.0 → 8.7

Time: 20.1s

Precision: 64

Math

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

↓

\[\begin{array}{l} \mathbf{if}\;b \le -1.6519381339788066 \cdot 10^{+37}:\\ \;\;\;\;\frac{\left(\frac{c}{b} - \frac{b}{a}\right) \cdot 2}{2}\\ \mathbf{elif}\;b \le -1.3761661522305357 \cdot 10^{-153}:\\ \;\;\;\;\frac{\frac{\sqrt{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b}}{\frac{a}{\sqrt{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b}}}}{2}\\ \mathbf{elif}\;b \le 6.555431533807236 \cdot 10^{+28}:\\ \;\;\;\;\frac{\frac{a \cdot \frac{-4 \cdot c}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} + b}}{a}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{-2 \cdot \frac{c}{b}}{2}\\ \end{array}\]

C

double f(double a, double b, double c) {
        double r2842197 = b;
        double r2842198 = -r2842197;
        double r2842199 = r2842197 * r2842197;
        double r2842200 = 4.0;
        double r2842201 = a;
        double r2842202 = r2842200 * r2842201;
        double r2842203 = c;
        double r2842204 = r2842202 * r2842203;
        double r2842205 = r2842199 - r2842204;
        double r2842206 = sqrt(r2842205);
        double r2842207 = r2842198 + r2842206;
        double r2842208 = 2.0;
        double r2842209 = r2842208 * r2842201;
        double r2842210 = r2842207 / r2842209;
        return r2842210;
}

↓

double f(double a, double b, double c) {
        double r2842211 = b;
        double r2842212 = -1.6519381339788066e+37;
        bool r2842213 = r2842211 <= r2842212;
        double r2842214 = c;
        double r2842215 = r2842214 / r2842211;
        double r2842216 = a;
        double r2842217 = r2842211 / r2842216;
        double r2842218 = r2842215 - r2842217;
        double r2842219 = 2.0;
        double r2842220 = r2842218 * r2842219;
        double r2842221 = r2842220 / r2842219;
        double r2842222 = -1.3761661522305357e-153;
        bool r2842223 = r2842211 <= r2842222;
        double r2842224 = -4.0;
        double r2842225 = r2842224 * r2842214;
        double r2842226 = r2842211 * r2842211;
        double r2842227 = fma(r2842225, r2842216, r2842226);
        double r2842228 = sqrt(r2842227);
        double r2842229 = r2842228 - r2842211;
        double r2842230 = sqrt(r2842229);
        double r2842231 = r2842216 / r2842230;
        double r2842232 = r2842230 / r2842231;
        double r2842233 = r2842232 / r2842219;
        double r2842234 = 6.555431533807236e+28;
        bool r2842235 = r2842211 <= r2842234;
        double r2842236 = r2842228 + r2842211;
        double r2842237 = r2842225 / r2842236;
        double r2842238 = r2842216 * r2842237;
        double r2842239 = r2842238 / r2842216;
        double r2842240 = r2842239 / r2842219;
        double r2842241 = -2.0;
        double r2842242 = r2842241 * r2842215;
        double r2842243 = r2842242 / r2842219;
        double r2842244 = r2842235 ? r2842240 : r2842243;
        double r2842245 = r2842223 ? r2842233 : r2842244;
        double r2842246 = r2842213 ? r2842221 : r2842245;
        return r2842246;
}

Error

Bits error versus a

Bits error versus b

Bits error versus c

Target

Original	33.0
Target	20.1
Herbie	8.7

\[\begin{array}{l} \mathbf{if}\;b \lt 0:\\ \;\;\;\;\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}\\ \mathbf{else}:\\ \;\;\;\;\frac{c}{a \cdot \frac{\left(-b\right) - \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}}\\ \end{array}\]

Derivation

1. Split input into 4 regimes

2. if b < -1.6519381339788066e+37

   1. Initial program 33.6

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

   2. Simplified 33.6

      \[\leadsto \color{blue}{\frac{\frac{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}{a}}{2}}\]

   3. Taylor expanded around -inf 6.5

      \[\leadsto \frac{\color{blue}{2 \cdot \frac{c}{b} - 2 \cdot \frac{b}{a}}}{2}\]

   4. Simplified 6.5

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

   if -1.6519381339788066e+37 < b < -1.3761661522305357e-153

   1. Initial program 5.2

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

   2. Simplified 5.2

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

   4. Applied clear-num 5.3

      \[\leadsto \frac{\color{blue}{\frac{1}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]
   5. Using strategy rm

   6. Applied *-un-lft-identity 5.3

      \[\leadsto \frac{\frac{1}{\frac{a}{\color{blue}{1 \cdot \left(\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b\right)}}}}{2}\]

   7. Applied *-un-lft-identity 5.3

      \[\leadsto \frac{\frac{1}{\frac{\color{blue}{1 \cdot a}}{1 \cdot \left(\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b\right)}}}{2}\]

   8. Applied times-frac 5.3

      \[\leadsto \frac{\frac{1}{\color{blue}{\frac{1}{1} \cdot \frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]

   9. Applied add-cube-cbrt 5.3

      \[\leadsto \frac{\frac{\color{blue}{\left(\sqrt[3]{1} \cdot \sqrt[3]{1}\right) \cdot \sqrt[3]{1}}}{\frac{1}{1} \cdot \frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}{2}\]

   10. Applied times-frac 5.3

       \[\leadsto \frac{\color{blue}{\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\frac{1}{1}} \cdot \frac{\sqrt[3]{1}}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]

   11. Simplified 5.3

       \[\leadsto \frac{\color{blue}{1} \cdot \frac{\sqrt[3]{1}}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}{2}\]

   12. Simplified 5.2

       \[\leadsto \frac{1 \cdot \color{blue}{\frac{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b}{a}}}{2}\]
   13. Using strategy rm

   14. Applied add-sqr-sqrt 5.6

       \[\leadsto \frac{1 \cdot \frac{\color{blue}{\sqrt{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b} \cdot \sqrt{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b}}}{a}}{2}\]

   15. Applied associate-/l* 5.6

       \[\leadsto \frac{1 \cdot \color{blue}{\frac{\sqrt{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b}}{\frac{a}{\sqrt{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b}}}}}{2}\]

   if -1.3761661522305357e-153 < b < 6.555431533807236e+28

   1. Initial program 24.5

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

   2. Simplified 24.5

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

   4. Applied clear-num 24.5

      \[\leadsto \frac{\color{blue}{\frac{1}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]
   5. Using strategy rm

   6. Applied *-un-lft-identity 24.5

      \[\leadsto \frac{\frac{1}{\frac{a}{\color{blue}{1 \cdot \left(\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b\right)}}}}{2}\]

   7. Applied *-un-lft-identity 24.5

      \[\leadsto \frac{\frac{1}{\frac{\color{blue}{1 \cdot a}}{1 \cdot \left(\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b\right)}}}{2}\]

   8. Applied times-frac 24.5

      \[\leadsto \frac{\frac{1}{\color{blue}{\frac{1}{1} \cdot \frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]

   9. Applied add-cube-cbrt 24.5

      \[\leadsto \frac{\frac{\color{blue}{\left(\sqrt[3]{1} \cdot \sqrt[3]{1}\right) \cdot \sqrt[3]{1}}}{\frac{1}{1} \cdot \frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}{2}\]

   10. Applied times-frac 24.5

       \[\leadsto \frac{\color{blue}{\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{\frac{1}{1}} \cdot \frac{\sqrt[3]{1}}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}}{2}\]

   11. Simplified 24.5

       \[\leadsto \frac{\color{blue}{1} \cdot \frac{\sqrt[3]{1}}{\frac{a}{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}}}{2}\]

   12. Simplified 24.5

       \[\leadsto \frac{1 \cdot \color{blue}{\frac{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b}{a}}}{2}\]
   13. Using strategy rm

   14. Applied flip-- 25.0

       \[\leadsto \frac{1 \cdot \frac{\color{blue}{\frac{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} \cdot \sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} - b \cdot b}{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} + b}}}{a}}{2}\]

   15. Simplified 17.1

       \[\leadsto \frac{1 \cdot \frac{\frac{\color{blue}{a \cdot \left(c \cdot -4\right)}}{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} + b}}{a}}{2}\]
   16. Using strategy rm

   17. Applied *-un-lft-identity 17.1

       \[\leadsto \frac{1 \cdot \frac{\frac{a \cdot \left(c \cdot -4\right)}{\color{blue}{1 \cdot \left(\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} + b\right)}}}{a}}{2}\]

   18. Applied times-frac 14.9

       \[\leadsto \frac{1 \cdot \frac{\color{blue}{\frac{a}{1} \cdot \frac{c \cdot -4}{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} + b}}}{a}}{2}\]

   19. Simplified 14.9

       \[\leadsto \frac{1 \cdot \frac{\color{blue}{a} \cdot \frac{c \cdot -4}{\sqrt{\mathsf{fma}\left(c \cdot -4, a, b \cdot b\right)} + b}}{a}}{2}\]

   if 6.555431533807236e+28 < b

   1. Initial program 56.3

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

   2. Simplified 56.3

      \[\leadsto \color{blue}{\frac{\frac{\sqrt{\mathsf{fma}\left(b, b, \left(a \cdot -4\right) \cdot c\right)} - b}{a}}{2}}\]

   3. Taylor expanded around inf 4.5

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

4. Final simplification 8.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \le -1.6519381339788066 \cdot 10^{+37}:\\ \;\;\;\;\frac{\left(\frac{c}{b} - \frac{b}{a}\right) \cdot 2}{2}\\ \mathbf{elif}\;b \le -1.3761661522305357 \cdot 10^{-153}:\\ \;\;\;\;\frac{\frac{\sqrt{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b}}{\frac{a}{\sqrt{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} - b}}}}{2}\\ \mathbf{elif}\;b \le 6.555431533807236 \cdot 10^{+28}:\\ \;\;\;\;\frac{\frac{a \cdot \frac{-4 \cdot c}{\sqrt{\mathsf{fma}\left(-4 \cdot c, a, b \cdot b\right)} + b}}{a}}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{-2 \cdot \frac{c}{b}}{2}\\ \end{array}\]

Reproduce

herbie shell --seed 2019164 +o rules:numerics
(FPCore (a b c)
  :name "The quadratic formula (r1)"

  :herbie-target
  (if (< b 0) (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) (/ c (* a (/ (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))))

  (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))