Average Error: 29.3 → 7.4
Time: 16.5s
Precision: 64
Internal Precision: 128
\[\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\]
↓
\[\sqrt{\left|1 - \frac{\frac{b}{a} \cdot b}{a}\right|}\]

Error

Bits error versus a

Bits error versus b

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 29.3

    \[\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\]

  2. Simplified7.4

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

  4. Applied associate-*l/7.4

    \[\leadsto \sqrt{\left|1 - \color{blue}{\frac{b \cdot \frac{b}{a}}{a}}\right|}\]

  5. Final simplification7.4

    \[\leadsto \sqrt{\left|1 - \frac{\frac{b}{a} \cdot b}{a}\right|}\]

Reproduce

herbie shell --seed 2019022 
(FPCore (a b)
  :name "Eccentricity of an ellipse"
  :pre (<= 0 b a 1)
  (sqrt (fabs (/ (- (* a a) (* b b)) (* a a)))))

Details

Time bar (total: 16.1s)Debug log

sample50.0ms

Algorithm
intervals
Results
29.0ms256×body80valid
10.0ms256×pre80true

simplify190.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
190.0ms
(sqrt (fabs (/ (- (* a a) (* b b)) (* a a))))

prune22.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 9.9b

localize26.0ms

Local error

Found 3 expressions with local error:

10.0b
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
0.4b
(* (/ b a) (/ b a))
0.0b
(- 1 (* (/ b a) (/ b a)))

rewrite10.0ms

Algorithm
rewrite-expression-head
Rules
8×pow1
7×add-sqr-sqrt
6×add-cube-cbrt
6×*-un-lft-identity
5×add-exp-log
5×add-cbrt-cube
4×associate-*l*
4×associate-*r*
3×add-log-exp
3×sqrt-prod
2×div-inv
1×flip--
1×cbrt-unprod
1×pow-prod-up
1×*-commutative
1×associate-*r/
1×prod-exp
1×associate-*l/
1×pow-prod-down
1×pow1/2
1×flip3--
1×pow-plus
1×sub-neg
1×rem-sqrt-square
1×frac-times
1×pow2
Counts
3 → 48
Calls
3 calls:
Slowest
6.0ms
(* (/ b a) (/ b a))
2.0ms
(- 1 (* (/ b a) (/ b a)))
1.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

series72.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
27.0ms
(- 1 (* (/ b a) (/ b a)))
23.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
23.0ms
(* (/ b a) (/ b a))

simplify1.3s

Counts
28 → 57
Calls
28 calls:
Slowest
835.0ms
(* (* (* (/ b a) (/ b a)) (/ b a)) (* (* (/ b a) (/ b a)) (/ b a)))
59.0ms
(sqrt (sqrt (fabs (- 1 (* (/ b a) (/ b a))))))
51.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
50.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
50.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))

prune537.0ms

Pruning

3 alts after pruning (2 fresh and 1 done)

Merged error: 9.9b

localize12.0ms

Local error

Found 4 expressions with local error:

10.0b
(sqrt (fabs (- 1 (/ (* b (/ b a)) a))))
6.9b
(/ (* b (/ b a)) a)
0.2b
(* b (/ b a))
0.0b
(- 1 (/ (* b (/ b a)) a))

rewrite10.0ms

Algorithm
rewrite-expression-head
Rules
7×add-sqr-sqrt
6×add-cube-cbrt
6×*-un-lft-identity
4×add-log-exp
4×add-exp-log
4×associate-*r*
4×add-cbrt-cube
4×pow1
3×sqrt-prod
2×div-inv
2×associate-*r/
1×associate-/l/
1×associate-/l*
1×flip--
1×*-commutative
1×pow1/2
1×flip3--
1×frac-2neg
1×sub-neg
1×clear-num
1×rem-sqrt-square
Counts
4 → 48
Calls
4 calls:
Slowest
5.0ms
(/ (* b (/ b a)) a)
2.0ms
(- 1 (/ (* b (/ b a)) a))
2.0ms
(* b (/ b a))
1.0ms
(sqrt (fabs (- 1 (/ (* b (/ b a)) a))))

series66.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
22.0ms
(sqrt (fabs (- 1 (/ (* b (/ b a)) a))))
18.0ms
(- 1 (/ (* b (/ b a)) a))
14.0ms
(/ (* b (/ b a)) a)
12.0ms
(* b (/ b a))

simplify426.0ms

Counts
21 → 60
Calls
21 calls:
Slowest
64.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
50.0ms
(sqrt (* (cbrt (fabs (- 1 (/ (* b (/ b a)) a)))) (cbrt (fabs (- 1 (/ (* b (/ b a)) a))))))
50.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
48.0ms
(sqrt (fabs (- 1 (/ (* b (/ b a)) a))))
39.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))

prune581.0ms

Pruning

3 alts after pruning (1 fresh and 2 done)

Merged error: 9.9b

localize26.0ms

Local error

Found 4 expressions with local error:

18.9b
(/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))
10.0b
(sqrt (fabs (/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))))
0.4b
(* (/ b a) (/ b a))
0.4b
(* (/ b a) (/ b a))

rewrite16.0ms

Algorithm
rewrite-expression-head
Rules
19×add-sqr-sqrt
18×add-cube-cbrt
18×*-un-lft-identity
14×pow1
12×times-frac
10×add-exp-log
10×add-cbrt-cube
8×associate-*l*
8×associate-*r*
5×div-inv
4×difference-of-squares
4×add-log-exp
4×associate-/l*
3×sqrt-prod
3×associate-/r*
2×associate-/l/
2×cbrt-unprod
2×pow-prod-up
2×*-commutative
2×associate-/r/
2×associate-*r/
2×prod-exp
2×associate-*l/
2×pow-prod-down
2×pow-plus
2×frac-times
2×pow2
1×flip--
1×flip-+
1×div-exp
1×div-sub
1×pow1/2
1×flip3--
1×flip3-+
1×frac-2neg
1×clear-num
1×rem-sqrt-square
1×cbrt-undiv
Counts
4 → 98
Calls
4 calls:
Slowest
8.0ms
(/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))
3.0ms
(* (/ b a) (/ b a))
3.0ms
(* (/ b a) (/ b a))
1.0ms
(sqrt (fabs (/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))))

series129.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
52.0ms
(/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))
34.0ms
(sqrt (fabs (/ (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))) (+ 1 (* (/ b a) (/ b a))))))
23.0ms
(* (/ b a) (/ b a))
19.0ms
(* (/ b a) (/ b a))

simplify10.4s

Counts
83 → 110
Calls
83 calls:
Slowest
781.0ms
(* (* (* (/ b a) (/ b a)) (/ b a)) (* (* (/ b a) (/ b a)) (/ b a)))
778.0ms
(* (* (* (/ b a) (/ b a)) (/ b a)) (* (* (/ b a) (/ b a)) (/ b a)))
751.0ms
(/ (cbrt (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a))))) (+ 1 (* (/ b a) (/ b a))))
496.0ms
(- (pow (* 1 1) 3) (pow (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a))) 3))
449.0ms
(/ (* (cbrt (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a))))) (cbrt (- (* 1 1) (* (* (/ b a) (/ b a)) (* (/ b a) (/ b a)))))) (sqrt (+ 1 (* (/ b a) (/ b a)))))

prune1.1s

Pruning

3 alts after pruning (0 fresh and 3 done)

Merged error: 9.9b

regimes52.0ms

Accuracy

0% (0.0b remaining)

Error of 7.4b against oracle of 7.4b and baseline of 7.4b

bsearch1.0ms

end0.0ms

sample1.2s

Algorithm
intervals
Results
819.0ms8000×body80valid
286.0ms8000×pre80true