On 7/23/20 10:42 AM, Pascal Kraft wrote:

I have Trillions compiled with support for complex numbers and also searched through the LinearAlgebra documentation.

I don't think I knew that one can compile Trilinos with complex numbers. How do you do that?

It does not greatly surprise me that we use TrilinosScalar and double interchangeably. If Trilinos can indeed be compiled with complex numbers, then we ought to find a way to (i) make TrilinosScalar dependent on whatever Trilinos was compiled for, (ii) ensure that all of the places that currently don't compile because we use double in place of TrilinosScalar are fixed.

Patches are, as always, very welcome!


I require GMRES as a solver (which should be possible, because the GMRES Versions all use a templated Vector which can take complex components) and MPI distribution of a sparse system. I have so far only seen FullMatrix to accept complex numbers.

I believe that GMRES could indeed be made to work for complex-valued problems, but I'm not sure any of us have every tried. When writing step-58, I toyed with the idea of looking up in the literature what one would need for a complex GMRES, but in the end decided to just make SparseDirectUMFPACK work instead. The issue is that for every matrix-vector and vector-vector operation that happens inside GMRES, you have to think about whether one or the other operand needs to be complex-conjugated. I'm certain that is possible, but would require an audit of a few hundred lines. It would probably be simpler to just use PETSc's (or Trilinos') GMRES implementation.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/0a11242b-16e1-7c2b-681e-398602d2c60c%40colostate.edu.

Reply via email to