Average Error: 18.7 → 13.4
Time: 15.6s
Precision: 64
Internal Precision: 128
\[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
\[\begin{array}{l} \mathbf{if}\;V \cdot \ell \le -1.125220159166195 \cdot 10^{-80}:\\ \;\;\;\;\sqrt{\sqrt{\frac{A}{V \cdot \ell}}} \cdot \left(c0 \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}\right)\\ \mathbf{elif}\;V \cdot \ell \le 2.2677613144113 \cdot 10^{-321}:\\ \;\;\;\;c0 \cdot \sqrt{\frac{\frac{A}{\ell}}{V}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt{A}}{\sqrt{V \cdot \ell}} \cdot c0\\ \end{array}\]

Error

Bits error versus c0

Bits error versus A

Bits error versus V

Bits error versus l

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 3 regimes
  2. if (* V l) < -1.125220159166195e-80

    1. Initial program 15.1

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm
    3. Applied div-inv15.1

      \[\leadsto c0 \cdot \sqrt{\color{blue}{A \cdot \frac{1}{V \cdot \ell}}}\]
    4. Using strategy rm
    5. Applied add-sqr-sqrt15.3

      \[\leadsto c0 \cdot \color{blue}{\left(\sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}} \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}\right)}\]
    6. Applied associate-*r*15.2

      \[\leadsto \color{blue}{\left(c0 \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}\right) \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}}\]
    7. Simplified15.2

      \[\leadsto \left(c0 \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}\right) \cdot \color{blue}{\sqrt{\sqrt{\frac{A}{V \cdot \ell}}}}\]

    if -1.125220159166195e-80 < (* V l) < 2.2677613144113e-321

    1. Initial program 33.8

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm
    3. Applied div-inv34.1

      \[\leadsto c0 \cdot \sqrt{\color{blue}{A \cdot \frac{1}{V \cdot \ell}}}\]
    4. Using strategy rm
    5. Applied *-un-lft-identity34.1

      \[\leadsto c0 \cdot \color{blue}{\left(1 \cdot \sqrt{A \cdot \frac{1}{V \cdot \ell}}\right)}\]
    6. Applied associate-*r*34.1

      \[\leadsto \color{blue}{\left(c0 \cdot 1\right) \cdot \sqrt{A \cdot \frac{1}{V \cdot \ell}}}\]
    7. Simplified25.7

      \[\leadsto \left(c0 \cdot 1\right) \cdot \color{blue}{\sqrt{\frac{\frac{A}{\ell}}{V}}}\]

    if 2.2677613144113e-321 < (* V l)

    1. Initial program 14.3

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm
    3. Applied div-inv14.6

      \[\leadsto c0 \cdot \sqrt{\color{blue}{A \cdot \frac{1}{V \cdot \ell}}}\]
    4. Using strategy rm
    5. Applied un-div-inv14.3

      \[\leadsto c0 \cdot \sqrt{\color{blue}{\frac{A}{V \cdot \ell}}}\]
    6. Applied sqrt-div6.4

      \[\leadsto c0 \cdot \color{blue}{\frac{\sqrt{A}}{\sqrt{V \cdot \ell}}}\]
  3. Recombined 3 regimes into one program.
  4. Final simplification13.4

    \[\leadsto \begin{array}{l} \mathbf{if}\;V \cdot \ell \le -1.125220159166195 \cdot 10^{-80}:\\ \;\;\;\;\sqrt{\sqrt{\frac{A}{V \cdot \ell}}} \cdot \left(c0 \cdot \sqrt{\sqrt{A \cdot \frac{1}{V \cdot \ell}}}\right)\\ \mathbf{elif}\;V \cdot \ell \le 2.2677613144113 \cdot 10^{-321}:\\ \;\;\;\;c0 \cdot \sqrt{\frac{\frac{A}{\ell}}{V}}\\ \mathbf{else}:\\ \;\;\;\;\frac{\sqrt{A}}{\sqrt{V \cdot \ell}} \cdot c0\\ \end{array}\]

Reproduce

herbie shell --seed 2018365 +o rules:numerics
(FPCore (c0 A V l)
  :name "Henrywood and Agarwal, Equation (3)"
  (* c0 (sqrt (/ A (* V l)))))

Details

Time bar (total: 9.5s)Debug log

start72.0ms

Algorithm
intervals

setup33.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 14.7b

localize28.0ms

Local error

Found 3 expressions with local error:

16.1b
(sqrt (/ A (* V l)))
3.6b
(/ A (* V l))
2.3b
(* c0 (sqrt (/ A (* V l))))

rewrite7.0ms

Algorithm
rewrite-expression-head
Counts
3 → 47
Calls

3 calls. Slowest were:

4.0ms
(* c0 (sqrt (/ A (* V l))))
2.0ms
(sqrt (/ A (* V l)))
1.0ms
(/ A (* V l))

