HP Prime complex division problem in CAS

[Albert Chan]

(11-28-2020 01:13 PM)Albert Chan Wrote:  I just noticed s itself might overflowed, if denominator is small.
We should do scaling in steps:
...

Perhaps multi-step scaling is overkill. Limiting s ≤ MAXREAL work too.

s = (MAXREAL/2) / max(0.5, |RE(x)|,|IM(x)|,|RE(y)|,|IM(y)|)
x/y = cdiv_smith(x*s, y*s)

Patching similar idea to previous cdiv Lua code:

Code:

< xr,xi,yr,yi = ldexp(xr,e),ldexp(xi,e),ldexp(yr,e),ldexp(yi,e)

> e = 2 ^ min(e,1023)
> xr,xi,yr,yi = xr*e, xi*e, yr*e, yi*e