[petsc-users] About compatibility with C++ complex type

[Jed Brown]

yuqing xia <yuqing.xia at colorado.edu> writes:

> Dear Developer and Users
>
>  I am new to SLEPC. So my PETSC is compiled with flag
> --with-scalar-type=complex --with-clanguage=c++, therefore SLEPC should
> also support c++ complex type. But when I tried to insert C++ complex type
> value into a PETSC complex matrix, It seems SLEPC's inside complex type is
> not compatible with C++ complex type according to the test result.
>
> Here is the matrix I want to insert their value to Petsc matrix
> (-0.131541,0)(0,0)(0,0)
> (0,0)(-0.188541,0)(0,0)
> (0,0)(0,0)(-0.245541,0)
> Here is the Petsc matrix after inserting the above matrix value into it
> Matrix Object: 1 MPI processes
> type: seqaij
> row 0: (0, -0.131541)  (1, 0)  (2, 0)
> row 1: (0, 0)  (1, -0.188541)  (2, 0)
> row 2: (0, 0)  (1, 0)  (2, -0.24554

These are (column_index, matrix_entry) pairs.

> Does anyone know how to avoid the problem.
>
> Best
> Yuqing Xia
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140426/393fd3cf/attachment.pgp>