[Dune] DUNE::SeqILU0 craches for some matrices

[Markus Blatt]

On Wed, Feb 06, 2008 at 02:30:55PM +0100, Thomas Witkowski wrote:
> 
> for some matrices DUNE::SeqILU0 crashes the whole machine with a 
> segmentation fault. I've made a very small example to show the problem. 
> Why it works for matrix B but not for matrix A? I'm using DUNE 1.0beta7
> 
>   A[0][0] = 1;
>   A[1][0] = 2;
>   A[2][0] = 3;
>   A[3][0] = 4;
>   A[3][1] = 5;
> 
A = [ 1 0 0 0
      2 0 0 0
      3 0 0 0
      4 5 0 0]

Matrix A has rank 2 and zeros on the diagonal! What would you expect?

Probably the code should throw a meaningfull exception here.

Markus