Dear developers, For the Navier-Stokes equations, there is a null space in the block matrix (right-bottom block) -- -- | -laplace grad | | div 0 | -- --
I use bjacobi+ilu as the preconditioner for GMRES in PETSc to solve the system. Sometimes it works, but most of time it does not converges. My code is just like the example: /systems_of_equations_ex2. I wonder if I try the way in example: /systems_of_equations_ex3, would it be better? Or should I use other solver/perconditioner ? I am appreciated for your suggestions. -- Li Luo ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports.http://sdm.link/zohodev2dev _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users