[Dune] [#928] breakdown of BiCGSTAB (Attachment added)

[Dune]

THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

A new Flyspray task has been opened.  Details are below. 

User who did this - Bernd Flemisch (Bernd) 

Attached to Project - Dune
Summary - breakdown of BiCGSTAB
Task Type - Bug Report
Category - Dune Core Modules
Status - Unconfirmed
Assigned To - 
Operating System - Unspecified / All
Severity - Low
Priority - Normal
Reported Version - SVN (pre2.2)
Due in Version - Undecided
Due Date - Undecided
Details - Dear Dune, 

in some cases, the BiCGSTABSolver breaks down with the message 
Dune reported error: ISTLError [apply:../../../dune/istl/solvers.hh:791]: h=0 in BiCGSTAB

I attach an example with a particular case, to be run with 
./test_bicgstab matrix.txt 

I checked the matrix with Matlab, it is a perfectly fine positive definite matrix. As can be seen from the example, the LoopSolver works (so does the CGSolver). The preconditioner does not seem to matter. 

What seems to matter, is the right hand side. If I change line 82 to 
b = 1;
everything is fine. 

It would be great if you could check. Thank you! 

BTW: The example contains a routine 
void readBCRSMatrixFromMatlabFile(const char *fileName, Matrix& A)
which I think could be interesting for inclusion in dune/istl/io.hh. It is the reverse to writeMatrixToMatlab.


One or more files have been attached.

More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=928

You are receiving this message because you have requested it from the Flyspray bugtracking system.  If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above.