


Bits error versus re



Bits error versus im
Results
if re < -5.087015824912348e+107Initial program 51.8
rmApplied add-sqr-sqrt51.8
Applied pow151.8
Applied log-pow51.8
Applied times-frac51.7
rmApplied div-inv51.7
Taylor expanded around -inf 8.8
Simplified8.8
if -5.087015824912348e+107 < re < 7.158781234828456e-138 or 1.1925302764880028e-107 < re < 4.3619401585985797e+95Initial program 21.4
rmApplied add-sqr-sqrt21.4
Applied pow121.4
Applied log-pow21.4
Applied times-frac21.3
rmApplied div-inv21.2
rmApplied add-sqr-sqrt21.2
Applied associate-*r*21.3
rmApplied associate-*l*21.2
if 7.158781234828456e-138 < re < 1.1925302764880028e-107Initial program 20.3
rmApplied add-sqr-sqrt20.3
Applied pow120.3
Applied log-pow20.3
Applied times-frac20.3
rmApplied div-inv20.1
rmApplied add-sqr-sqrt20.1
Applied associate-*r*20.2
Taylor expanded around 0 36.8
if 4.3619401585985797e+95 < re Initial program 50.3
rmApplied add-sqr-sqrt50.3
Applied pow150.3
Applied log-pow50.3
Applied times-frac50.3
rmApplied div-inv50.2
Taylor expanded around inf 9.0
Final simplification17.5
herbie shell --seed 2019007
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))
Time bar (total: 21.5s)Debug log
1 calls. Slowest were:
| 11.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.6b
Found 3 expressions with local error:
| 31.1b | (sqrt (+ (* re re) (* im im))) |
| 0.6b | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 0.0b | (+ (* re re) (* im im)) |
3 calls. Slowest were:
| 9.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 2.0ms | (+ (* re re) (* im im)) |
| 2.0ms | (sqrt (+ (* re re) (* im im))) |
3 calls. Slowest were:
| 147.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 27.0ms | (+ (* re re) (* im im)) |
| 26.0ms | (sqrt (+ (* re re) (* im im))) |
55 calls. Slowest were:
| 141.0ms | (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))) |
| 133.0ms | (/ (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (log 10) (log 10)) (log 10))) |
| 104.0ms | (/ (log (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (log 10)) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 7.1b
Found 4 expressions with local error:
| 31.1b | (sqrt (+ (* re re) (* im im))) |
| 0.4b | (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10))) |
| 0.2b | (* (/ 1 (sqrt (log 10))) (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10)))) |
| 0.0b | (+ (* re re) (* im im)) |
4 calls. Slowest were:
| 9.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10))) |
| 8.0ms | (* (/ 1 (sqrt (log 10))) (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10)))) |
| 2.0ms | (+ (* re re) (* im im)) |
4 calls. Slowest were:
| 199.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10))) |
| 130.0ms | (* (/ 1 (sqrt (log 10))) (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10)))) |
| 33.0ms | (sqrt (+ (* re re) (* im im))) |
| 24.0ms | (+ (* re re) (* im im)) |
138 calls. Slowest were:
| 346.0ms | (* (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))) (* (* (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10))) (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10)))) (/ (log (sqrt (+ (* re re) (* im im)))) (sqrt (log 10))))) |
| 155.0ms | (* (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))) (/ (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (sqrt (log 10)) (sqrt (log 10))) (sqrt (log 10))))) |
| 151.0ms | (/ (log (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (sqrt (log 10))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 7.0b
Found 4 expressions with local error:
| 31.1b | (sqrt (+ (* re re) (* im im))) |
| 0.2b | (* (/ 1 (sqrt (log 10))) (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10))))) |
| 0.2b | (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10)))) |
| 0.0b | (+ (* re re) (* im im)) |
4 calls. Slowest were:
| 24.0ms | (* (/ 1 (sqrt (log 10))) (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10))))) |
| 11.0ms | (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10)))) |
| 4.0ms | (+ (* re re) (* im im)) |
4 calls. Slowest were:
| 191.0ms | (* (/ 1 (sqrt (log 10))) (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10))))) |
| 176.0ms | (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10)))) |
| 33.0ms | (sqrt (+ (* re re) (* im im))) |
| 15.0ms | (+ (* re re) (* im im)) |
51 calls. Slowest were:
| 512.0ms | (* (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))) (* (* (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10)))) (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10))))) (* (log (sqrt (+ (* re re) (* im im)))) (/ 1 (sqrt (log 10)))))) |
| 439.0ms | (* (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10))))) |
| 308.0ms | (* (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))) (* (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))))) |
13 alts after pruning (13 fresh and 0 done)
Merged error: 6.9b
Found 4 expressions with local error:
| 31.1b | (sqrt (+ (* re re) (* im im))) |
| 0.3b | (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) |
| 0.2b | (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10))))) |
| 0.2b | (* (/ 1 (sqrt (log 10))) (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10)))))) |
4 calls. Slowest were:
| 33.0ms | (* (/ 1 (sqrt (log 10))) (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10)))))) |
| 22.0ms | (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10))))) |
| 8.0ms | (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) |
4 calls. Slowest were:
| 308.0ms | (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) |
| 201.0ms | (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10))))) |
| 160.0ms | (* (/ 1 (sqrt (log 10))) (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10)))))) |
| 31.0ms | (sqrt (+ (* re re) (* im im))) |
74 calls. Slowest were:
| 637.0ms | (* (* (* (/ 1 (sqrt (log 10))) (/ 1 (sqrt (log 10)))) (/ 1 (sqrt (log 10)))) (* (* (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10))))) (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10)))))) (* (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10))))))) |
| 386.0ms | (+ (log (/ 1 (sqrt (log 10)))) (+ (log (* (log (sqrt (+ (* re re) (* im im)))) (sqrt (/ 1 (sqrt (log 10)))))) (log (sqrt (/ 1 (sqrt (log 10))))))) |
| 352.0ms | (* (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (sqrt (/ 1 (sqrt (log 10)))) (sqrt (/ 1 (sqrt (log 10))))) (sqrt (/ 1 (sqrt (log 10)))))) |
14 alts after pruning (14 fresh and 0 done)
Merged error: 6.9b
57.9% (10.0b remaining)
Error of 17.5b against oracle of 7.5b and baseline of 31.2b