[SciPy-User] iterative matrix methods seem slow

[Pauli Virtanen]

On Sat, 25 Jun 2011 13:04:46 +0300, eat wrote:
[clip]
> Anyway it will produce similar looking results than original

You can also look at the norm of the residual to check the solution,
no need to rely on visual inspection. The error goes down to
machine precision.

> C:\Python27\lib\site-packages\pyamg\util\linalg.py:233: ComplexWarning:
> Casting
> complex values to real discards the imaginary part
>  H[i,j] = numpy.dot(numpy.conjugate(numpy.ravel(v)), numpy.ravel(w))

No idea. PyAMG's a black box to me.