Hi there, you may also like to check out the "PetscExt" project that a colleague of mine has developed over the last couple of years for solving block matrix systems in PETSc: http://jupiter.ethz.ch/~dmay/Research/PetscExt/index.html -- Patrick. On Sat, Oct 2, 2010 at 3:28 PM, Jed Brown <[email protected]> wrote:
2010/10/2 Gong Ding <[email protected]>:
Dear all, I had take a look with the manual of petsc.It has refered schur complement. However, I can not find more useful information.
A couple useful references, there are others depending on your application.
@article{benzi2005nss, title={{Numerical solution of saddle point problems}}, author={Benzi, M. and Golub, G.H. and Liesen, J.}, journal={Acta Numerica}, volume={14}, pages={1--137}, year={2005}, publisher={Cambridge Univ Press} }
@article{elman2008tcp, title={{A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier-Stokes equations}}, author={Elman, H.C. and Howle, V.E. and Shadid, J. and Shuttleworth, R. and Tuminaro, R.}, journal={Journal of Computational Physics}, volume={227}, number={1}, pages={1790--1808}, year={2008}, publisher={Academic Press} }
Many of these methods are easy to implement using PCFieldSplit, when you find one that is likely to work well for your application, ask and we can suggest a good way to implement it.
Jed
-- Patrick Sunter VPAC Senior Computational Software Developer - AuScope Monash SAM Project Monash University Adjunct Research Associate (School of Mathematical Sciences) Room 301, Building 28 Monash University VIC 3800 Australia Ph: +61 (0)3 9905 4468 For VPAC/AuScope admin issues email: [email protected] For AuScope software development email: [email protected]