[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