dolfin team mailing list archive

["Matthew Knepley" <knepley@xxxxxxxxx>]

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 <ndl@xxxxxxxxxxxxxx> 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:ndl@xxxxxxxxxxxxxx        (FCUL/CMAF)
>           nlopes@xxxxxxxxxxxxxxx    (ISEL)
> http://ptmat.ptmat.fc.ul.pt/%7Endl/
>
> Thu Sep 25 12:30:54 WEST 2008
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> 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