


Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus Omc
Results
if (/ t l) < 2.1709653990804656e+23Initial program 7.2
Initial simplification7.2
rmApplied add-sqr-sqrt7.8
rmApplied add-sqr-sqrt7.3
if 2.1709653990804656e+23 < (/ t l) Initial program 21.4
Initial simplification21.4
rmApplied sqrt-div21.4
Taylor expanded around -inf 1.1
Final simplification5.8
herbie shell --seed 2018365
(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: 57.4s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 12.0b
Found 4 expressions with local error:
| 10.9b | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 0.5b | (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 0.4b | (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 0.2b | (* (/ t l) (/ t l)) |
4 calls. Slowest were:
| 4.0ms | (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 3.0ms | (* (/ t l) (/ t l)) |
| 2.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
4 calls. Slowest were:
| 90.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 51.0ms | (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 14.0ms | (* (/ t l) (/ t l)) |
| 11.0ms | (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
75 calls. Slowest were:
| 910.0ms | (* (* (* (/ t l) (/ t l)) (/ t l)) (* (* (/ t l) (/ t l)) (/ t l))) |
| 587.0ms | (+ (* 1 1) (- (* (* 2 (* (/ t l) (/ t l))) (* 2 (* (/ t l) (/ t l)))) (* 1 (* 2 (* (/ t l) (/ t l)))))) |
| 492.0ms | (sqrt (/ 1 (* (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 11.9b
Found 4 expressions with local error:
| 11.0b | (sqrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 0.6b | (/ (sqrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (sqrt (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 0.5b | (asin (/ (sqrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (sqrt (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 0.2b | (* (/ t l) (/ t l)) |
4 calls. Slowest were:
| 5.0ms | (/ (sqrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (sqrt (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 3.0ms | (* (/ t l) (/ t l)) |
| 1.0ms | (sqrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
4 calls. Slowest were:
| 123.0ms | (/ (sqrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (sqrt (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 57.0ms | (sqrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 16.0ms | (asin (/ (sqrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (sqrt (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 14.0ms | (* (/ t l) (/ t l)) |
123 calls. Slowest were:
| 1.0s | (* (* (* (/ t l) (/ t l)) (/ t l)) (* (* (/ t l) (/ t l)) (/ t l))) |
| 526.0ms | (asin (sqrt (/ (- 1 (/ (pow Om 2) (pow Omc 2))) (+ (* 2 (/ (pow t 2) (pow l 2))) 1)))) |
| 487.0ms | (sqrt (+ (* 1 1) (- (* (* 2 (* (/ t l) (/ t l))) (* 2 (* (/ t l) (/ t l)))) (* 1 (* 2 (* (/ t l) (/ t l))))))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 7.3b
Found 4 expressions with local error:
| 18.3b | (cbrt (+ (* (* (/ t l) (/ t l)) 2) 1)) |
| 18.3b | (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 2.0b | (sqrt (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 0.5b | (asin (* (fabs (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ (* (* (/ t l) (/ t l)) 2) 1)))) (sqrt (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))))))) |
4 calls. Slowest were:
| 17.0ms | (cbrt (+ (* (* (/ t l) (/ t l)) 2) 1)) |
| 5.0ms | (sqrt (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 1.0ms | (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
4 calls. Slowest were:
| 216.0ms | (sqrt (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))))) |
| 140.0ms | (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))) |
| 134.0ms | (cbrt (+ (* (* (/ t l) (/ t l)) 2) 1)) |
| 39.0ms | (asin (* (fabs (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ (* (* (/ t l) (/ t l)) 2) 1)))) (sqrt (/ (cbrt (- 1 (* (/ Om Omc) (/ Om Omc)))) (cbrt (+ 1 (* 2 (* (/ t l) (/ t l))))))))) |
70 calls. Slowest were:
| 830.0ms | (exp (* 1/6 (- (+ (* 2 (log (/ 1 t))) (+ (* 2 (log (/ 1 Omc))) (log -1/2))) (+ (* 2 (log (/ 1 Om))) (* 2 (log (/ 1 l))))))) |
| 473.0ms | (asin (* (pow (/ (- 1 (/ (pow Om 2) (pow Omc 2))) (+ (* 2 (/ (pow t 2) (pow l 2))) 1)) 1/6) (fabs (pow (/ (- 1 (/ (pow Om 2) (pow Omc 2))) (+ (* 2 (/ (pow t 2) (pow l 2))) 1)) 1/3)))) |
| 465.0ms | (cbrt (+ (pow (* (* (/ t l) (/ t l)) 2) 3) (pow 1 3))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 7.3b
Found 4 expressions with local error:
| 10.9b | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 10.9b | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 9.0b | (sqrt (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
| 9.0b | (sqrt (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
4 calls. Slowest were:
| 4.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 4.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 3.0ms | (sqrt (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
4 calls. Slowest were:
| 112.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 111.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))) |
| 32.0ms | (sqrt (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
| 30.0ms | (sqrt (asin (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ 1 (* 2 (* (/ t l) (/ t l)))))))) |
56 calls. Slowest were:
| 549.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ (pow 1 3) (pow (* 2 (* (/ t l) (/ t l))) 3)))) |
| 517.0ms | (sqrt (/ (- 1 (* (/ Om Omc) (/ Om Omc))) (+ (pow 1 3) (pow (* 2 (* (/ t l) (/ t l))) 3)))) |
| 485.0ms | (sqrt (asin (sqrt (/ (- 1 (/ (pow Om 2) (pow Omc 2))) (+ (* 2 (/ (pow t 2) (pow l 2))) 1))))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 7.3b
96.6% (0.2b remaining)
Error of 5.8b against oracle of 5.7b and baseline of 10.6b