Average Error: 47.4 → 19.6
Time: 1.5m
Precision: 64
Internal Precision: 128
\[\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
\[\begin{array}{l} \mathbf{if}\;\ell \cdot \ell \le 1.1997133799852223 \cdot 10^{+295}:\\ \;\;\;\;\frac{2}{\frac{1}{\frac{{\ell}^{2} \cdot \cos k}{\left(\left(k \cdot t\right) \cdot k\right) \cdot {\left(\sin k\right)}^{2}}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{\left(\tan k \cdot \left(\left(\frac{t \cdot t}{\ell} \cdot \frac{t}{\ell}\right) \cdot \sin k\right)\right) \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + 1\right) - 1\right)}\\ \end{array}\]

Error

Bits error versus t

Bits error versus l

Bits error versus k

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if (* l l) < 1.1997133799852223e+295

    1. Initial program 44.6

      \[\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
    2. Taylor expanded around -inf 17.2

      \[\leadsto \frac{2}{\color{blue}{\frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{{\ell}^{2} \cdot \cos k}}}\]
    3. Using strategy rm
    4. Applied associate-*r*15.7

      \[\leadsto \frac{2}{\frac{\color{blue}{\left(t \cdot {k}^{2}\right) \cdot {\left(\sin k\right)}^{2}}}{{\ell}^{2} \cdot \cos k}}\]
    5. Using strategy rm
    6. Applied unpow215.7

      \[\leadsto \frac{2}{\frac{\left(t \cdot \color{blue}{\left(k \cdot k\right)}\right) \cdot {\left(\sin k\right)}^{2}}{{\ell}^{2} \cdot \cos k}}\]
    7. Applied associate-*r*13.2

      \[\leadsto \frac{2}{\frac{\color{blue}{\left(\left(t \cdot k\right) \cdot k\right)} \cdot {\left(\sin k\right)}^{2}}{{\ell}^{2} \cdot \cos k}}\]
    8. Using strategy rm
    9. Applied clear-num13.2

      \[\leadsto \frac{2}{\color{blue}{\frac{1}{\frac{{\ell}^{2} \cdot \cos k}{\left(\left(t \cdot k\right) \cdot k\right) \cdot {\left(\sin k\right)}^{2}}}}}\]

    if 1.1997133799852223e+295 < (* l l)

    1. Initial program 61.2

      \[\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
    2. Using strategy rm
    3. Applied add-cube-cbrt61.2

      \[\leadsto \frac{2}{\left(\left(\frac{\color{blue}{\left(\sqrt[3]{{t}^{3}} \cdot \sqrt[3]{{t}^{3}}\right) \cdot \sqrt[3]{{t}^{3}}}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
    4. Applied times-frac54.9

      \[\leadsto \frac{2}{\left(\left(\color{blue}{\left(\frac{\sqrt[3]{{t}^{3}} \cdot \sqrt[3]{{t}^{3}}}{\ell} \cdot \frac{\sqrt[3]{{t}^{3}}}{\ell}\right)} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
    5. Simplified54.9

      \[\leadsto \frac{2}{\left(\left(\left(\color{blue}{\frac{t \cdot t}{\ell}} \cdot \frac{\sqrt[3]{{t}^{3}}}{\ell}\right) \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
    6. Simplified51.5

      \[\leadsto \frac{2}{\left(\left(\left(\frac{t \cdot t}{\ell} \cdot \color{blue}{\frac{t}{\ell}}\right) \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification19.6

    \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \cdot \ell \le 1.1997133799852223 \cdot 10^{+295}:\\ \;\;\;\;\frac{2}{\frac{1}{\frac{{\ell}^{2} \cdot \cos k}{\left(\left(k \cdot t\right) \cdot k\right) \cdot {\left(\sin k\right)}^{2}}}}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{\left(\tan k \cdot \left(\left(\frac{t \cdot t}{\ell} \cdot \frac{t}{\ell}\right) \cdot \sin k\right)\right) \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + 1\right) - 1\right)}\\ \end{array}\]

Reproduce

herbie shell --seed 2019007 
(FPCore (t l k)
  :name "Toniolo and Linder, Equation (10-)"
  (/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))))

Details

Time bar (total: 1.5m)Debug log

sample1.2s

Algorithm
intervals

simplify161.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

161.0ms
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1)))

prune19.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 47.2b

localize92.0ms

Local error

Found 4 expressions with local error:

33.8b
(* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))
14.3b
(/ (pow t 3) (* l l))
1.3b
(- (+ 1 (pow (/ k t) 2)) 1)
1.0b
(* (* (/ (pow t 3) (* l l)) (sin k)) (tan k))

rewrite124.0ms

Algorithm
rewrite-expression-head
Counts
4 → 103
Calls

4 calls. Slowest were:

74.0ms
(* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))
27.0ms
(* (* (/ (pow t 3) (* l l)) (sin k)) (tan k))
15.0ms
(- (+ 1 (pow (/ k t) 2)) 1)

series586.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

387.0ms
(* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))
155.0ms
(* (* (/ (pow t 3) (* l l)) (sin k)) (tan k))
29.0ms
(/ (pow t 3) (* l l))
15.0ms
(- (+ 1 (pow (/ k t) 2)) 1)

simplify18.0s

Counts
86 → 115
Calls

86 calls. Slowest were:

950.0ms
(* (cos k) (+ (* (+ 1 (pow (/ k t) 2)) (+ 1 (pow (/ k t) 2))) (+ (* 1 1) (* (+ 1 (pow (/ k t) 2)) 1))))
924.0ms
(* (* (* (pow t 3) (sin k)) (sin k)) (- (+ 1 (pow (/ k t) 2)) 1))
837.0ms
(* (* (* (pow t 3) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))

prune2.0s

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 18.4b

localize31.0ms

Local error

Found 4 expressions with local error:

18.6b
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
10.3b
(* t (* (pow k 2) (pow (sin k) 2)))
0.5b
(/ 2 (/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k))))
0.4b
(pow (sin k) 2)

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
4 → 57
Calls

4 calls. Slowest were:

9.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
8.0ms
(* t (* (pow k 2) (pow (sin k) 2)))
1.0ms
(pow (sin k) 2)

series592.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

252.0ms
(/ 2 (/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k))))
222.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
93.0ms
(* t (* (pow k 2) (pow (sin k) 2)))
25.0ms
(pow (sin k) 2)

simplify5.7s

Counts
32 → 69
Calls

32 calls. Slowest were:

560.0ms
(/ (* (* (* t (* (pow k 2) (pow (sin k) 2))) (* t (* (pow k 2) (pow (sin k) 2)))) (* t (* (pow k 2) (pow (sin k) 2)))) (* (* (* (pow l 2) (pow l 2)) (pow l 2)) (* (* (cos k) (cos k)) (cos k))))
529.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
528.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))

prune1.4s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 15.0b

localize26.0ms

Local error

Found 4 expressions with local error:

18.6b
(/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k)))
9.2b
(* t (pow k 2))
0.5b
(/ 2 (/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k))))
0.4b
(pow (sin k) 2)

