<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Nov 22, 2013 at 12:24 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Michael Povolotskyi <<a href="mailto:mpovolot@purdue.edu">mpovolot@purdue.edu</a>> writes:<br>
<br>
> This sounds very reasonable for me and my colleagues to have<br>
> MatSetTransposeType(A,MAT_TRANSPOSE_HERMITIAN).<br>
> Our problems always require Hermitian transpose but I assume there are<br>
> problems where it is not the case.<br>
<br>
</div>Complex symmetric shows up in frequency-domain wave problems.  Methods<br>
like Sou-Cheng's CS-MINRES take advantage of that structure.<br>
<div class="im"><br>
> Actually my goal is to solve the transposed problem with PARDISO<br>
> package. It looks like it is the fastest option for the problems we<br>
> usually solve.<br>
><br>
> By the way, it there a chance to have interface to Pardiso in the next<br>
> PETSc release? My colleague Jose David Bermeo sent the modules for this<br>
> some time ago for your approval.<br>
<br>
</div>I was buried in a proposal deadline and then SC13, but I will go through<br>
the pull requests that nobody has gotten to.  I think the Pardiso<br>
contribution is more-or-less ready to merge.  For MatMatSolveTranspose,<br>
I would recommend doing it in a new branch starting from 'master'.<br>
Please include a reference implementation so that we can test without<br>
Pardiso (which is not free software).  The code should be similar to<br>
MatMatSolve_SeqAIJ.  Once that works (and has a test), you'll merge the<br>
Pardiso branch and add MatMatSolveTranspose_Pardiso.<br></blockquote><div>By analogy with MatTransposeMatMult() etc., should this be MatTransposeMatSolve()?</div><div>Since 'Transpose' applies to the first Mat.</div>

<div>Dmitry. </div></div><br></div><br clear="all"><div><br></div></div>