<div dir="ltr"><div dir="ltr">On Sat, Apr 9, 2022 at 7:41 PM Jean Marques <<a href="mailto:jeanmarques.eng@gmail.com">jeanmarques.eng@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div></div><div>This may be a naive question, and I hope this is the right place to ask about it.</div><div>I need to solve a direct linear system with a sparse matrix R, then an adjoint system the hermitian of R.</div><div><br></div><div>I use a petsc4py, so what I do is this:</div><div>self.R.setUp()</div><div>to set up the PETSc KSP variable R, then I do:</div><div>self.R.solve(f, q)</div><div>and later:</div><div>self.R.solveTranspose(f, q)<br></div><div><br></div><div>However, 'solveTranspose()' only works when I use MUMPS. If I try STRUMPACK or SUPERLU_DIST it fails, it seems that 'solveTranpose()' is not defined for them?</div></div></blockquote><div><br></div><div>Yes, they do not support transpose solves.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div> or is there a specific way to call them with these libraries?</div><div>Maybe the approach would be to define another 'self.R' variable but then set it as a transpose().conjugate() before setUp()?</div><div>I was trying STRUMPACK because it has a low-memory approach and that's my main bottleneck.</div></div></blockquote><div><br></div><div>What kind of a system are you solving?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>I appreciate any help you can provide.</div><div><br></div><div>Best,</div><div>Jean<br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px"><div dir="ltr"><div style="font-size:12.8px"><div><b><i>Jean Helder Marques Ribeiro</i></b><br></div><font color="#000000">Ph.D. Candidate</font></div><div><font color="#000000">University of California, Los Angeles</font></div><font color="#000000">420 Westwood Plaza, Los Angeles, CA 90095<br></font></div><div><span style="font-size:12.8px"><font color="#000000">phone: (310) 689-6593</font></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>