This is definitely wrong. Dave and I have been doing this all week here and those times are way way off.
This looks to me like ML on the entire Stokes system. There is no way that is ML on Poisson (the 0,0 block). Matt On Thu, Sep 25, 2008 at 6:34 AM, Nuno David Lopes <[EMAIL PROTECTED]> wrote: >> Have you tried with ML (-pc_type >> ml -mat_no_inode)? Multigrid should work very well for this problem. > It works! > _________________________________________________ > Solving linear system of size 566769 x 566769 (Krylov solver). > PETSc Krylov solver (gmres, ml) converged in 9 iterations. > End velocity Init Step after 113.390000 seconds. > _________________________________________________ > > SOR still is faster > ___________________________________________________ > Solving linear system of size 566769 x 566769 (Krylov solver). > PETSc Krylov solver (gmres, sor) converged in 35 iterations. > End velocity Init Step after 69.530000 seconds. > ___________________________________________________ > >> Jed > > > > -- > Nuno David Lopes > > e-mail:[EMAIL PROTECTED] (FCUL/CMAF) > [EMAIL PROTECTED] (ISEL) > http://ptmat.ptmat.fc.ul.pt/%7Endl/ > > Thu Sep 25 12:30:54 WEST 2008 > _______________________________________________ > DOLFIN-dev mailing list > [email protected] > http://www.fenics.org/mailman/listinfo/dolfin-dev > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
