



Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
| Original | 26.1 |
|---|---|
| Target | 0.5 |
| Herbie | 13.4 |
if (* (sqrt (/ 1 (hypot c d))) (* (sqrt (/ 1 (hypot c d))) (/ (- (* c b) (* a d)) (hypot c d)))) < -1.7761824129243311e+308 or 1.7795226370784402e+308 < (* (sqrt (/ 1 (hypot c d))) (* (sqrt (/ 1 (hypot c d))) (/ (- (* c b) (* a d)) (hypot c d)))) Initial program 62.5
rmApplied add-sqr-sqrt62.5
Applied *-un-lft-identity62.5
Applied times-frac62.5
Applied simplify62.5
Applied simplify62.0
rmApplied associate-*r/62.0
Applied simplify62.0
Taylor expanded around 0 62.8
Applied simplify47.7
if -1.7761824129243311e+308 < (* (sqrt (/ 1 (hypot c d))) (* (sqrt (/ 1 (hypot c d))) (/ (- (* c b) (* a d)) (hypot c d)))) < 1.7795226370784402e+308Initial program 13.5
rmApplied add-sqr-sqrt13.5
Applied *-un-lft-identity13.5
Applied times-frac13.5
Applied simplify13.5
Applied simplify1.6
rmApplied associate-*r/1.6
Applied simplify1.5
Time bar (total: 1.1m)Debug logProfile
herbie shell --seed 2018178 +o rules:numerics
(FPCore (a b c d)
:name "Complex division, imag part"
:herbie-target
(if (< (fabs d) (fabs c)) (/ (- b (* a (/ d c))) (+ c (* d (/ d c)))) (/ (+ (- a) (* b (/ c d))) (+ d (* c (/ c d)))))
(/ (- (* b c) (* a d)) (+ (* c c) (* d d))))