rewrite28.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

20.0ms
(/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k)))
2.0ms
(/ 2 (/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k))))
2.0ms
(* t (pow k 2))

series578.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

287.0ms
(/ 2 (/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k))))
248.0ms
(/ (* (* t (pow k 2)) (pow (sin k) 2)) (* (pow l 2) (cos k)))
26.0ms
(* t (pow k 2))
16.0ms
(pow (sin k) 2)

simplify7.1s

Counts
46 → 83
Calls

46 calls. Slowest were:

674.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
631.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
593.0ms
(- (+ (log (* t (pow k 2))) (* (log (sin k)) 2)) (log (* (pow l 2) (cos k))))

prune1.3s

Pruning

15 alts after pruning (15 fresh and 0 done)

Merged error: 12.6b

localize28.0ms

Local error

Found 4 expressions with local error:

18.6b
(/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k)))
0.5b
(/ 2 (/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k))))
0.4b
(pow (sin k) 2)
0.3b
(* (* (* t k) k) (pow (sin k) 2))

rewrite50.0ms

Algorithm
rewrite-expression-head
Counts
4 → 79
Calls

4 calls. Slowest were:

24.0ms
(/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k)))
20.0ms
(* (* (* t k) k) (pow (sin k) 2))
3.0ms
(/ 2 (/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k))))

series624.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

301.0ms
(/ 2 (/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k))))
235.0ms
(/ (* (* (* t k) k) (pow (sin k) 2)) (* (pow l 2) (cos k)))
66.0ms
(* (* (* t k) k) (pow (sin k) 2))
22.0ms
(pow (sin k) 2)

simplify7.0s

Counts
55 → 91
Calls

55 calls. Slowest were:

558.0ms
(/ (* t (* (pow k 2) (pow (sin k) 2))) (* (pow l 2) (cos k)))
537.0ms
(/ (* (* (* (* (* t k) k) (* (* t k) k)) (* (* t k) k)) (* (* (pow (sin k) 2) (pow (sin k) 2)) (pow (sin k) 2))) (* (* (* (pow l 2) (pow l 2)) (pow l 2)) (* (* (cos k) (cos k)) (cos k))))
528.0ms
(/ (* (* (* (* (* t k) k) (pow (sin k) 2)) (* (* (* t k) k) (pow (sin k) 2))) (* (* (* t k) k) (pow (sin k) 2))) (* (* (* (pow l 2) (cos k)) (* (pow l 2) (cos k))) (* (pow l 2) (cos k))))

prune1.4s

Pruning

17 alts after pruning (17 fresh and 0 done)

Merged error: 11.9b

regimes531.0ms

Accuracy

16.9% (6.7b remaining)

Error of 19.6b against oracle of 13.0b and baseline of 21.0b

bsearch7.0ms

end0.0ms

sample40.9s

Algorithm
intervals