<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Dave,<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 19 Oct 2018, at 11:10 PM, Dave May <<a href="mailto:dave.mayhem23@gmail.com" class="">dave.mayhem23@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class=""><div class=""><div dir="auto" class="">Hi Pierre,</div></div><div class=""><br class=""><div class="gmail_quote"></div></div></div><div class=""><div dir="ltr" class="">On Fri, 19 Oct 2018 at 21:28, Pierre Jolivet <<a href="mailto:pierre.jolivet@enseeiht.fr" target="_blank" class="">pierre.jolivet@enseeiht.fr</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br class="">
I’m trying to call MatTransposeGetMat on a Mat created with MatCreateHermitianTranspose, but this fails.<br class="">
Why is there no MatHermitianTransposeGetMat?</blockquote><div dir="auto" class=""><br class=""></div></div><div class=""><div dir="auto" class="">Looks like an oversight to me.</div></div></div></blockquote><div><br class=""></div><div>That is what I was afraid of. I’ve proposed a pull request as you suggested.</div><div><br class=""></div><div>Thanks,</div><div>Pierre</div><br class=""><blockquote type="cite" class=""><div class=""><div class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> Is there a work around?</blockquote><div dir="auto" class=""><br class=""></div></div><div class=""><div dir="auto" class="">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</div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><div class=""><a href="https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/transpose/htransm.c.html" target="_blank" class="">https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/transpose/htransm.c.html</a></div></div><div dir="auto" class=""><br class=""></div><div dir="auto" class="">Support for XXXGetMat() needs to go into master. Essentially one needs to just copy the code from </div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><div class=""><a href="https://www.mcs.anl.gov/pets/petsc-current/src/mat/impls/transpose/transm.c.html" target="_blank" class="">https://www.mcs.anl.gov/pets/petsc-current/src/mat/impls/transpose/transm.c.html</a></div><br class=""></div><div dir="auto" class="">associated with MatTransposeGetMat and change the type from Mat_Transpose to Mat_HT.</div><div dir="auto" class=""><br class=""></div><div dir="auto" class="">Want to make a PR?</div><div dir="auto" class=""><br class=""></div><div dir="auto" class="">Cheers,</div><div dir="auto" class=""> Dave</div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><br class=""></div><div dir="auto" class=""><br class=""></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br class="">
<br class="">
Thanks,<br class="">
Pierre</blockquote>
</div>
</div></blockquote></div><br class=""></body></html>