Average Error: 33.5 → 10.3
Time: 18.5s
Precision: 64
Internal Precision: 128
\[\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\]
\[\begin{array}{l} \mathbf{if}\;b_2 \le -2.3604286795685533 \cdot 10^{-29}:\\ \;\;\;\;\frac{-1}{2} \cdot \frac{c}{b_2}\\ \mathbf{elif}\;b_2 \le 3.328175344273382 \cdot 10^{+63}:\\ \;\;\;\;\left(-\frac{b_2}{a}\right) - \frac{\sqrt{b_2 \cdot b_2 - c \cdot a}}{a}\\ \mathbf{else}:\\ \;\;\;\;\frac{c}{b_2} \cdot \frac{1}{2} - 2 \cdot \frac{b_2}{a}\\ \end{array}\]

Error

Bits error versus a

Bits error versus b_2

Bits error versus c

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 3 regimes
  2. if b_2 < -2.3604286795685533e-29

    1. Initial program 54.6

      \[\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\]
    2. Using strategy rm
    3. Applied div-sub55.4

      \[\leadsto \color{blue}{\frac{-b_2}{a} - \frac{\sqrt{b_2 \cdot b_2 - a \cdot c}}{a}}\]
    4. Taylor expanded around -inf 6.8

      \[\leadsto \color{blue}{\frac{-1}{2} \cdot \frac{c}{b_2}}\]

    if -2.3604286795685533e-29 < b_2 < 3.328175344273382e+63

    1. Initial program 15.2

      \[\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\]
    2. Using strategy rm
    3. Applied div-sub15.2

      \[\leadsto \color{blue}{\frac{-b_2}{a} - \frac{\sqrt{b_2 \cdot b_2 - a \cdot c}}{a}}\]

    if 3.328175344273382e+63 < b_2

    1. Initial program 37.3

      \[\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\]
    2. Taylor expanded around inf 5.0

      \[\leadsto \color{blue}{\frac{1}{2} \cdot \frac{c}{b_2} - 2 \cdot \frac{b_2}{a}}\]
  3. Recombined 3 regimes into one program.
  4. Final simplification10.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;b_2 \le -2.3604286795685533 \cdot 10^{-29}:\\ \;\;\;\;\frac{-1}{2} \cdot \frac{c}{b_2}\\ \mathbf{elif}\;b_2 \le 3.328175344273382 \cdot 10^{+63}:\\ \;\;\;\;\left(-\frac{b_2}{a}\right) - \frac{\sqrt{b_2 \cdot b_2 - c \cdot a}}{a}\\ \mathbf{else}:\\ \;\;\;\;\frac{c}{b_2} \cdot \frac{1}{2} - 2 \cdot \frac{b_2}{a}\\ \end{array}\]

Reproduce

herbie shell --seed 2019016 
(FPCore (a b_2 c)
  :name "quad2m (problem 3.2.1, negative)"
  (/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))

Details

Time bar (total: 17.6s)Debug log

sample224.0ms

Algorithm
intervals
Results
368×(pre true 80)
182×(body real 80)
91×(body nan 80)
35×(body real 1280)
32×(body real 2560)
13×(body real 640)
(body real 320)
(body real 160)
(body real 5120)

simplify41.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
41.0ms
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a)

prune12.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 35.0b

localize57.0ms

Local error

Found 4 expressions with local error:

23.7b
(sqrt (- (* b_2 b_2) (* a c)))
1.2b
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
1.2b
(- (* b_2 b_2) (* a c))
0.1b
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a)

rewrite53.0ms

Algorithm
rewrite-expression-head
Rules
15×add-sqr-sqrt
10×*-un-lft-identity
add-log-exp
associate-/l*
add-cube-cbrt
sqrt-prod
difference-of-squares
flip--
add-exp-log
add-cbrt-cube
flip3--
pow1
associate-/l/
distribute-lft-out--
diff-log
sqrt-div
sub-neg
div-inv
div-sub
pow1/2
frac-2neg
clear-num
rem-sqrt-square
neg-sub0
associate--l-
Counts
4 → 61
Calls
4 calls:
Slowest
29.0ms
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a)
11.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
6.0ms
(sqrt (- (* b_2 b_2) (* a c)))
4.0ms
(- (* b_2 b_2) (* a c))

series233.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
76.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
71.0ms
(sqrt (- (* b_2 b_2) (* a c)))
55.0ms
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a)
31.0ms
(- (* b_2 b_2) (* a c))

simplify1.5s

Counts
32 → 73
Calls
32 calls:
Slowest
267.0ms
(sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3)))
238.0ms
(sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c))))
233.0ms
(- (pow (- b_2) 3) (pow (sqrt (- (* b_2 b_2) (* a c))) 3))
175.0ms
(- (* (- b_2) (- b_2)) (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))))
157.0ms
(- b_2 (* 1/2 (/ (* a c) b_2)))

prune798.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 4.7b

localize6.0ms

Local error

Found 4 expressions with local error:

