[Dune] [#729] Singularity detection in FieldMatrix' LU decomposition buggy

Dune flyspray at dune-project.org
Thu Jan 28 14:21:49 CET 2010


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

The following task has a new comment added:

FS#729 - Singularity detection in FieldMatrix' LU decomposition buggy
User who did this - Sven Marnach (sven)

----------
Andreas, you are right, there seem to be some issues with the pivoting in the LU decomposition.  In particular. it would be nice to have full pivoting instead of row pivoting only.  Furthermore, the singularity threshold should depend on the data type of the matrix entries (at the moment, it is more or less hard-coded).

But your matrix is very ill-conditioned (condition ~ 10^19) -- it is not a complete surprise it is rejected by the algorithm.

-- Christian and Sven
----------

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

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.




More information about the Dune mailing list