[SciPy-user] conjugate gradients solver - operator adjoint
Wed Apr 22 03:00:48 CDT 2009
Wed, 22 Apr 2009 07:46:49 +0000, Pauli Virtanen kirjoitti:
>> (which should be, if i understand the method correctly, necessary for
>> getting the result, since it's impossible to invert something that big
> This is probably debatable. There are many iterative algorithms that can
> invert non-symmetric matrices not needing to know the adjoint.
Oops, I'm mistaken here, scipy.sparse.linalg.bicg and .qmr require the
adjoint. (However, it is automatically computed if you give in a sparse
So the algorithm you described probably is scipy.sparse.linalg.bicg or
some related method.
More information about the SciPy-user