I want to use "RDF" linear algebra. As far as i understand, operations are implemented using numpy/scipy. But many things differ from a direct call to scipy modules; for example:
(computing the Singular Value Decomposition of a matrix):
sage: A=matrix(RDF,[[1,3,2],[1,2,3],[0,5,2],[1, 1, 1]])
sage: U,Sig,V=A.SVD()
but, a direct call to scipy (even in sage) is:
from scipy import *
from scipy.linalg import *
A=mat('[1 3 2; 1 2 3; 0 5 2; 1, 1, 1]')
M,N=A.shape
U,s,Vh=svd(A)
-Qusetion: is there some "dictionary", some documentation about how
the scipy functions are "mapped" to sage? and what can be directly used?
One can use A.LU(), A.QR()... witeh A a matrix(RDF,...), all this
without importing explicitly anything.
-Remark: In the preceding examples, if you compare as given by U,Sig,V=A.SVD() and U given by U,s,Vh=svd(A) they are transposed :-( Yours, sincerely t. -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org To unsubscribe, reply using "remove me" as the subject.
<<attachment: tdumont.vcf>>
