Average Error: 1.5 → 1.5
Time: 42.4s
Precision: 64
Internal Precision: 128
\[\sqrt{\frac{1}{2} \cdot \left(1 + \frac{1}{\sqrt{1 + {\left(\frac{2 \cdot \ell}{Om}\right)}^{2} \cdot \left({\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}\right)}}\right)}\]
\[\sqrt{\frac{1}{2} + \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} \cdot \left(\sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} \cdot \left(\sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} \cdot \sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}\right)}}\right)}\]

Error

Bits error versus l

Bits error versus Om

Bits error versus kx

Bits error versus ky

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 1.5

    \[\sqrt{\frac{1}{2} \cdot \left(1 + \frac{1}{\sqrt{1 + {\left(\frac{2 \cdot \ell}{Om}\right)}^{2} \cdot \left({\left(\sin kx\right)}^{2} + {\left(\sin ky\right)}^{2}\right)}}\right)}\]
  2. Initial simplification1.5

    \[\leadsto \sqrt{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} + \frac{1}{2}}\]
  3. Using strategy rm
  4. Applied add-cube-cbrt1.5

    \[\leadsto \sqrt{\color{blue}{\left(\sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}}\right) \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}}} + \frac{1}{2}}\]
  5. Using strategy rm
  6. Applied add-cube-cbrt1.5

    \[\leadsto \sqrt{\left(\sqrt[3]{\frac{\frac{1}{2}}{\color{blue}{\left(\sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} \cdot \sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}\right) \cdot \sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}}}} \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}}\right) \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} + \frac{1}{2}}\]
  7. Final simplification1.5

    \[\leadsto \sqrt{\frac{1}{2} + \sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} \cdot \left(\sqrt[3]{\frac{\frac{1}{2}}{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}} \cdot \sqrt[3]{\frac{\frac{1}{2}}{\sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} \cdot \left(\sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}} \cdot \sqrt[3]{\sqrt{\left(\sin ky \cdot \sin ky + \sin kx \cdot \sin kx\right) \cdot \left(\frac{\ell}{\frac{Om}{2}} \cdot \frac{\ell}{\frac{Om}{2}}\right) + 1}}\right)}}\right)}\]

Reproduce

herbie shell --seed 2018362 
(FPCore (l Om kx ky)
  :name "Toniolo and Linder, Equation (3a)"
  (sqrt (* (/ 1 2) (+ 1 (/ 1 (sqrt (+ 1 (* (pow (/ (* 2 l) Om) 2) (+ (pow (sin kx) 2) (pow (sin ky) 2))))))))))

Details

Time bar (total: 36.0s)Debug log

start158.0ms

Algorithm
intervals

setup332.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 1.1b

localize81.0ms

Local error

Found 4 expressions with local error:

9.1b
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
3.3b
(* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2))))
1.0b
(/ 1/2 (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
0.2b
(* (sin ky) (sin ky))

rewrite93.0ms

Algorithm
rewrite-expression-head
Counts
4 → 87
Calls

4 calls. Slowest were:

40.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
35.0ms
(* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2))))
12.0ms
(/ 1/2 (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))

series934.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

371.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
314.0ms
(/ 1/2 (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
227.0ms
(* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2))))
22.0ms
(* (sin ky) (sin ky))

simplify15.9s

Counts
72 → 99
Calls

72 calls. Slowest were:

2.5s
(* (+ (* (- (cos (- ky ky)) (cos (+ ky ky))) 2) (* 2 (- (cos (- kx kx)) (cos (+ kx kx))))) (* (/ l (/ Om 2)) (/ l (/ Om 2))))
1.8s
(* (+ (* (- (cos (- ky ky)) (cos (+ ky ky))) 2) (* 2 (- (cos (- kx kx)) (cos (+ kx kx))))) (* (/ l (/ Om 2)) l))
1.5s
(* (+ (* (- (cos (- ky ky)) (cos (+ ky ky))) 2) (* 2 (- (cos (- kx kx)) (cos (+ kx kx))))) (* l (/ l (/ Om 2))))

prune1.9s

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.6b

localize18.0ms

Local error

Found 4 expressions with local error:

9.1b
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
9.1b
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
9.1b
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
8.9b
(cbrt (/ 1/2 (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))))

rewrite59.0ms

Algorithm
rewrite-expression-head
Counts
4 → 57
Calls

4 calls. Slowest were:

18.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
16.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
16.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))

series988.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

266.0ms
(cbrt (/ 1/2 (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))))
241.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
240.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
240.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))

simplify3.3s

Counts
37 → 69
Calls

37 calls. Slowest were:

234.0ms
(sqrt (* (cbrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)) (cbrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))))
232.0ms
(sqrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
217.0ms
(sqrt (* (cbrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)) (cbrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))))

prune1.9s

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.6b

localize56.0ms

Local error

Found 4 expressions with local error:

10.0b
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
10.0b
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
10.0b
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
9.1b
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))

rewrite77.0ms

Algorithm
rewrite-expression-head
Counts
4 → 62
Calls

4 calls. Slowest were:

28.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
21.0ms
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
15.0ms
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))

series1.1s

Counts
4 → 12
Calls

4 calls. Slowest were:

308.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
297.0ms
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
270.0ms
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))
241.0ms
(cbrt (sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1)))

simplify6.6s

Counts
42 → 74
Calls

42 calls. Slowest were:

578.0ms
(exp (* 1/6 (- (+ (log (+ (* 4 (pow (sin ky) 2)) (* 4 (pow (sin kx) 2)))) (* 2 (log (/ -1 Om)))) (* 2 (log (/ -1 l))))))
306.0ms
(sqrt (+ (* (+ (* (sin ky) (sin ky)) (* (sin kx) (sin kx))) (* (/ l (/ Om 2)) (/ l (/ Om 2)))) 1))
303.0ms
(exp (* 1/6 (- (+ (* 2 (log (/ 1 Om))) (log (+ (* 4 (pow (sin ky) 2)) (* 4 (pow (sin kx) 2))))) (* 2 (log (/ 1 l))))))

prune2.3s

Pruning

2 alts after pruning (1 fresh and 1 done)

Merged error: 0.6b

localize10.0ms

Local error

Found 0 expressions with local error:

rewrite0.0ms

Algorithm
rewrite-expression-head
Counts
0 → 0
Calls

0 calls. Slowest were:

series0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

simplify0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

prune1.0ms

Pruning

2 alts after pruning (0 fresh and 2 done)

Merged error: 0.6b

regimes317.0ms

Accuracy

0% (0.7b remaining)

Error of 1.5b against oracle of 0.8b and baseline of 1.5b

bsearch4.0ms