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

Thomas Witkowski thomas.witkowski at tu-dresden.de
Tue Mar 4 14:31:38 CET 2008


Markus Blatt wrote:
> On Tue, Mar 04, 2008 at 01:43:29PM +0100, Thomas Witkowski wrote:
>   
>> Is it possible to use the dune solvers to solve a linear system of 
>> equations where A is a small (2x2 or 3x3) matrix of large sparse 
>> matrices? I've created a matrix of type Dune::Matrix<  Dune::BCRSMatrix< 
>> Dune::FieldMatrix<double, 1, 1> > >, but the program does not compile. 
>> It seems that the type Dune::Matrix<T> is never thought to be used in 
>> Dune::Preconditioner and Dune::InverseOperator. Is there another 
>> matrix-type I can use for this purpose?
>>
>>     
>
> 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?

Thomas




More information about the Dune mailing list