[SciPy-User] scipy.sparse.linalg.bicg for Non Symmetric Matrix.
Thu Dec 9 19:42:28 CST 2010
Am using the iterative methods of scipy.sparse.linalg for solving a linear
system of equations Ax = b. My matrix A is non symmetric. I've been using
the scipy.sparse.linalg.cg() function multiplying both matrix "A" and "b"
with the transpose of A so the matrix will become symmetric:
Asym = matrix(dot(A.T,A))
bsym = matrix(dot(A.T,b))
sol = cg(A,b,tol = 1e-10,maxiter=30)
Also i've been reading about the biconjugate gradient method, and if i am
not mistaken the literature says that this method works on
non-symmetric matrix, but when i try to use scipy.sparse.linalg.bicg() it
sol = bicg(A,b,tol = 1e-10,maxiter=30)
line 74, in bicg
A,M,x,b,postprocess = make_system(A,M,x0,b,xtype)
line 65, in make_system
raise ValueError('expected square matrix (shape=%s)' % shape)
NameError: global name 'shape' is not defined
Any help will be greatly appreciated, am comparing this methods for my data
with an important emphasis on speed for my master thesis, so any discrepancy
with the theory would be a great problem for me.
Thanks in advance,
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User