Pascal,

The wrapped Trilinos matrices are based on Epetra which only supports
double AFAICT. That's why you can replace TrilinosScalar easily.
On the other hand, you should be able to compile PETSc with complex scalar
type and use that with MPI.

Best,
Daniel

Am Do., 23. Juli 2020 um 12:42 Uhr schrieb Pascal Kraft <
kraft.pas...@gmail.com>:

> Dear Deal.II devs and users,
>
> In the latest release a lot of (great) work has been done to make complex
> numbers more of a first-class citizen in deal, which has made my code a lot
> more readable. Currently, I am stuck with one problem, though. Are there
> any distributed datatypes for matrices that accept complex numbers?
>
> The dealii sparse matrix implementation is a template and allows complex
> numbers - however that implementation has no MPI functionality, which I
> need.
>
> The Petsc Sparse Matrix and Trilinos Sparse Matrix are no templates. In
> the types header I found the declaration of TrilinosScalar as double but
> changing it and recompiling dealii with the changed header threw an error.
>
> I have Trillions compiled with support for complex numbers and also
> searched through the LinearAlgebra documentation.
>
> 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.
>
> Can anyone give me a pointer on what is possible?
>
> Kind regards,
> Pascal Kraft
>
> --
> 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/4f4af1e0-4020-4dbf-aa9a-3a3fb7297d90o%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/4f4af1e0-4020-4dbf-aa9a-3a3fb7297d90o%40googlegroups.com?utm_medium=email&utm_source=footer>
> .
>

-- 
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/CAOYDWb%2BZZn-X6ZROtzYVX6bR1Qe8FXdxcueL%3DPL7qOFFmVhDHA%40mail.gmail.com.

Reply via email to