[Dune] Complex systems

[Christian Engwer]

Dear Dimitrios,

> I am currently looking for a library that handles adaptive mesh
> refinement fem and solves complex linear systems with iterative
> methods, preferably cg type methods. Does anyone know if DUNE
> supports the numerical solution of complex linear systems? (complex
> system refers to a linear system where the matrix, the right hand
> side and the solution are matrices with complex numbers).

Dune has interfaces to a range of FEM meshed, e.g. the alugrid which
is a locally adaptive hexahedral mesh. In general you can interface
any LA library you want. As part of Dune we develop dune-istl, a
library for hierarchic matrices and iterative solvers. These matrices
do also support complex and most of the solvers should work with
complex values. There are only very few people using Dune with complex
and a few bugs had to be fixed recently, but the idea is to have
everything working with complex and the usual stuff does work.

Christian

-- 
Prof. Dr. Christian Engwer 
Institut für Numerische und Angewandte Mathematik
Fachbereich Mathematik und Informatik der Universität Münster
Einsteinstrasse 62
48149 Münster

E-Mail	christian.engwer at uni-muenster.de
Telefon	+49 251 83-35067
FAX		+49 251 83-32729