[Dune] Solving Ax=b with A being a matrix of sparse matrices

Markus Blatt Markus.Blatt at ipvs.uni-stuttgart.de
Tue Mar 4 15:04:35 CET 2008


On Tue, Mar 04, 2008 at 02:31:38PM +0100, Thomas Witkowski wrote:
> Markus Blatt wrote:
> > The problem is that BCRSMatrix has no solve method. One has to set the
> > blocklevel of the preconditioners to l=2 instead of 1.
> >   
> And how can I do this?
> 

I would have loved to RTFM, but it does not seem as fine as it ought
to be.

E.g. use
SeqSSOR<Matrix,Vector,Vector,2>
instead of 
SeqSSOR<Matrix,Vector,Vector>
(which defaults to
SeqSSOR<Matrix,Vector,Vector,1>)
as a preconditioner.

Markus




More information about the Dune mailing list