Average Error: 25.0 → 21.7
Time: 15.8s
Precision: 64
Internal Precision: 128
\[\frac{x.im \cdot y.re - x.re \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}\]
\[\frac{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}} - x.re \cdot \frac{y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]

Error

Bits error versus x.re

Bits error versus x.im

Bits error versus y.re

Bits error versus y.im

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 25.0

    \[\frac{x.im \cdot y.re - x.re \cdot y.im}{y.re \cdot y.re + y.im \cdot y.im}\]
  2. Using strategy rm
  3. Applied add-sqr-sqrt25.0

    \[\leadsto \frac{x.im \cdot y.re - x.re \cdot y.im}{\color{blue}{\sqrt{y.re \cdot y.re + y.im \cdot y.im} \cdot \sqrt{y.re \cdot y.re + y.im \cdot y.im}}}\]
  4. Applied associate-/r*24.9

    \[\leadsto \color{blue}{\frac{\frac{x.im \cdot y.re - x.re \cdot y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}\]
  5. Using strategy rm
  6. Applied div-sub25.0

    \[\leadsto \frac{\color{blue}{\frac{x.im \cdot y.re}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}} - \frac{x.re \cdot y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]
  7. Using strategy rm
  8. Applied associate-/l*23.5

    \[\leadsto \frac{\color{blue}{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}}} - \frac{x.re \cdot y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]
  9. Using strategy rm
  10. Applied *-un-lft-identity23.5

    \[\leadsto \frac{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}} - \frac{x.re \cdot y.im}{\color{blue}{1 \cdot \sqrt{y.re \cdot y.re + y.im \cdot y.im}}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]
  11. Applied times-frac21.7

    \[\leadsto \frac{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}} - \color{blue}{\frac{x.re}{1} \cdot \frac{y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]
  12. Simplified21.7

    \[\leadsto \frac{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}} - \color{blue}{x.re} \cdot \frac{y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]
  13. Final simplification21.7

    \[\leadsto \frac{\frac{x.im}{\frac{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}{y.re}} - x.re \cdot \frac{y.im}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}}{\sqrt{y.re \cdot y.re + y.im \cdot y.im}}\]

Reproduce

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))))

Details

Time bar (total: 12.2s)Debug log

start78.0ms

Algorithm
intervals

setup55.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 22.7b

localize37.0ms

Local error

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))

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
3 → 55
Calls

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))

series119.0ms

Counts
3 → 9
Calls

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))

simplify1.3s

Counts
41 → 64
Calls

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))))

prune786.0ms

Pruning

9 alts after pruning (9 fresh and 0 done)

Merged error: 22.3b

localize16.0ms

Local error

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))

rewrite28.0ms

Algorithm
rewrite-expression-head
Counts
4 → 84
Calls

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))

series198.0ms

Counts
4 → 12
Calls

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))

simplify2.2s

Counts
76 → 96
Calls

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)))))

prune1.2s

Pruning

14 alts after pruning (14 fresh and 0 done)

Merged error: 20.1b

localize11.0ms

Local error

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))))

rewrite10.0ms

Algorithm
rewrite-expression-head
Counts
4 → 69
Calls

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)))

series127.0ms

Counts
4 → 12
Calls

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)))

simplify1.5s

Counts
52 → 81
Calls

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)))))

prune1.1s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 19.7b

localize20.0ms

Local error

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))))

rewrite11.0ms

Algorithm
rewrite-expression-head
Counts
4 → 69
Calls

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)))

series132.0ms

Counts
4 → 12
Calls

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)))

simplify1.5s

Counts
52 → 81
Calls

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)))))

prune1.1s

Pruning

11 alts after pruning (11 fresh and 0 done)

Merged error: 19.2b

regimes523.0ms

Accuracy

0% (0.4b remaining)

Error of 21.7b against oracle of 21.3b and baseline of 21.7b

bsearch3.0ms