On Wed, Jun 11, 2008 at 1:08 PM, Ondrej Certik <[EMAIL PROTECTED]> wrote: > On Wed, Jun 11, 2008 at 11:32 AM, Ondrej Certik <[EMAIL PROTECTED]> wrote: >> On Tue, Jun 10, 2008 at 11:04 PM, mabshoff <[EMAIL PROTECTED]> wrote: >>> >>> Riccardo Gori wrote: >>>> Hello, >>>> >>> >>> Hi Riccardo, >>> >>> thanks for the bug report. I have forwarded it to sage-devel. >>> >>> Cheers, >>> >>> Michael >>>> In SAGE-3.0.2 with a Intel Mac OSX 10.5 I found the following bug: >>>> >>>> Trying to evaluate a sympy function gives a wrong result. With sympy in a >>>> python shell everything works well. >>>> >>>> To reproduce: >>>> >>>> sage: import sympy >>>> sage: sympy.sin(1.1) >>>> sin(1) >>>> sage: sympy.sin(float(1.2)) >>>> sin(1) >>>> sage: sympy.sin(1.3) >>>> sin(1.300000000000000044408920985) >>>> sage: sympy.sin(float(1.4)) >>>> sin(1.39999999999999991118215803) >>>> >>>> I didn't find a way to get a correct result even with float(). >> >> I created an issue for that: >> >> http://code.google.com/p/sympy/issues/detail?id=881 >> >> It's a serious bug. However, so far I cannot reproduce it nor in >> sympy, nor in sage: >> >> sage: import sympy >> sage: sympy.sin(float(1.1)) >> sin(1.10000000000000008881784197) >> sage: sympy.sin(1.1) >> sin(1.10000000000000008881784197) >> >> I'll dig further and fix this in the next release of sympy. > > Ok, I reproduced it, one needs to just start sage and do: > > sage: import sympy > sage: sympy.sin(1.1) > sin(1) > > if you start sage and do: > > sage: import sympy > sage: sympy.sin(float(1.1)) > sin(1.10000000000000008881784197) > sage: sympy.sin(1.1) > sin(1.10000000000000008881784197) > > it works. I'll look into it and fix it.
Ok, the problem was in sympify(). I fixed that, wrote tests and sent it for a review here: http://code.google.com/p/sympy/issues/detail?id=881 it will end up in the next sympy release. Ondrej --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to [EMAIL PROTECTED] For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---