2011/6/6 Dürrwang, Jürgen <span dir="ltr"><<a href="mailto:Juergen.Duerrwang@iosb.fraunhofer.de">Juergen.Duerrwang@iosb.fraunhofer.de</a>></span><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div lang="DE" link="blue" vlink="purple"><div><p class="MsoNormal">Hi,</p><p class="MsoNormal"> </p><p class="MsoNormal"><span lang="EN-US">I am working on a NVIDIA CUDA improved cg solver, with fast preconditioning. Since I have tested PETSC PCJACOBI with ILU(0) on each block I have some nice results. Unfortunately the pc is on CPU, so I try to port them on GPU. Now my question is, has body some material about this PETSC implementation?</span></p>
</div></div></blockquote><div><br></div><div>The ILU(0) algorithm is described in Yousef Saad's book on Iterative Methods. It is notoriously hard to parallelize.</div><div>ASM and Block-Jacobi jsut decompose the system into blocks. They are also described in the book.</div>
<div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div lang="DE" link="blue" vlink="purple"><div><p class="MsoNormal"> </p>
<p class="MsoNormal">Kind regards,</p><p class="MsoNormal">Jürgen Dürrwang</p></div></div></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>