


Bits error versus re



Bits error versus im
Results
if re < -1.98901646800801e+52Initial program 42.7
rmApplied add-cube-cbrt42.7
rmApplied add-exp-log42.7
Taylor expanded around -inf 11.2
Simplified11.2
if -1.98901646800801e+52 < re < 2.7167323166120062e+127Initial program 21.4
rmApplied pow1/221.4
Applied log-pow21.4
Applied associate-/l*21.4
rmApplied *-un-lft-identity21.4
Applied add-sqr-sqrt21.4
Applied times-frac21.5
Applied associate-/r*21.3
if 2.7167323166120062e+127 < re Initial program 55.2
rmApplied add-cube-cbrt55.2
Taylor expanded around inf 7.6
Final simplification17.3
herbie shell --seed 2019021
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))
Time bar (total: 19.5s)Debug log
| 42.0ms | 256× | body | 80 | valid |
| 8.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 28.7b
Found 3 expressions with local error:
| 27.4b | (sqrt (+ (* re re) (* im im))) |
| 0.6b | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 0.0b | (+ (* re re) (* im im)) |
| 15× | times-frac |
| 15× | add-sqr-sqrt |
| 14× | add-cube-cbrt |
| 14× | *-un-lft-identity |
| 8× | log-pow |
| 7× | pow1 |
| 5× | add-log-exp |
| 5× | associate-/l* |
| 5× | add-exp-log |
| 5× | pow1/2 |
| 5× | add-cbrt-cube |
| 4× | sqrt-div |
| 3× | flip-+ |
| 3× | sqrt-prod |
| 3× | flip3-+ |
| 3× | associate-/r* |
| 2× | div-sub |
| 2× | log-div |
| 1× | div-inv |
| 1× | div-exp |
| 1× | frac-2neg |
| 1× | sum-log |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 14.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 6.0ms | (sqrt (+ (* re re) (* im im))) |
| 2.0ms | (+ (* re re) (* im im)) |
| 137.0ms | (/ (log (sqrt (+ (* re re) (* im im)))) (log 10)) |
| 30.0ms | (sqrt (+ (* re re) (* im im))) |
| 16.0ms | (+ (* re re) (* im im)) |
| 172.0ms | (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))) |
| 154.0ms | (/ (log (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (log 10)) |
| 151.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))) |
| 70.0ms | (* -1 (/ (log (/ 1 re)) (log 10))) |
| 70.0ms | (/ (log (sqrt (+ (pow (* re re) 3) (pow (* im im) 3)))) (log 10)) |
15 alts after pruning (15 fresh and 0 done)
Merged error: 6.5b
Found 4 expressions with local error:
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 27.4b | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (cbrt (+ (* re re) (* im im))))) |
| 60× | cbrt-div |
| 30× | flip-+ |
| 30× | flip3-+ |
| 28× | frac-times |
| 26× | sqrt-div |
| 14× | associate-*l/ |
| 9× | cbrt-prod |
| 8× | associate-*r/ |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 4× | add-log-exp |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | pow1/3 |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 10.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (cbrt (+ (* re re) (* im im))))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 120.0ms | (cbrt (+ (* re re) (* im im))) |
| 99.0ms | (cbrt (+ (* re re) (* im im))) |
| 94.0ms | (cbrt (+ (* re re) (* im im))) |
| 21.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (cbrt (+ (* re re) (* im im))))) |
| 201.0ms | (sqrt (* (* (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (cbrt (+ (* re re) (* im im))))) |
| 197.0ms | (sqrt (* (* (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))) (cbrt (+ (* re re) (* im im)))) (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))))) |
| 193.0ms | (sqrt (* (* (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))))) |
| 190.0ms | (sqrt (* (* (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))))) |
| 187.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))))) |
16 alts after pruning (15 fresh and 1 done)
Merged error: 6.5b
Found 4 expressions with local error:
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 28.3b | (cbrt (+ (* re re) (* im im))) |
| 27.4b | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (exp (log (cbrt (+ (* re re) (* im im))))))) |
| 60× | cbrt-div |
| 30× | flip-+ |
| 30× | flip3-+ |
| 28× | frac-times |
| 26× | sqrt-div |
| 18× | exp-diff |
| 18× | log-div |
| 14× | associate-*l/ |
| 9× | cbrt-prod |
| 8× | associate-*r/ |
| 7× | add-cube-cbrt |
| 7× | *-un-lft-identity |
| 7× | add-sqr-sqrt |
| 4× | add-log-exp |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | pow1 |
| 3× | pow1/3 |
| 1× | sqrt-prod |
| 1× | pow1/2 |
| 17.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (exp (log (cbrt (+ (* re re) (* im im))))))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 2.0ms | (cbrt (+ (* re re) (* im im))) |
| 135.0ms | (cbrt (+ (* re re) (* im im))) |
| 115.0ms | (cbrt (+ (* re re) (* im im))) |
| 90.0ms | (cbrt (+ (* re re) (* im im))) |
| 21.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (+ (* re re) (* im im)))) (exp (log (cbrt (+ (* re re) (* im im))))))) |
| 202.0ms | (sqrt (* (* (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))) (cbrt (+ (pow (* re re) 3) (pow (* im im) 3)))) (exp (log (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))))))) |
| 191.0ms | (sqrt (* (* (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (exp (log (cbrt (+ (* re re) (* im im))))))) |
| 189.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (exp (log (cbrt (+ (* re re) (* im im))))))) |
| 187.0ms | (sqrt (* (* (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))) (cbrt (+ (pow (* re re) 3) (pow (* im im) 3)))) (exp (log (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))))))) |
| 185.0ms | (sqrt (* (* (cbrt (+ (* re re) (* im im))) (cbrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (exp (log (cbrt (+ (pow (* re re) 3) (pow (* im im) 3))))))) |
16 alts after pruning (14 fresh and 2 done)
Merged error: 6.5b
Found 4 expressions with local error:
| 28.4b | (log (+ (* re re) (* im im))) |
| 0.6b | (/ (log 10) (log (+ (* re re) (* im im)))) |
| 0.3b | (/ 1/2 (/ (log 10) (log (+ (* re re) (* im im))))) |
| 0.0b | (+ (* re re) (* im im)) |
| 24× | times-frac |
| 22× | add-cube-cbrt |
| 22× | *-un-lft-identity |
| 22× | add-sqr-sqrt |
| 20× | associate-/r* |
| 12× | pow1 |
| 8× | log-pow |
| 7× | add-exp-log |
| 6× | add-log-exp |
| 6× | add-cbrt-cube |
| 3× | associate-/l* |
| 3× | div-inv |
| 3× | log-prod |
| 2× | flip-+ |
| 2× | flip3-+ |
| 2× | frac-2neg |
| 2× | clear-num |
| 2× | log-div |
| 1× | associate-/r/ |
| 1× | rem-log-exp |
| 1× | div-exp |
| 1× | sum-log |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 6.0ms | (/ (log 10) (log (+ (* re re) (* im im)))) |
| 3.0ms | (log (+ (* re re) (* im im))) |
| 2.0ms | (+ (* re re) (* im im)) |
| 1.0ms | (/ 1/2 (/ (log 10) (log (+ (* re re) (* im im))))) |
| 132.0ms | (/ (log 10) (log (+ (* re re) (* im im)))) |
| 131.0ms | (/ 1/2 (/ (log 10) (log (+ (* re re) (* im im))))) |
| 43.0ms | (log (+ (* re re) (* im im))) |
| 15.0ms | (+ (* re re) (* im im)) |
| 134.0ms | (/ (* (* (log 10) (log 10)) (log 10)) (* (* (log (+ (* re re) (* im im))) (log (+ (* re re) (* im im)))) (log (+ (* re re) (* im im))))) |
| 111.0ms | (log (- (* (* re re) (* re re)) (* (* im im) (* im im)))) |
| 91.0ms | (* -1 (/ (log (/ -1 re)) (log 10))) |
| 83.0ms | (* -1 (/ (log (/ 1 re)) (log 10))) |
| 47.0ms | (log (+ (pow (* re re) 3) (pow (* im im) 3))) |
19 alts after pruning (18 fresh and 1 done)
Merged error: 6.5b
56.5% (10.3b remaining)
Error of 17.3b against oracle of 7.0b and baseline of 30.6b
| 1.4s | 8000× | body | 80 | valid |