Quadratic roots, wide range

Average Error: 52.6 → 0.3

Time: 4.4s

Precision: binary64

\[\]

Math

\[\]

↓

\[\]

Error

Bits error versus a

Bits error versus b

Bits error versus c

Derivation

1. Initial program 52.6

   \[\]
2. Using strategy rm

3. Applied flip-+ 52.6

   \[\leadsto \]

4. Simplified 0.4

   \[\leadsto \]

5. Simplified 0.4

   \[\leadsto \]
6. Using strategy rm

7. Applied *-un-lft-identity 0.4

   \[\leadsto \]

8. Applied times-frac 0.4

   \[\leadsto \]

9. Applied times-frac 0.4

   \[\leadsto \]

10. Simplified 0.4

    \[\leadsto \]

11. Simplified 0.1

    \[\leadsto \]
12. Using strategy rm

13. Applied div-inv 0.3

    \[\leadsto \]

14. Final simplification 0.3

    \[\leadsto \]

Reproduce

herbie shell --seed 2020181 
(FPCore (a b c)
  :name "Quadratic roots, wide range"
  :precision binary64
  :pre (and (< 4.930380657631324e-32 a 2.028240960365167e+31) (< 4.930380657631324e-32 b 2.028240960365167e+31) (< 4.930380657631324e-32 c 2.028240960365167e+31))
  (/ (+ (neg b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))