Average Error: 33.7 → 10.0

Time: 24.6s

Precision: 64

Internal Precision: 128

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

↓

\[\begin{array}{l} \mathbf{if}\;b \le -9.2206393300312 \cdot 10^{+145}:\\ \;\;\;\;\frac{c}{b} - \frac{b}{a}\\ \mathbf{elif}\;b \le 1.5126665316535186 \cdot 10^{-110}:\\ \;\;\;\;\frac{\sqrt{{b}^{2} - \left(c \cdot a\right) \cdot 4} - b}{a \cdot 2}\\ \mathbf{else}:\\ \;\;\;\;-\frac{c}{b}\\ \end{array}\]

Error

The X axis uses an exponential scale — Bits error versus `a`

Try it out

In
Out

Enter valid numbers for all inputs

Target

Original	33.7
Target	21.0
Herbie	10.0

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

Derivation

Split input into 3 regimes
if b < -9.2206393300312e+145
1. Initial program 58.1
  \[\frac{\left(-b\right) + \sqrt{b \cdot b - 4 \cdot \left(a \cdot c\right)}}{2 \cdot a}\]
2. Initial simplification58.1
  \[\leadsto \frac{\sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c} - b}{2 \cdot a}\]
3. Taylor expanded around 0 58.1
  \[\leadsto \frac{\sqrt{\color{blue}{{b}^{2} - 4 \cdot \left(a \cdot c\right)}} - b}{2 \cdot a}\]
4. Taylor expanded around -inf 2.2
  \[\leadsto \color{blue}{\frac{c}{b} - \frac{b}{a}}\]
if -9.2206393300312e+145 < b < 1.5126665316535186e-110
1. Initial program 11.3
  \[\frac{\left(-b\right) + \sqrt{b \cdot b - 4 \cdot \left(a \cdot c\right)}}{2 \cdot a}\]
2. Initial simplification11.3
  \[\leadsto \frac{\sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c} - b}{2 \cdot a}\]
3. Taylor expanded around 0 11.3
  \[\leadsto \frac{\sqrt{\color{blue}{{b}^{2} - 4 \cdot \left(a \cdot c\right)}} - b}{2 \cdot a}\]
if 1.5126665316535186e-110 < b
1. Initial program 51.2
  \[\frac{\left(-b\right) + \sqrt{b \cdot b - 4 \cdot \left(a \cdot c\right)}}{2 \cdot a}\]
2. Initial simplification51.3
  \[\leadsto \frac{\sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c} - b}{2 \cdot a}\]
3. Taylor expanded around 0 51.2
  \[\leadsto \frac{\sqrt{\color{blue}{{b}^{2} - 4 \cdot \left(a \cdot c\right)}} - b}{2 \cdot a}\]
4. Taylor expanded around inf 10.9
  \[\leadsto \color{blue}{-1 \cdot \frac{c}{b}}\]
5. Simplified10.9
  \[\leadsto \color{blue}{\frac{-c}{b}}\]
Recombined 3 regimes into one program.
Final simplification10.0
\[\leadsto \begin{array}{l} \mathbf{if}\;b \le -9.2206393300312 \cdot 10^{+145}:\\ \;\;\;\;\frac{c}{b} - \frac{b}{a}\\ \mathbf{elif}\;b \le 1.5126665316535186 \cdot 10^{-110}:\\ \;\;\;\;\frac{\sqrt{{b}^{2} - \left(c \cdot a\right) \cdot 4} - b}{a \cdot 2}\\ \mathbf{else}:\\ \;\;\;\;-\frac{c}{b}\\ \end{array}\]

Runtime

	Baseline	Herbie	Oracle	Span	%
Regimes	33.7	10.0	5.6	28.2	84.1%

herbie shell --seed 2018355 
(FPCore (a b c)
  :name "quadp (p42, positive)"

  :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)))

Error

Try it out

Target

Derivation

`if b < -9.2206393300312e+145`

`if -9.2206393300312e+145 < b < 1.5126665316535186e-110`

`if 1.5126665316535186e-110 < b`

Runtime