<div dir="ltr"><div><br></div><div>At a practical level you may be looking for parallel ways to use PETSc's native (but serial) ILU.  The simplest parallel extension is processor-block Jacobi which applies ILU on the processor-owned degrees of freedom but with no overlap:</div><div><br></div><div>-pc_type bjacobi -sub_pc_type ilu</div><div><br></div><div>The overlapping version is with ASM:</div><div><br></div><div>-pc_type asm -sub_pc_type ilu</div><div><br></div><div>where -pc_asm_overlap controls the overlap distance.  Finally one can use the serial ILU as a smoother in geometric multigrid, for example in combination with overlapping ASM:</div><div><br></div><div>-pc_type mg -mg_levels_pc_type asm -mg_levels_sub_pc_type ilu<br></div><div><br></div>Apologies if this was all assumed to be known in your original question.  Indeed, application of parallel direct solvers in PETSc often uses hyper/mumps/superlu, though I'm not familiar with ILU in these packages.<div><br></div><div>Ed</div><div><br clear="all"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Professor,<br>I'm writing about this mail about the ILU algorithm in PETSc. I have noticed that the ILU in PETSs cannot be paralleled, and the Hypre providing the parallel implement of ILU.<br>My questions is : why the PETSc does not provide the parallel implement of ILU directly, instead of linking to the Hypre to obtain the Parallel ILU algorithm?<br>PS: the attachment is the Parallel ILU algorithm in Hypre.<br>I'm looking forward to hear you!<br>Best,<br>Gang Chen<br>Sichuan Univertsity,<br>School of Mathematics,<br>Chengdu, China</blockquote><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Ed Bueler<br>Dept of Mathematics and Statistics<br>University of Alaska Fairbanks<br>Fairbanks, AK 99775-6660<br>306C Chapman<br></div></div></div></div></div></div></div></div></div>