



Bits error versus x
Results
| Original | 61.0 |
|---|---|
| Target | 0.3 |
| Herbie | 0.0 |
Initial program 61.0
Simplified60.1
rmApplied sub-neg60.1
Applied log1p-def0.0
rmApplied add-log-exp0.0
Final simplification0.0
herbie shell --seed 2019010 +o rules:numerics
(FPCore (x)
:name "qlog (example 3.10)"
:pre (and (< -1 x) (< x 1))
:herbie-target
(- (+ (+ (+ 1 x) (/ (* x x) 2)) (* 5/12 (pow x 3))))
(/ (log (- 1 x)) (log (+ 1 x))))
Time bar (total: 18.1s)Debug log
| 8.0ms | (/ (log (- 1 x)) (log (+ 1 x))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 60.3b
Found 2 expressions with local error:
| 5.2b | (log (- 1 x)) |
| 0.0b | (/ (log (- 1 x)) (log1p x)) |
| 12× | add-cube-cbrt |
| 12× | times-frac |
| 12× | *-un-lft-identity |
| 12× | add-sqr-sqrt |
| 7× | pow1 |
| 5× | log-pow |
| 5× | add-exp-log |
| 4× | associate-/l* |
| 4× | add-cbrt-cube |
| 4× | log-div |
| 3× | associate-/r* |
| 3× | log-prod |
| 2× | add-log-exp |
| 2× | flip-- |
| 2× | log1p-expm1-u |
| 2× | div-sub |
| 2× | flip3-- |
| 2× | log1p-def |
| 2× | expm1-log1p-u |
| 1× | div-inv |
| 1× | rem-log-exp |
| 1× | div-exp |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 6.0ms | (/ (log (- 1 x)) (log1p x)) |
| 1.0ms | (log (- 1 x)) |
| 100.0ms | (/ (log (- 1 x)) (log1p x)) |
| 38.0ms | (log (- 1 x)) |
| 339.0ms | (- (+ x (+ (* 1/3 (pow x 3)) (* 1/2 (pow x 2))))) |
| 291.0ms | (- (log -1) (+ (log (/ 1 x)) (+ (* 1/2 (/ 1 (pow x 2))) (/ 1 x)))) |
| 268.0ms | (- (+ (log (/ -1 x)) (+ (* 1/2 (/ 1 (pow x 2))) (/ 1 x)))) |
| 204.0ms | (- (+ (* 2 (/ 1 (* (pow (log (/ 1 x)) 2) (pow x 2)))) (+ (* 2 (/ 1 (* (log (/ 1 x)) x))) (+ (* 1/2 (/ (log -1) (* (pow (log (/ 1 x)) 2) (pow x 2)))) 1))) (+ (/ (log -1) (* (pow (log (/ 1 x)) 3) (pow x 2))) (+ (/ (log -1) (log (/ 1 x))) (/ (log -1) (* (pow (log (/ 1 x)) 2) x))))) |
| 164.0ms | (- (+ (/ (log (/ -1 x)) (* x (pow (- (log -1) (log (/ -1 x))) 2))) (/ 1 (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (+ (/ (log (/ -1 x)) (- (log -1) (log (/ -1 x)))) (+ (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 3))) (+ (/ 1 (* x (- (log -1) (log (/ -1 x))))) (+ (* 1/2 (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (* 1/2 (/ 1 (* (pow x 2) (- (log -1) (log (/ -1 x))))))))))) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.0b
Found 1 expressions with local error:
| 0.0b | (/ (log1p (- x)) (log1p x)) |
| 9× | add-cube-cbrt |
| 9× | times-frac |
| 9× | *-un-lft-identity |
| 9× | add-sqr-sqrt |
| 3× | associate-/l* |
| 3× | add-exp-log |
| 3× | add-cbrt-cube |
| 3× | associate-/r* |
| 1× | add-log-exp |
| 1× | div-inv |
| 1× | log1p-expm1-u |
| 1× | div-exp |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 1× | pow1 |
| 1× | expm1-log1p-u |
| 7.0ms | (/ (log1p (- x)) (log1p x)) |
| 86.0ms | (/ (log1p (- x)) (log1p x)) |
| 211.0ms | (- (+ (/ (log (/ -1 x)) (* x (pow (- (log -1) (log (/ -1 x))) 2))) (/ 1 (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (+ (/ (log (/ -1 x)) (- (log -1) (log (/ -1 x)))) (+ (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 3))) (+ (/ 1 (* x (- (log -1) (log (/ -1 x))))) (+ (* 1/2 (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (* 1/2 (/ 1 (* (pow x 2) (- (log -1) (log (/ -1 x))))))))))) |
| 164.0ms | (- (+ (* 2 (/ 1 (* (pow (log (/ 1 x)) 2) (pow x 2)))) (+ (* 2 (/ 1 (* (log (/ 1 x)) x))) (+ (* 1/2 (/ (log -1) (* (pow (log (/ 1 x)) 2) (pow x 2)))) 1))) (+ (/ (log -1) (* (pow (log (/ 1 x)) 3) (pow x 2))) (+ (/ (log -1) (log (/ 1 x))) (/ (log -1) (* (pow (log (/ 1 x)) 2) x))))) |
| 159.0ms | (- (+ x (+ (* 1/2 (pow x 2)) 1))) |
| 139.0ms | (/ (* (* (log1p (- x)) (log1p (- x))) (log1p (- x))) (* (* (log1p x) (log1p x)) (log1p x))) |
| 16.0ms | (/ (* (cbrt (log1p (- x))) (cbrt (log1p (- x)))) 1) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0b
Found 2 expressions with local error:
| 0.0b | (/ (log1p (- x)) (log1p x)) |
| 0.0b | (exp (/ (log1p (- x)) (log1p x))) |
| 18× | times-frac |
| 17× | add-cube-cbrt |
| 17× | *-un-lft-identity |
| 17× | add-sqr-sqrt |
| 13× | exp-prod |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 3× | add-log-exp |
| 3× | associate-/l* |
| 3× | associate-/r* |
| 2× | div-inv |
| 2× | log1p-expm1-u |
| 2× | pow1 |
| 2× | expm1-log1p-u |
| 1× | rem-exp-log |
| 1× | div-exp |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | cbrt-undiv |
| 7.0ms | (/ (log1p (- x)) (log1p x)) |
| 6.0ms | (exp (/ (log1p (- x)) (log1p x))) |
| 111.0ms | (exp (/ (log1p (- x)) (log1p x))) |
| 80.0ms | (/ (log1p (- x)) (log1p x)) |
| 424.0ms | (- (+ (* 1/2 (/ (* (log -1) (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x)))))) (* (pow (log (/ 1 x)) 2) (pow x 2)))) (+ (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x))))) (+ (* 4 (/ (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x))))) (* (pow (log (/ 1 x)) 2) (pow x 2)))) (+ (* 2 (/ (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x))))) (* (log (/ 1 x)) x))) (* 1/2 (/ (* (pow (log -1) 2) (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x)))))) (* (pow (log (/ 1 x)) 4) (pow x 2)))))))) (+ (* 3 (/ (* (log -1) (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x)))))) (* (pow (log (/ 1 x)) 3) (pow x 2)))) (/ (* (log -1) (exp (* -1 (/ (- (log -1) (log (/ 1 x))) (log (/ 1 x)))))) (* (pow (log (/ 1 x)) 2) x)))) |
| 287.0ms | (- (exp -1) (+ (* 1/12 (* (pow x 3) (exp -1))) (* x (exp -1)))) |
| 238.0ms | (- (+ (* 2 (/ 1 (* (pow (log (/ 1 x)) 2) (pow x 2)))) (+ (* 2 (/ 1 (* (log (/ 1 x)) x))) (+ (* 1/2 (/ (log -1) (* (pow (log (/ 1 x)) 2) (pow x 2)))) 1))) (+ (/ (log -1) (* (pow (log (/ 1 x)) 3) (pow x 2))) (+ (/ (log -1) (log (/ 1 x))) (/ (log -1) (* (pow (log (/ 1 x)) 2) x))))) |
| 182.0ms | (- (+ (/ (log (/ -1 x)) (* x (pow (- (log -1) (log (/ -1 x))) 2))) (/ 1 (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (+ (/ (log (/ -1 x)) (- (log -1) (log (/ -1 x)))) (+ (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 3))) (+ (/ 1 (* x (- (log -1) (log (/ -1 x))))) (+ (* 1/2 (/ (log (/ -1 x)) (* (pow x 2) (pow (- (log -1) (log (/ -1 x))) 2)))) (* 1/2 (/ 1 (* (pow x 2) (- (log -1) (log (/ -1 x))))))))))) |
| 153.0ms | (- (+ x (+ (* 1/2 (pow x 2)) 1))) |
2 alts after pruning (1 fresh and 1 done)
Merged error: 0b
Found 3 expressions with local error:
| 0.3b | (* 1/12 x) |
| 0.0b | (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1)) |
| 0.0b | (log (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1))) |
| 4× | add-cube-cbrt |
| 4× | add-exp-log |
| 4× | *-un-lft-identity |
| 4× | pow1 |
| 4× | add-sqr-sqrt |
| 3× | add-log-exp |
| 3× | log1p-expm1-u |
| 3× | add-cbrt-cube |
| 3× | log-prod |
| 3× | expm1-log1p-u |
| 1× | fma-udef |
| 1× | *-commutative |
| 1× | log-pow |
| 1× | rem-log-exp |
| 1.0ms | (log (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1))) |
| 0.0ms | (* 1/12 x) |
| 0.0ms | (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1)) |
| 125.0ms | (log (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1))) |
| 76.0ms | (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1)) |
| 7.0ms | (* 1/12 x) |
| 345.0ms | (log (* (cbrt (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1))) (cbrt (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1))))) |
| 284.0ms | (log (sqrt (fma (fma (* 1/12 x) (* x x) x) (- (exp -1)) (exp -1)))) |
| 280.0ms | (- (exp -1) (+ (* 1/12 (* (pow x 3) (exp -1))) (* x (exp -1)))) |
| 270.0ms | (- (exp -1) (+ (* 1/12 (* (pow x 3) (exp -1))) (* x (exp -1)))) |
| 232.0ms | (- (+ (log (* 1/12 (exp -1))) (* 12 (/ 1 (pow x 2)))) (+ (* 3 (log (/ -1 x))) (* 12 (/ 1 (pow x 3))))) |
2 alts after pruning (0 fresh and 2 done)
Merged error: 0b
0% (0.0b remaining)
Error of 0.0b against oracle of 0.0b and baseline of 0.0b