On Feb 24, 10:56 am, "[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
wrote:
> Support team:
>
> Here is a puzzle (to me if not to you).
>
> In sage 2.10.2 (although the problem occurs in earlier versions)
> on a CoreDuo MacIntosh running under MacOS X.5.2 we have:
>
> a = gp.sin(1)
> octave.eval('format long g')
> b = octave.sin(1)
> c = math.sin(1)
> eps = 1. - (4./3.-1.)*3.
> d = a*eps
> a
>         0.8414709848078965066525023216
> b
>         0.841470984807897
> c
>         0.8414709848078965
> (a-b)/d
>         -2.640423314436302686
> (a-c)/d
>         -553.9021661362176753
> (b-c)/d
>         -553.7921192926074518167939202
>
> Notice that a, b, c all seem to agree to about 16 digits, but when
> the 'native' sage value is mixed with the gp and octave values a
> horrendously large error occurs.

Thanks for reporting this!  This is a bug when transferring floating-
point values from Sage to interpreters such as gp and octave.  We are
now tracking this bug here: http://trac.sagemath.org/sage_trac/ticket/2294;
I expect to have a fix shortly.

Carl

--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to sage-support@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-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---

Reply via email to