series116.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

61.0ms
(* c0 (sqrt (/ A (* V l))))
33.0ms
(sqrt (/ A (* V l)))
22.0ms
(/ A (* V l))

simplify79.0ms

Counts
23 → 56
Calls

23 calls. Slowest were:

6.0ms
(sqrt (* (cbrt (/ A (* V l))) (cbrt (/ A (* V l)))))
6.0ms
(sqrt (cbrt (/ A (* V l))))
6.0ms
(sqrt (sqrt (/ A (* V l))))

prune593.0ms

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 12.4b

localize12.0ms

Local error

Found 4 expressions with local error:

16.1b
(sqrt (* A (/ 1 (* V l))))
3.7b
(* A (/ 1 (* V l)))
2.3b
(* c0 (sqrt (* A (/ 1 (* V l)))))
0.2b
(/ 1 (* V l))

rewrite5.0ms

Algorithm
rewrite-expression-head
Counts
4 → 61
Calls

4 calls. Slowest were:

2.0ms
(* c0 (sqrt (* A (/ 1 (* V l)))))
1.0ms
(sqrt (* A (/ 1 (* V l))))
1.0ms
(* A (/ 1 (* V l)))

series111.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

56.0ms
(* c0 (sqrt (* A (/ 1 (* V l)))))
23.0ms
(* A (/ 1 (* V l)))
20.0ms
(sqrt (* A (/ 1 (* V l))))
11.0ms
(/ 1 (* V l))

simplify232.0ms

Counts
28 → 73
Calls

28 calls. Slowest were:

79.0ms
(sqrt (sqrt (* A (/ 1 (* V l)))))
47.0ms
(cbrt (sqrt (* A (/ 1 (* V l)))))
39.0ms
(sqrt (* A (/ 1 (* V l))))

prune959.0ms

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 12.4b

localize13.0ms

Local error

Found 4 expressions with local error:

16.1b
(sqrt (* A (/ 1 (* V l))))
16.1b
(sqrt (/ A (* V l)))
3.7b
(* A (/ 1 (* V l)))
3.6b
(/ A (* V l))

rewrite4.0ms

Algorithm
rewrite-expression-head
Counts
4 → 58
Calls

4 calls. Slowest were:

1.0ms
(sqrt (* A (/ 1 (* V l))))
1.0ms
(* A (/ 1 (* V l)))
1.0ms
(sqrt (/ A (* V l)))

series83.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

24.0ms
(sqrt (/ A (* V l)))
23.0ms
(/ A (* V l))
23.0ms
(sqrt (* A (/ 1 (* V l))))
13.0ms
(* A (/ 1 (* V l)))

simplify80.0ms

Counts
27 → 70
Calls

27 calls. Slowest were:

7.0ms
(sqrt (* (cbrt (/ A (* V l))) (cbrt (/ A (* V l)))))
6.0ms
(sqrt (/ 1 (* V l)))
6.0ms
(cbrt (/ 1 (* V l)))

prune1.1s

Pruning

13 alts after pruning (11 fresh and 2 done)

Merged error: 12.4b

localize34.0ms

Local error

Found 4 expressions with local error:

20.0b
(cbrt (/ A (* V l)))
20.0b
(cbrt (/ A (* V l)))
20.0b
(cbrt (/ A (* V l)))
10.1b
(sqrt (* (cbrt (/ A (* V l))) (cbrt (/ A (* V l)))))

rewrite7.0ms

Algorithm
rewrite-expression-head
Counts
4 → 60
Calls

4 calls. Slowest were:

4.0ms
(sqrt (* (cbrt (/ A (* V l))) (cbrt (/ A (* V l)))))
1.0ms
(cbrt (/ A (* V l)))
1.0ms
(cbrt (/ A (* V l)))

series400.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

121.0ms
(cbrt (/ A (* V l)))
109.0ms
(cbrt (/ A (* V l)))
85.0ms
(sqrt (* (cbrt (/ A (* V l))) (cbrt (/ A (* V l)))))
85.0ms
(cbrt (/ A (* V l)))

simplify3.4s

Counts
32 → 72
Calls

32 calls. Slowest were:

510.0ms
(exp (* 1/3 (- (+ (log (/ 1 l)) (log (/ 1 V))) (log (/ 1 A)))))
469.0ms
(exp (* 1/3 (- (+ (log (/ 1 l)) (log (/ 1 V))) (log (/ 1 A)))))
391.0ms
(exp (* 1/3 (- (+ (log (/ 1 l)) (log (/ 1 V))) (log (/ 1 A)))))

prune1.0s

Pruning

11 alts after pruning (10 fresh and 1 done)

Merged error: 12.4b

regimes1.1s

Accuracy

64% (3.2b remaining)

Error of 13.4b against oracle of 10.3b and baseline of 19.0b

bsearch8.0ms