> In theory if you use petsc-3.3 (release any day) switching to use > the PCFIELDSPLIT solver (with appropriate options) should be able > to reproduce your algorithm (I think) with very little coding on > your part. This would make it easier for us to help you understand > the issue since we know exactly what PCFIELDSPLIT does and how it > is suppose to work and the -ksp_view option (or -snes_view) tells > exactly the algorithm being used. > > Barry
Yes, I've been thinking about that, it should be possible without too much effort. At this point I just don't have any time but I'll keep it in mind. Thanks to Jed and to you for your answers! dr. ir. Christiaan Klaij CFD Researcher Research & Development E mailto:C.Klaij at marin.nl T +31 317 49 33 44 MARIN 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl
