Average Error: 29.7 → 7.5
Time: 6.2s
Precision: 64
Internal Precision: 128
\[\sqrt{\left|\frac{a \cdot a - b \cdot b}{a \cdot a}\right|}\]
↓
\[\sqrt{\left|1 - \frac{b}{a} \cdot \frac{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.7

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

  2. Simplified7.5

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

  3. Final simplification7.5

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

Reproduce

herbie shell --seed 2019010 
(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: 5.5s)Debug log

sample36.0ms

Algorithm
intervals

simplify169.0ms

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

prune9.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 9.2b

localize13.0ms

Local error

Found 3 expressions with local error:

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

rewrite6.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
3.0ms
(* (/ b a) (/ b a))
1.0ms
(- 1 (* (/ b a) (/ b a)))
1.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

series52.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
24.0ms
(- 1 (* (/ b a) (/ b a)))
14.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
13.0ms
(* (/ b a) (/ b a))

simplify1.4s

Counts
28 → 57
Calls
28 calls:
Slowest
925.0ms
(* (* (* (/ b a) (/ b a)) (/ b a)) (* (* (/ b a) (/ b a)) (/ b a)))
67.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
63.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
55.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
53.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

prune624.0ms

Pruning

2 alts after pruning (1 fresh and 1 done)

Merged error: 9.2b

localize20.0ms

Local error

Found 4 expressions with local error:

13.9b
(cbrt (* (* (sqrt (fabs (- 1 (* (/ b a) (/ b a))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))))
9.2b
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
9.2b
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
9.2b
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

rewrite20.0ms

Algorithm
rewrite-expression-head
Rules
10×add-sqr-sqrt
9×sqrt-prod
7×add-cube-cbrt
7×*-un-lft-identity
4×add-log-exp
4×add-exp-log
4×add-cbrt-cube
4×pow1
3×pow1/2
3×rem-sqrt-square
1×rem-cbrt-cube
1×pow3
1×pow1/3
1×cbrt-prod
Counts
4 → 46
Calls
4 calls:
Slowest
16.0ms
(cbrt (* (* (sqrt (fabs (- 1 (* (/ b a) (/ b a))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))))
1.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
1.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
1.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

series82.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
24.0ms
(cbrt (* (* (sqrt (fabs (- 1 (* (/ b a) (/ b a))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))) (sqrt (fabs (- 1 (* (/ b a) (/ b a)))))))
23.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
22.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))
13.0ms
(sqrt (fabs (- 1 (* (/ b a) (/ b a)))))

simplify1.2s

Counts
26 → 58
Calls
26 calls:
Slowest
68.0ms
(sqrt (* (cbrt (fabs (- 1 (* (/ b a) (/ b a))))) (cbrt (fabs (- 1 (* (/ b a) (/ b a)))))))
67.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
67.0ms
(sqrt (fabs (- 1 (/ (pow b 2) (pow a 2)))))
65.0ms
(sqrt (sqrt (fabs (- 1 (* (/ b a) (/ b a))))))
64.0ms
(sqrt (sqrt (fabs (- 1 (* (/ b a) (/ b a))))))

prune735.0ms

Pruning

2 alts after pruning (0 fresh and 2 done)

Merged error: 9.2b

regimes54.0ms

Accuracy

0% (0.0b remaining)

Error of 7.5b against oracle of 7.5b and baseline of 7.5b

bsearch2.0ms

end0.0ms

sample1.1s

Algorithm
intervals