23.7b
(sqrt (- (* b_2 b_2) (* a c)))
1.2b
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
1.2b
(- (* b_2 b_2) (* a c))
0.2b
(/ 1 (/ a (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))

rewrite21.0ms

Algorithm
rewrite-expression-head
Rules
11×add-sqr-sqrt
add-log-exp
*-un-lft-identity
add-cube-cbrt
associate-/r*
add-exp-log
pow1
flip--
sqrt-prod
add-cbrt-cube
flip3--
associate-/r/
difference-of-squares
div-inv
diff-log
sqrt-div
sub-neg
inv-pow
distribute-lft-out--
pow-flip
pow1/2
frac-2neg
clear-num
rec-exp
rem-sqrt-square
neg-sub0
associate--l-
Counts
4 → 62
Calls
4 calls:
Slowest
11.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
4.0ms
(- (* b_2 b_2) (* a c))
3.0ms
(sqrt (- (* b_2 b_2) (* a c)))
2.0ms
(/ 1 (/ a (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))

series205.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
80.0ms
(/ 1 (/ a (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))
52.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
40.0ms
(sqrt (- (* b_2 b_2) (* a c)))
34.0ms
(- (* b_2 b_2) (* a c))

simplify1.6s

Counts
33 → 74
Calls
33 calls:
Slowest
688.0ms
(sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3)))
334.0ms
(+ (* (- b_2) (- b_2)) (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))
167.0ms
(sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c))))
113.0ms
(- b_2 (* 1/2 (/ (* a c) b_2)))
41.0ms
(- (* 1/2 (/ (* a c) b_2)) (* 2 b_2))

prune897.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 4.7b

localize13.0ms

Local error

Found 4 expressions with local error:

23.7b
(sqrt (- (* b_2 b_2) (* a c)))
1.2b
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
1.2b
(- (* b_2 b_2) (* a c))
0.2b
(* (/ 1 a) (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))

rewrite44.0ms

Algorithm
rewrite-expression-head
Rules
16×add-sqr-sqrt
11×*-un-lft-identity
add-log-exp
add-cube-cbrt
add-exp-log
associate-*r*
add-cbrt-cube
pow1
flip--
sqrt-prod
flip3--
difference-of-squares
associate-*l*
sub-neg
distribute-lft-out--
associate-*r/
diff-log
sqrt-div
distribute-lft-in
frac-times
distribute-rgt-in
div-inv
cbrt-unprod
*-commutative
prod-exp
associate-*l/
pow-prod-down
pow1/2
rem-sqrt-square
neg-sub0
associate--l-
Counts
4 → 72
Calls
4 calls:
Slowest
24.0ms
(* (/ 1 a) (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))
11.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
4.0ms
(sqrt (- (* b_2 b_2) (* a c)))
4.0ms
(- (* b_2 b_2) (* a c))

series197.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
73.0ms
(* (/ 1 a) (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))))
62.0ms
(sqrt (- (* b_2 b_2) (* a c)))
40.0ms
(- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))
21.0ms
(- (* b_2 b_2) (* a c))

simplify3.7s

Counts
47 → 84
Calls
47 calls:
Slowest
861.0ms
(* (* (* (/ 1 a) (/ 1 a)) (/ 1 a)) (* (* (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))) (- (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))
356.0ms
(+ (* (- b_2) (- b_2)) (+ (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* (- b_2) (sqrt (- (* b_2 b_2) (* a c))))))
352.0ms
(* 1 (- (* (- b_2) (- b_2)) (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c))))))
276.0ms
(* 1 (- (pow (- b_2) 3) (pow (sqrt (- (* b_2 b_2) (* a c))) 3)))
247.0ms
(sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3)))

prune822.0ms

Pruning

9 alts after pruning (8 fresh and 1 done)

Merged error: 4.7b

localize8.0ms

Local error

Found 4 expressions with local error:

23.7b
(sqrt (- (* b_2 b_2) (* a c)))
7.9b
(- (/ (- b_2) a) (/ (sqrt (- (* b_2 b_2) (* a c))) a))
1.2b
(- (* b_2 b_2) (* a c))
0.1b
(/ (sqrt (- (* b_2 b_2) (* a c))) a)

rewrite17.0ms

Algorithm
rewrite-expression-head
Rules
10×add-sqr-sqrt
*-un-lft-identity
add-log-exp
add-cube-cbrt
associate-/l*
sqrt-prod
flip--
add-exp-log
sqrt-div
add-cbrt-cube
flip3--
pow1
div-inv
associate-/l/
diff-log
sub-neg
difference-of-squares
sub-div
distribute-lft-out--
frac-sub
div-sub
pow1/2
distribute-rgt-out--
frac-2neg
clear-num
rem-sqrt-square
neg-sub0
associate--l-
Counts
4 → 62
Calls
4 calls:
Slowest
8.0ms
(- (/ (- b_2) a) (/ (sqrt (- (* b_2 b_2) (* a c))) a))
5.0ms
(/ (sqrt (- (* b_2 b_2) (* a c))) a)
2.0ms
(sqrt (- (* b_2 b_2) (* a c)))
2.0ms
(- (* b_2 b_2) (* a c))

series187.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
68.0ms
(/ (sqrt (- (* b_2 b_2) (* a c))) a)
49.0ms
(- (/ (- b_2) a) (/ (sqrt (- (* b_2 b_2) (* a c))) a))
39.0ms
(sqrt (- (* b_2 b_2) (* a c)))
32.0ms
(- (* b_2 b_2) (* a c))

simplify1.6s

Counts
36 → 74
Calls
36 calls:
Slowest
308.0ms
(sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3)))
247.0ms
(sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c))))
246.0ms
(sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3)))
213.0ms
(sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c))))
162.0ms
(- b_2 (* 1/2 (/ (* a c) b_2)))

prune795.0ms

Pruning

9 alts after pruning (7 fresh and 2 done)

Merged error: 4.7b

regimes375.0ms

Accuracy

83.2% (4.7b remaining)

Error of 10.3b against oracle of 5.7b and baseline of 33.6b

bsearch151.0ms

end0.0ms

sample4.0s

Algorithm
intervals
Results
11456×(pre true 80)
5804×(body real 80)
2736×(body nan 80)
1109×(body real 2560)
946×(body real 1280)
435×(body real 640)
215×(body real 320)
110×(body real 5120)
101×(body real 160)