[Dune] DUNE::SeqILU0 craches for some matrices

Thomas Witkowski thomas.witkowski at tu-dresden.de
Fri Feb 8 15:56:40 CET 2008


Okay, this was a very stuipd example :)  But it helped me to find a 
problem in my program. So there is no bug in DUNE at this point!

Regards,

Thomas

Markus Blatt wrote:
> 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
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
>
>
>   





More information about the Dune mailing list