[Dune] Applying an ISTL InverseOperator to multiple right-hand sides

[Markus Blatt]

Hi,

On Mon, May 03, 2010 at 05:52:58PM +0200, Bård Skaflestad wrote:
> 
> > Hi Bård!
> > 
> > Assembling the matrix once and then solving for all rhs at once (e.g.
> > with a BCRSMatrix<ScaledIdMatrix>) may(!) still be faster, because you
> > get a better cache usage.
> 
> Hi Oliver,
> 
> The ScaledIdentityMatrix is indeed a possible solution.  I'll
> investigate its performance characteristics in my case.  

While AMG compiles if SuperLU is deactivated, not using a direct
solver on the coarse level might present a problem concerning
convergence.

Therefore one should consider the time to adapt SuperLU to the matrix
type.

Cheers,

Markus