[GiNaC-list] simple patch for Li_evalf
Jens Vollinga
vollinga at thep.physik.uni-mainz.de
Thu Sep 23 16:09:03 CEST 2004
Hi,
On Thu, Sep 23, 2004 at 11:27:09AM +0400, Sheplyakov Alexei wrote:
> gives
>
> Li(5,exp((1/3*I)*Pi)) = Li(5,exp((1/3*I)*Pi))
>
> instead of expected
>
> Li(5,exp((1/3*I)*Pi)) = 0.48005914589970829926+0.89203268884598503196*I
>
> Attached patch seems to solve the problem.
GiNaC functions do only return floating point numbers, if their arguments
already contain floating point numbers. So actually in the example above
everything is as expected. If you want a floating point result, you
either have to call evalf() explicitly or write something like
...(1/3.0*I)*...
Regards,
Jens
More information about the GiNaC-list
mailing list