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 -~----------~----~----~----~------~----~------~--~---