Complex Number Compatibility?

[Heatherms27]

Are there any plans to add complex number capabilities for different solvers/preconditioners? More specifically, I was interested in reading in a file containing a CSR matrix with complex values to use Euclid.

[victorapm]

Hi @Heatherms27, there aren't plans now for adding support of complex numbers to preconditioners in the IJ interface. Sorry about that! Currently, many compilation warnings are issued when --enable-complex is used and, in general, they point to locations where algorithms need to be rethought to make them work with complex numbers. As a matter of fact, the ILU implementations in hypre (Euclid and hypre_ILU) are not supposed to deal with complex numbers.

As for reading complex matrices from file, If you are interested, I can point to the modifications needed in HYPRE_IJMatrixRead to make it work for this case.

Thanks!

[lucasbanting]

@victorapm Do you know if Hypre will/currently support complex matrices with zero imaginary part?

I would like to use AMS with my existing complex PETSc code, but want to avoid having to manually split my complex petsc matrix into multiple hypre real matrices.