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 <jed at 59a2.org> wrote:
> 2010/10/2 Gong Ding <gdiso at ustc.edu>:
>> 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: patrick at vpac.org
For AuScope software development email: patdevelop at gmail.com
