


Bits error versus l



Bits error versus Om



Bits error versus kx



Bits error versus ky
Initial program 1.6
Initial simplification1.6
Taylor expanded around inf 16.8
Simplified0.6
Final simplification0.6
herbie shell --seed 2018365 +o rules:numerics
(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))))))))))
Time bar (total: 25.2s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 1.6b
Found 4 expressions with local error:
| 9.1b | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 2.4b | (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1) |
| 0.4b | (/ 1/2 (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1))) |
| 0.3b | (* (sin kx) (sin kx)) |
4 calls. Slowest were:
| 3.0ms | (* (sin kx) (sin kx)) |
| 1.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 1.0ms | (/ 1/2 (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1))) |
4 calls. Slowest were:
| 307.0ms | (/ 1/2 (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1))) |
| 259.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 137.0ms | (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1) |
| 25.0ms | (* (sin kx) (sin kx)) |
35 calls. Slowest were:
| 348.0ms | (+ (* 4 (/ (* (pow l 2) (pow (sin ky) 2)) (pow Om 2))) (+ (* 4 (/ (* (pow (sin kx) 2) (pow l 2)) (pow Om 2))) 1)) |
| 344.0ms | (+ (* 4 (/ (* (pow l 2) (pow (sin ky) 2)) (pow Om 2))) (+ (* 4 (/ (* (pow (sin kx) 2) (pow l 2)) (pow Om 2))) 1)) |
| 296.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.7b
Found 4 expressions with local error:
| 9.1b | (sqrt (fma 4 (fma (/ (sin kx) (/ Om l)) (/ (sin kx) (/ Om l)) (* (/ (sin ky) (/ Om l)) (/ (sin ky) (/ Om l)))) 1)) |
| 3.9b | (/ (sin ky) (/ Om l)) |
| 3.9b | (/ (sin ky) (/ Om l)) |
| 2.3b | (/ (sin kx) (/ Om l)) |
4 calls. Slowest were:
| 6.0ms | (/ (sin ky) (/ Om l)) |
| 6.0ms | (/ (sin ky) (/ Om l)) |
| 6.0ms | (/ (sin kx) (/ Om l)) |
4 calls. Slowest were:
| 286.0ms | (sqrt (fma 4 (fma (/ (sin kx) (/ Om l)) (/ (sin kx) (/ Om l)) (* (/ (sin ky) (/ Om l)) (/ (sin ky) (/ Om l)))) 1)) |
| 93.0ms | (/ (sin ky) (/ Om l)) |
| 86.0ms | (/ (sin ky) (/ Om l)) |
| 62.0ms | (/ (sin kx) (/ Om l)) |
112 calls. Slowest were:
| 264.0ms | (- (/ (* ky l) Om) (* 1/6 (/ (* (pow ky 3) l) Om))) |
| 243.0ms | (- (/ (* ky l) Om) (* 1/6 (/ (* (pow ky 3) l) Om))) |
| 222.0ms | (/ (* (* (sin ky) (sin ky)) (sin ky)) (* (* (/ Om l) (/ Om l)) (/ Om l))) |
3 alts after pruning (2 fresh and 1 done)
Merged error: 0.7b
Found 4 expressions with local error:
| 3.9b | (/ (sin ky) (/ Om l)) |
| 3.9b | (/ (sin ky) (/ Om l)) |
| 2.3b | (/ (sin kx) (/ Om l)) |
| 2.3b | (/ (sin kx) (/ Om l)) |
4 calls. Slowest were:
| 6.0ms | (/ (sin ky) (/ Om l)) |
| 6.0ms | (/ (sin kx) (/ Om l)) |
| 6.0ms | (/ (sin ky) (/ Om l)) |
4 calls. Slowest were:
| 72.0ms | (/ (sin kx) (/ Om l)) |
| 68.0ms | (/ (sin ky) (/ Om l)) |
| 61.0ms | (/ (sin kx) (/ Om l)) |
| 58.0ms | (/ (sin ky) (/ Om l)) |
140 calls. Slowest were:
| 273.0ms | (/ (* (* (sin kx) (sin kx)) (sin kx)) (* (* (/ Om l) (/ Om l)) (/ Om l))) |
| 262.0ms | (- (/ (* ky l) Om) (* 1/6 (/ (* (pow ky 3) l) Om))) |
| 245.0ms | (/ (* (* (sin kx) (sin kx)) (sin kx)) (* (* (/ Om l) (/ Om l)) (/ Om l))) |
3 alts after pruning (1 fresh and 2 done)
Merged error: 0.7b
Found 4 expressions with local error:
| 9.1b | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 9.1b | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 9.1b | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 8.4b | (cbrt (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1))) |
4 calls. Slowest were:
| 1.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 1.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 1.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
4 calls. Slowest were:
| 321.0ms | (cbrt (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1))) |
| 264.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 248.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
| 224.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
30 calls. Slowest were:
| 278.0ms | (sqrt (* (cbrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) (cbrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)))) |
| 272.0ms | (cbrt (sqrt (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)))) |
| 260.0ms | (sqrt (fma (* (* l (/ 2 Om)) (* l (/ 2 Om))) (fma (sin ky) (sin ky) (* (sin kx) (sin kx))) 1)) |
3 alts after pruning (0 fresh and 3 done)
Merged error: 0.7b
0% (0.0b remaining)
Error of 0.6b against oracle of 0.6b and baseline of 0.6b