[Dune] [#1018] RFC Patch: allow different convergence criteria in dune-istl (Attachment added)
Dune
flyspray at dune-project.org
Mon Jan 16 17:23:19 CET 2012
THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.
A new Flyspray task has been opened. Details are below.
User who did this - Andreas Lauser (andlaus)
Attached to Project - Dune
Summary - RFC Patch: allow different convergence criteria in dune-istl
Task Type - Feature Request
Category - Solver
Status - Unconfirmed
Assigned To -
Operating System - Unspecified / All
Severity - Low
Priority - Normal
Reported Version - 2.1
Due in Version - Undecided
Due Date - Undecided
Details - Hi,
Although we at LH2 are generally really happy with the solvers of Dune-ISTL, we would like to use a different convergence criterion than the reduction of the two-norm of the residual. I've prepared a proof-of patch which makes the criterion pluggable and did some measurements with Dumux. In the end, the whole simulation was about 25 % faster (see the attached logfiles).
note, that this patch is currently more a proof of concept against dune 2.1 and I would appreciate some feedback. the patch description is the following:
it works by adding a template argument 'ConvergenceCriterion' to the
solver. The default criterion is to look at the reduction of the
defect for the current iterative solution. This is the behaviour of
the ISTL solvers before this patch, so no code needs to be changed if
the current behaviour ought to be kept. The patch adds a second
criterion which checks whether the weighted difference between two
iterations is below a given threshold. This second criterion is called
'FixPointCriterion' and shows better performance (at least for Dumux,
where the Newton-Raphson solver uses a relative convergence criterion).
So far, only the BiCGStab solver has been modified. If this patch has
any chance of getting merged into the Dune trunk, I will also modify
the other solvers.
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=1018
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