[Dune] [#928] breakdown of BiCGSTAB

[Markus Blatt]

On Thu, Jun 09, 2011 at 02:42:54PM +0200, Christian Engwer wrote:
> Hi Bernd,
> 
> did you try reducing the breakdown criterion? In other implementations
> it is '== 0' instead of '< epsilon'. Does it work then, read "it is
> just the choice of the epsilon", or does it still fail, read "there
> is a real bug in the BiCGStab"?

h is actually zero! And even if it were just really small, rounding
errors would be dramatic.

Markus

-- 
Markus Blatt, Interdisciplinary Center for Scientific Computing,
University Heidelberg, Im Neuenheimer Feld 368, D-69120 Heidelberg

'Mathematics knows no races or geographic boundaries; for
mathematics, the cultural world is one country' - David Hilbert