<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? 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><br></div><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" class="gmail_signature" data-smartmail="gmail_signature"><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>