[petsc-dev] MatHermitianTransposeGetMat

[Dave May]

Hi Pierre,

On Fri, 19 Oct 2018 at 21:28, Pierre Jolivet <pierre.jolivet at enseeiht.fr>
wrote:

> Hello,
> I’m trying to call MatTransposeGetMat on a Mat created with
> MatCreateHermitianTranspose, but this fails.
> Why is there no MatHermitianTransposeGetMat?


Looks like an oversight to me.


Is there a work around?


Doesn't seem to be a hack workaround possible as the object holding the
reference to the original matrix is declared directly in the C file. See
typedef struct Mat_HT

https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/transpose/htransm.c.html

Support for XXXGetMat() needs to go into master. Essentially one needs to
just copy the code from

https://www.mcs.anl.gov/pets/petsc-current/src/mat/impls/transpose/transm.c.html

associated with MatTransposeGetMat and change the type from Mat_Transpose
to Mat_HT.

Want to make a PR?

Cheers,
  Dave






>
> Thanks,
> Pierre
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20181019/eae2b1ec/attachment.html>