[SciPy-user] conjugate gradients solver - operator adjoint
Lubos Vrbka
lists at vrbka.net
Wed Apr 22 04:00:07 EDT 2009
hi,
> The conjugate gradient method assumes that the matrix is hermitian, hence
> no need for specifying adjoint separately.
>
> The algorithm you refer to is probably something different. A reference
> to the paper would come handy.
the paper is journal of computational physics 61 (1985) 280-285. in case
of need/interest i could forward it to you. for obvious reasons i don't
attach it to this message...
would be glad to know your opinion on this problem... will check also
the other options. i think the problem is, that in order for the
algorithm to work (it's newton-raphson + conjugate gradients), it's
necessary to have the exact adjoint in the solution...
best,
--
Lubos _ at _"
http://www.lubos.vrbka.net
More information about the SciPy-User
mailing list