Pascal,
Originally, I only knew Trilinos because I used the distributed matrices and
vectors in my application. I also knew that there is a configuration of
trilinos to make complex numbers available in all packages that support it.
However, from what I can tell, that only effects Tpetra
The documentation states that Tpetra supports
- MPI
- Shared memory parallelization: (OpenMP, CUDA, Posix)
and:
Scalar: A Scalar is the type of values in the sparse matrix or dense
vector. This is the type most likely to be changed by many users. The most
common use cases
are float, double,
Hi Wolfgang,
here is what I found out about the topic:
Originally, I only knew Trilinos because I used the distributed matrices
and vectors in my application. I also knew that there is a configuration of
trilinos to make complex numbers available in all packages that support it.
However, from
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
Hi Daniel,
oh, I'm really sorry for asking if that works. I had seen that neither
PETSc nor Trilinos Sparse Matrices are templated and assumed that if the
more modern version (Trilinos) doesn't work with complex numbers, trying
PETSc wouldn't be very promising. But you are right, I will try
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
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