>
> If I read this correctly, gamma(1+1/2*%i),numer won't evaluate
> numerically, because (numberp ji) is false. My patch does. But with
> my patch, gamma(1+0.5*%i) does not give a number unless you specify
> numer. So both are deficient in some way.
>
Standard Maxima semantics are that operations on exact numbers are exact,
except with numer=true, so gamma(1+1/2*%i) should not return a float. Floats
and bfloats are inexact.
Currently, Maxima does not have the notion of an inexact complex in the same
way, so though sin(1.0) returns a float, sin(1.0+2.0*%i) does not.
It seems reasonable to define inexact+%i*exact and exact+%i*inexact as
inexact, so that sin(1+2.0*%i), gamma(1+0.5*%i), and gamma(1.0+%i) return
floats.
-s