#6115: make symbolic matrices use pynac symbolics
----------------------------+-----------------------------------------------
Reporter: jason | Owner: was
Type: enhancement | Status: new
Priority: major | Milestone: sage-4.0.1
Component: linear algebra | Keywords:
----------------------------+-----------------------------------------------
Comment(by jason):
Almost positive review. Here is a problem:
{{{
sage: a=matrix(SR, [[1,2],[3,4]])
sage: a.eig
a.eigenmatrix_left a.eigenspaces a.eigenspaces_right
a.eigenvectors_left
a.eigenmatrix_right a.eigenspaces_left a.eigenvalues
a.eigenvectors_right
sage: a.eigenm
a.eigenmatrix_left a.eigenmatrix_right
sage: a.eigenmatrix_right()
/home/grout/.sage/temp/good/26004/_home_grout__sage_init_sage_0.py:1:
UserWarning: Using generic algorithm for an inexact ring, which may result
in garbarge from numerical precision issues.
# -*- coding: utf-8 -*-
/home/grout/.sage/temp/good/26004/_home_grout__sage_init_sage_0.py:1:
UserWarning: Using generic algorithm for an inexact ring, which will
probably give incorrect results due to numerical precision issues.
# -*- coding: utf-8 -*-
---------------------------------------------------------------------------
TypeError Traceback (most recent call
last)
/home/grout/.sage/temp/good/26004/_home_grout__sage_init_sage_0.py in
<module>()
/home/grout/download/sage-sage-4.0.alpha0.5/local/lib/python2.5/site-
packages/sage/matrix/matrix2.so in
sage.matrix.matrix2.Matrix.eigenmatrix_right
(sage/matrix/matrix2.c:21183)()
3505 True
3506 """
-> 3507 D,P=self.transpose().eigenmatrix_left()
3508 return D,P.transpose()
3509
/home/grout/download/sage-sage-4.0.alpha0.5/local/lib/python2.5/site-
packages/sage/matrix/matrix2.so in
sage.matrix.matrix2.Matrix.eigenmatrix_left
(sage/matrix/matrix2.c:20868)()
3440 """
3441 from sage.misc.flatten import flatten
-> 3442 evecs = self.eigenvectors_left()
3443 D =
sage.matrix.constructor.diagonal_matrix(flatten([[e[0]]*e[2] for e in
evecs]))
3444 rows = []
/home/grout/download/sage-sage-4.0.alpha0.5/local/lib/python2.5/site-
packages/sage/matrix/matrix2.so in
sage.matrix.matrix2.Matrix.eigenvectors_left
(sage/matrix/matrix2.c:20199)()
3324 from sage.rings.qqbar import QQbar
3325 from sage.categories.homset import hom
-> 3326 eigenspaces =
self.eigenspaces_left(algebraic_multiplicity=True)
3327 evec_list=[]
3328 n = self._nrows
/home/grout/download/sage-sage-4.0.alpha0.5/local/lib/python2.5/site-
packages/sage/matrix/matrix2.so in
sage.matrix.matrix2.Matrix.eigenspaces_left
(sage/matrix/matrix2.c:18811)()
3061 i = 0
3062 for h, e in G:
-> 3063 if h.degree() == 1:
3064 alpha = -h[0]/h[1]
3065 F = alpha.parent()
TypeError: degree() takes exactly one argument (0 given)
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/6115#comment:2>
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
-~----------~----~----~----~------~----~------~--~---
