I am solving 2-d euler equations using an unstructured grid fvm code. I want to use GMRES for convergence acceleration. This needs a good preconditioner. My tests with a block-diagonal preconditioner are not very encouraging. This has been noted in the literature already. I want to try ILU(0) but I am not able to get any good reference on how to implement this. I am looking for some paper/textbook that explains how to implement ILU(0) for a sparse block matrix. Also are there any libraries in fortran available for doing this ?

Dear Praveen,

There is some information available in the following doctoral thesis.

Phillipe Geuzaine, "An Implicit Upwind Finite Volume method for compressible turbulent flows on unstrucutred meshes" , Univ. of Liege, April 1999.

Hope this helps.

Regards,

Ganesh

visit Y. Saad homepage: http://www-users.cs.umn.edu/~saad/

first version of his book is freely avalible contains your needed materials.

his codes e.g. SPARSEKIT have implementation of ILU type preconditioner e.g. ILU(0)

there are several free opensources that contains ILU(0) precond, it is sufficient to look for opensource iterative solvers.

Thanks to Ganesh and RT. The SPARSEKIT package seems to be ideal for me; it has a lot of preconditioners.