June 30, 2006
8:25 a.m.
On Fri, 30 Jun 2006, Nils Wagner wrote: [...]
Did you check general band matrices ?
No.
If you need a test for real nonsymmetric band matrices
Not this time - the wrapped routines are for symmetric or Hermetian routines only.
from scipy import linalg, diag, arange, linspace, sort def sylvester(n): # # Sylvester matrix # tmp = diag(arange(1,n),1)+diag(arange(n-1,0,-1),-1) return tmp
n = 9 A = sylvester(n) w1=linalg.eigvals(A) w1=sort(w1.real) # # Exact eigenvalues # w_ex = linspace(-(n-1),n-1,n,endpoint=True) print linalg.norm(w1-w_ex)
Nils
_______________________________________________ Scipy-dev mailing list Scipy-dev@scipy.org http://projects.scipy.org/mailman/listinfo/scipy-dev