Wolfgang, I have convergence to the correct solution for very coarse meshes (I have an analytic free-space solution that I can compare to), and I have also been able to solve the problem with a direct solver. The issue of non-convergence crops up when I use a finer mesh, and as a result, have a larger system to solve. I have tried to estimate the condition number of this larger system with ksp_monitor_singular_value, and it seems reasonable; so it still seems like the preconditioner is just not good enough.
There are papers that already address preconditioning for PML, but since multigrid preconditioning is a black-box method for me, I wanted to ask if anyone can suggest a good approach that is implementable in deal. Artur -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
