


Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus Omc
Results
if (/ t l) < -4.0647357290577327e+34 or 1.6327068826043366e+85 < (/ t l) Initial program 23.4
rmApplied expm1-log1p-u23.4
rmApplied *-un-lft-identity23.4
Applied add-cube-cbrt23.4
Applied times-frac23.4
Applied sqrt-prod23.4
Simplified23.4
rmApplied add-sqr-sqrt23.5
Applied add-sqr-sqrt23.5
Applied cbrt-prod23.5
Applied times-frac23.5
Applied rem-sqrt-square23.5
Taylor expanded around -inf 0.9
if -4.0647357290577327e+34 < (/ t l) < 1.6327068826043366e+85Initial program 1.1
rmApplied expm1-log1p-u1.1
rmApplied *-un-lft-identity1.1
Applied add-cube-cbrt1.1
Applied times-frac1.1
Applied sqrt-prod1.1
Simplified1.1
rmApplied add-sqr-sqrt1.2
Applied add-sqr-sqrt1.2
Applied cbrt-prod1.2
Applied times-frac1.2
Applied rem-sqrt-square1.2
Final simplification1.1
herbie shell --seed 2019007 +o rules:numerics
(FPCore (t l Om Omc)
:name "Toniolo and Linder, Equation (2)"
(asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))))
Time bar (total: 37.3s)Debug log
1 calls. Slowest were:
| 228.0ms | (asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 10.4b
Found 4 expressions with local error:
| 9.7b | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 0.3b | (asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
| 0.3b | (pow (/ t l) 2) |
| 0.3b | (pow (/ Om Omc) 2) |
4 calls. Slowest were:
| 4.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 1.0ms | (pow (/ Om Omc) 2) |
| 1.0ms | (pow (/ t l) 2) |
4 calls. Slowest were:
| 106.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 20.0ms | (asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
| 17.0ms | (pow (/ Om Omc) 2) |
| 17.0ms | (pow (/ t l) 2) |
42 calls. Slowest were:
| 209.0ms | (sqrt (* (cbrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) (cbrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))))) |
| 198.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ (pow 1 3) (pow (* 2 (pow (/ t l) 2)) 3)))) |
| 193.0ms | (sqrt (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 10.4b
Found 4 expressions with local error:
| 9.7b | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 0.3b | (asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
| 0.3b | (pow (/ t l) 2) |
| 0.3b | (pow (/ Om Omc) 2) |
4 calls. Slowest were:
| 3.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 1.0ms | (pow (/ t l) 2) |
| 1.0ms | (pow (/ Om Omc) 2) |
4 calls. Slowest were:
| 111.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 17.0ms | (asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
| 17.0ms | (pow (/ Om Omc) 2) |
| 16.0ms | (pow (/ t l) 2) |
42 calls. Slowest were:
| 253.0ms | (sqrt (* (cbrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))) (cbrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2))))))) |
| 248.0ms | (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ (pow 1 3) (pow (* 2 (pow (/ t l) 2)) 3)))) |
| 216.0ms | (sqrt (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 10.4b
Found 4 expressions with local error:
| 9.7b | (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 0.3b | (asin (* (fabs (cbrt (- 1 (* (/ Om Omc) (/ Om Omc))))) (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (+ 1 (* 2 (pow (/ t l) 2))))))) |
| 0.3b | (pow (/ t l) 2) |
| 0.3b | (* (/ Om Omc) (/ Om Omc)) |
4 calls. Slowest were:
| 13.0ms | (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 7.0ms | (* (/ Om Omc) (/ Om Omc)) |
| 2.0ms | (pow (/ t l) 2) |
4 calls. Slowest were:
| 291.0ms | (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (+ 1 (* 2 (pow (/ t l) 2))))) |
| 28.0ms | (asin (* (fabs (cbrt (- 1 (* (/ Om Omc) (/ Om Omc))))) (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (+ 1 (* 2 (pow (/ t l) 2))))))) |
| 15.0ms | (pow (/ t l) 2) |
| 14.0ms | (* (/ Om Omc) (/ Om Omc)) |
61 calls. Slowest were:
| 880.0ms | (* (* (* (/ Om Omc) (/ Om Omc)) (/ Om Omc)) (* (* (/ Om Omc) (/ Om Omc)) (/ Om Omc))) |
| 169.0ms | (sqrt (/ (cbrt (- 1 (pow (/ Om Omc) 2))) (- (* 1 1) (* (* 2 (pow (/ t l) 2)) (* 2 (pow (/ t l) 2)))))) |
| 143.0ms | (asin (* (pow (- 1 (/ (pow Om 2) (pow Omc 2))) 1/6) (* (sqrt (/ 1 (+ (* 2 (/ (pow t 2) (pow l 2))) 1))) (fabs (pow (- 1 (/ (pow Om 2) (pow Omc 2))) 1/3))))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 10.4b
Found 4 expressions with local error:
| 9.6b | (sqrt (+ 1 (* 2 (pow (/ t l) 2)))) |
| 0.6b | (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2))))) |
| 0.3b | (asin (* (fabs (cbrt (- 1 (* (/ Om Omc) (/ Om Omc))))) (fabs (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2)))))))) |
| 0.3b | (pow (/ t l) 2) |
4 calls. Slowest were:
| 7.0ms | (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2))))) |
| 1.0ms | (pow (/ t l) 2) |
| 1.0ms | (asin (* (fabs (cbrt (- 1 (* (/ Om Omc) (/ Om Omc))))) (fabs (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2)))))))) |
4 calls. Slowest were:
| 217.0ms | (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2))))) |
| 36.0ms | (sqrt (+ 1 (* 2 (pow (/ t l) 2)))) |
| 18.0ms | (asin (* (fabs (cbrt (- 1 (* (/ Om Omc) (/ Om Omc))))) (fabs (/ (cbrt (sqrt (- 1 (pow (/ Om Omc) 2)))) (sqrt (+ 1 (* 2 (pow (/ t l) 2)))))))) |
| 10.0ms | (pow (/ t l) 2) |
152 calls. Slowest were:
| 809.0ms | (cbrt (sqrt (- (* 1 1) (* (pow (/ Om Omc) 2) (pow (/ Om Omc) 2))))) |
| 295.0ms | (cbrt (sqrt (- (pow 1 3) (pow (pow (/ Om Omc) 2) 3)))) |
| 224.0ms | (asin (* (fabs (* (pow (- 1 (/ (pow Om 2) (pow Omc 2))) 1/6) (sqrt (/ 1 (+ (* 2 (/ (pow t 2) (pow l 2))) 1))))) (fabs (pow (- 1 (/ (pow Om 2) (pow Omc 2))) 1/3)))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 0.9b
98.4% (0.2b remaining)
Error of 1.1b against oracle of 0.9b and baseline of 10.4b