#4502: [with patch; needs work] numerical noise in matrix_double_dense on intel
mac os X 10.5: inverting a singular matrix
-------------------------------------+--------------------------------------
Reporter: jhpalmieri | Owner: somebody
Type: defect | Status: new
Priority: major | Milestone: sage-3.2
Component: linear algebra | Resolution:
Keywords: numerical noise, matrix |
-------------------------------------+--------------------------------------
Comment (by mabshoff):
With
{{{
diff -r 50c6651a1286 sage/matrix/matrix_double_dense.pyx
--- a/sage/matrix/matrix_double_dense.pyx Tue Nov 18 08:58:27 2008
-0800
+++ b/sage/matrix/matrix_double_dense.pyx Tue Nov 18 13:49:26 2008
-0800
@@ -441,7 +441,7 @@
sage: A.determinant() < 10e-12
True
- sage: ~A
+ sage: ~A # random
Traceback (most recent call last):
...
LinAlgError: singular matrix
}}}
I get
{{{
sage -t devel/sage/sage/matrix/matrix_double_dense.pyx
**********************************************************************
File "/scratch/mabshoff/release-
cycle/sage-3.2.rc2/devel/sage/sage/matrix/matrix_double_dense.pyx", line
444:
sage: print "ignore this"; ~A # random
Exception raised:
Traceback (most recent call last):
File "/scratch/mabshoff/release-
cycle/sage-3.2.rc2/local/bin/ncadoctest.py", line 1231, in run_one_test
self.run_one_example(test, example, filename, compileflags)
File "/scratch/mabshoff/release-
cycle/sage-3.2.rc2/local/bin/sagedoctest.py", line 38, in run_one_example
OrigDocTestRunner.run_one_example(self, test, example, filename,
compileflags)
File "/scratch/mabshoff/release-
cycle/sage-3.2.rc2/local/bin/ncadoctest.py", line 1172, in run_one_example
compileflags, 1) in test.globs
File "<doctest __main__.example_15[10]>", line 1, in <module>
print "ignore this"; ~A # random###line 444:
sage: print "ignore this"; ~A # random
File "matrix_double_dense.pyx", line 460, in
sage.matrix.matrix_double_dense.Matrix_double_dense.__invert__
(sage/matrix/matrix_double_dense.c:3878)
File "/scratch/mabshoff/release-
cycle/sage-3.2.rc2/local/lib/python2.5/site-
packages/scipy/linalg/basic.py", line 306, in inv
if info>0: raise LinAlgError, "singular matrix"
LinAlgError: singular matrix
**********************************************************************
1 items had failures:
1 of 11 in __main__.example_15
***Test Failed*** 1 failures.
For whitespace errors, see the file /scratch/mabshoff/release-
cycle/sage-3.2.rc2/tmp/.doctest_matrix_double_dense.py
[2.4 s]
exit code: 1024
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4502#comment:12>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
