#9689: Numerical noise on devel/sage-main/sage/symbolic/expression.pyx computing
sinh(1)
------------------------+---------------------------------------------------
Reporter: drkirkby | Owner: mvngu
Type: defect | Status: needs_review
Priority: major | Milestone: sage-4.5.3
Component: doctest | Keywords:
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
------------------------+---------------------------------------------------
Changes (by drkirkby):
* status: new => needs_review
Old description:
> Whilst there is no complete 64-bit build of Sage on Solaris x86, a
> sufficiently large part of Sage does build (with a few changes) on
> Solaris 10 x86. When built on 'fulvia', a Dell Optiplex with Xeon
> processors, there was a numerical noise issue - see #9099
>
> {{{
> sage -t -long devel/sage/sage/symbolic/expression.pyx
> **********************************************************************
> File "/home/palmieri/fulvia/sage-4.5.2.rc0/devel/sage-
> main/sage/symbolic/expression.pyx", line 498\
> 3:
> sage: maxima('sinh(1.0)')
> Expected:
> 1.175201193643801
> Got:
> 1.175201193643802
> }}}
>
> A computation with Mathematica, using 60 digits of precision gives
>
> {{{
> In[2]:= N[Sinh[1],60]
>
> Out[2]= 1.17520119364380145688238185059560081515571798133409587022957
> }}}
>
> The absolute error on Solaris x86 is slighly higher than seen on some
> other systems, but this is still a perfectly acceptable result.
>
> This should be fairly easy to fix. I'll make a patch later today
>
> Dave
New description:
Whilst there is no complete 64-bit build of Sage on Solaris x86, a
sufficiently large part of Sage does build (with a few changes) on Solaris
10 x86. When built on 'fulvia', a Dell Optiplex with Xeon processors,
there was a numerical noise issue - see #9099
{{{
sage -t -long devel/sage/sage/symbolic/expression.pyx
**********************************************************************
File "/home/palmieri/fulvia/sage-4.5.2.rc0/devel/sage-
main/sage/symbolic/expression.pyx", line 498\
3:
sage: maxima('sinh(1.0)')
Expected:
1.175201193643801
Got:
1.175201193643802
}}}
A computation with Mathematica, using 60 digits of precision gives
{{{
In[2]:= N[Sinh[1],60]
Out[2]= 1.17520119364380145688238185059560081515571798133409587022957
}}}
The absolute error on Solaris x86 is slightly higher than seen on some
other systems, but this is still a perfectly acceptable result.
This should be fairly easy to fix. I'll make a patch later today
Dave
--
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/9689#comment:2>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
