Hi Pierre, On Fri, 19 Oct 2018 at 21:28, Pierre Jolivet <[email protected]> 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
