[petsc-users] MatIsSymmetric / MatIsHermitian issues for MPI+complex
Barry Smith
bsmith at mcs.anl.gov
Mon Mar 7 13:14:32 CST 2016
If you know the matrix is symmetric or Hermitian you can call MatSetOption() to tell the matrix. Then it doesn't need to try to check.
Barry
> On Mar 7, 2016, at 4:26 AM, Denis Davydov <davydden at gmail.com> wrote:
>
> Nevermind, i should have check the error code:
>
> 45: [1]PETSC ERROR: No support for this operation for this object type
> 45: [1]PETSC ERROR: Matrix of type <mpiaij> does not support checking for symmetric
>
>
>> On 7 Mar 2016, at 11:21, Denis Davydov <davydden at gmail.com> wrote:
>>
>> Dear all,
>>
>> I call MatIsSymmetric and MatIsHermitian for MPI complex-valued matrix in PETSc.
>> For a moment, i store a mass matrix (real-valued, symmetric) in a matrix.
>> However, the result of both functions is NOT PETSC_TRUE.
>> The same program works fine with a single MPI core.
>>
>> Are there any known issues of MatIsSymmetric / MatIsHermitian for COMPLEX+MPI matrices?
>> I am quite certain that my matrix is indeed symmetric (i also test vmult and Tvmult on a vector and the results is the same).
>>
>> Kind regards,
>> Denis
>>
>
More information about the petsc-users
mailing list