Hi Karl,

Karl Rupp <r...@iue.tuwien.ac.at> writes:
> back to the original problem: The RHS is not passed correctly. Use the
> sample system attached, it is just 4x4 and should converge nicely. If
> I print the RHS vector passed to the iterative solver, it is all
> zero. Therefore, the solver doesn't even start to iterate, but instead
> returns the exact solution (vector with all zeros) directly.

How strange, but you are right! I put a print into the C++, and for some
reason the vector is 0.

> I haven't searched for the exact spot where the error occurs, but I'm
> pretty sure you will find this quickly :-)

I wish I could be so sure! Somehow, the contents of the vector get lost
somewhere during the call to the solver function. The vector is correct
immediately before the call in Python, but during the call it seems to
change. Very strange. At least I know where the bug is, though!

Thanks,


Toby


------------------------------------------------------------------------------
Managing the Performance of Cloud-Based Applications
Take advantage of what the Cloud has to offer - Avoid Common Pitfalls.
Read the Whitepaper.
http://pubads.g.doubleclick.net/gampad/clk?id=121054471&iu=/4140/ostg.clktrk
_______________________________________________
ViennaCL-devel mailing list
ViennaCL-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/viennacl-devel

Reply via email to