#10042: Doctest failure in sage/rings/polynomial/polynomial_element.pyx
-----------------------+----------------------------------------------------
Reporter: mpatel | Owner: mvngu
Type: defect | Status: needs_review
Priority: blocker | Milestone: sage-4.6
Component: doctest | Keywords:
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
-----------------------+----------------------------------------------------
Comment(by zimmerma):
> thanks Dima for giving the instructions on how to test DGEEV. I'll give
it a try.
on a 64-bit Core 2 I get:
{{{
tarte% ./a.out < inp.txt
DGEEV Example Program Results
Eigenvalue( 1) = 1.77245385090551588191940000000000000E+00
Eigenvector( 1)
8.7095E-01
4.9138E-01
Eigenvalue( 2) = -1.77245385090551588191940000000000000E+00
Eigenvector( 2)
-8.7095E-01
4.9138E-01
}}}
On Dave machine I get:
{{{
zimme...@hawk:~$ ./a.out < inp.txt
DGEEV Example Program Results
Eigenvalue( 1) = 1.77245385090551588191940000000000000E+00
Eigenvector( 1)
8.7095E-01
4.9138E-01
Eigenvalue( 2) = -1.77245385090551610396400000000000000E+00
Eigenvector( 2)
-8.7095E-01
4.9138E-01
}}}
The second eigenvalue differs. Now we have to investigate inside the DGEEV
call.
If it performs only additions, subtractions, multiplications, divisions
and square
roots, there should be no difference according to IEEE 754 (unless
different formats are used).
Paul
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10042#comment:71>
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.
