[petsc-users] Binary format in real vs. complex scalar type configurations

[Mark Adams]

On Tue, May 11, 2021 at 9:30 AM Frederico Teixeira <teixeira at zmt.swiss>
wrote:

> Dear fellows,
>
> I hope this message finds you safe and well.
>
> I have a complex-valued matrix and its real/imaginary components in binary
> format. They were extracted from a solver that only works with
> "scalar-type=complex" configuration.
> I am getting weird results when I load them into a small test program
> that's configured with "scalar-type=real", but I believe this is expected.
> At the end of the day, I would like to have both real and imaginary
> components as real-valued matrices.
> Is it possible to do it? I want to test preconditioners that are tailored
> for this sort of problem.
>

Do you mean you want what is called "equivalent real form" where the real
and complex parts are stored as type 'double' for example and operations
like multiply take two pairs of doubles, do a complex multiply manually,
and return a real/complex pair of doubles?


>
> Regards,
>  Frederico.
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210511/ae69fb4b/attachment.html>