Hi Timo and Wolfgang It turns out the SparseILU needs to be tuned very carefully. I set the AdditioalData to (0.03, 0) the solver now works fine. Although it is still slower than UMFPACK. Maybe using just the right parameters in AdditionalData can further improve it. I observed similar thing in step-57 too. Using ILU + GMRES for A inverse won't even converge in reasonable time with default AdditionalData. Playing with the AdditionalData helps... In a word, SparseILU seems to be sensitive to the AdditionalData and how to choose the right parameter is unclear.
Wolfgang, would you like me to contribute this code to the code gallery? Is there anything in this code that I should polish? Thank you Jie -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
