float and bfloat issues



These were sent to me by Jim FitzSimons, who wrote the bigfloat
implementation of elliptic integrals that I translated into Lisp.

bfloat(sqrt(%i)) -> 1.0b0*(-1)^(1/4)

Should this be more like .707+.707*%i?

Another issue is

float(polarform(sqrt(1+%i))) -> 1.18*2.718^(%i*%pi/8)

But

bfloat(polarform(sqrt(1+%i))) -> 1.18*(3.82b-1*%i+9.23b-1)

Not sure which is the correct answer.  In some ways, the float result is
better since it preserves the polarform, but doesn't float the phase
part.  The bfloat result is better since everything is a bfloat, but the
polarform is converted to rectangular.

Ray