


Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
Initial program 25.0
rmApplied add-sqr-sqrt25.0
Applied associate-/r*24.9
rmApplied div-sub25.0
rmApplied associate-/l*23.5
rmApplied *-un-lft-identity23.5
Applied times-frac21.7
Simplified21.7
Final simplification21.7
herbie shell --seed 2018362
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, imaginary part"
(/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))))
Time bar (total: 12.2s)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 22.7b
Found 3 expressions with local error:
| 22.7b | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 0.5b | (- (* x.im y.re) (* x.re y.im)) |
| 0.0b | (+ (* y.re y.re) (* y.im y.im)) |
3 calls. Slowest were:
| 13.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 4.0ms | (+ (* y.re y.re) (* y.im y.im)) |
| 3.0ms | (- (* x.im y.re) (* x.re y.im)) |
3 calls. Slowest were:
| 58.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))) |
| 38.0ms | (- (* x.im y.re) (* x.re y.im)) |
| 23.0ms | (+ (* y.re y.re) (* y.im y.im)) |
41 calls. Slowest were:
| 333.0ms | (- (pow (* x.im y.re) 3) (pow (* x.re y.im) 3)) |
| 272.0ms | (- (* (* x.im y.re) (* x.im y.re)) (* (* x.re y.im) (* x.re y.im))) |
| 150.0ms | (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im)))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 22.3b
Found 4 expressions with local error:
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 15.1b | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 0.5b | (- (* x.im y.re) (* x.re y.im)) |
4 calls. Slowest were:
| 16.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 3.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 3.0ms | (- (* x.im y.re) (* x.re y.im)) |
4 calls. Slowest were:
| 111.0ms | (/ (- (* x.im y.re) (* x.re y.im)) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 32.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 29.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 26.0ms | (- (* x.im y.re) (* x.re y.im)) |
76 calls. Slowest were:
| 347.0ms | (- (* (* x.im y.re) (* x.im y.re)) (* (* x.re y.im) (* x.re y.im))) |
| 282.0ms | (- (pow (* x.im y.re) 3) (pow (* x.re y.im) 3)) |
| 207.0ms | (/ (* (* (- (* x.im y.re) (* x.re y.im)) (- (* x.im y.re) (* x.re y.im))) (- (* x.im y.re) (* x.re y.im))) (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im))))) |
14 alts after pruning (14 fresh and 0 done)
Merged error: 20.1b
Found 4 expressions with local error:
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 9.9b | (/ (* x.im y.re) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
4 calls. Slowest were:
| 4.0ms | (/ (* x.im y.re) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
4 calls. Slowest were:
| 38.0ms | (/ (* x.im y.re) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 36.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 31.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 22.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
52 calls. Slowest were:
| 634.0ms | (/ (* (* (* x.im y.re) (* x.im y.re)) (* x.im y.re)) (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im))))) |
| 182.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
| 162.0ms | (sqrt (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im))))) |
12 alts after pruning (12 fresh and 0 done)
Merged error: 19.7b
Found 4 expressions with local error:
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 28.3b | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 9.5b | (/ (* x.re y.im) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
4 calls. Slowest were:
| 5.0ms | (/ (* x.re y.im) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 2.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
4 calls. Slowest were:
| 53.0ms | (/ (* x.re y.im) (sqrt (+ (* y.re y.re) (* y.im y.im)))) |
| 31.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 27.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
| 22.0ms | (sqrt (+ (* y.re y.re) (* y.im y.im))) |
52 calls. Slowest were:
| 678.0ms | (/ (* (* (* x.re y.im) (* x.re y.im)) (* x.re y.im)) (* (* (sqrt (+ (* y.re y.re) (* y.im y.im))) (sqrt (+ (* y.re y.re) (* y.im y.im)))) (sqrt (+ (* y.re y.re) (* y.im y.im))))) |
| 150.0ms | (sqrt (- (* (* y.re y.re) (* y.re y.re)) (* (* y.im y.im) (* y.im y.im)))) |
| 119.0ms | (sqrt (+ (* (* y.re y.re) (* y.re y.re)) (- (* (* y.im y.im) (* y.im y.im)) (* (* y.re y.re) (* y.im y.im))))) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 19.2b
0% (0.4b remaining)
Error of 21.7b against oracle of 21.3b and baseline of 21.7b