[petsc-users] [KSP] solveTranspose fails with Strumpack and SuperLU_dist

Matthew Knepley knepley at gmail.com
Sat Apr 9 18:49:52 CDT 2022

On Sat, Apr 9, 2022 at 7:41 PM Jean Marques <jeanmarques.eng at gmail.com>

> Hi all,
> This may be a naive question, and I hope this is the right place to ask
> about it.
> I need to solve a direct linear system with a sparse matrix R, then an
> adjoint system the hermitian of R.
> I use a petsc4py, so what I do is this:
> self.R.setUp()
> to set up the PETSc KSP variable R, then I do:
> self.R.solve(f, q)
> and later:
> self.R.solveTranspose(f, q)
> 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?

Yes, they do not support transpose solves.

> or is there a specific way to call them with these libraries?
> Maybe the approach would be to define another 'self.R' variable but then
> set it as a transpose().conjugate() before setUp()?
> I was trying STRUMPACK because it has a low-memory approach and that's my
> main bottleneck.

What kind of a system are you solving?



> I appreciate any help you can provide.
> Best,
> Jean
> --
> *Jean Helder Marques Ribeiro*
> Ph.D. Candidate
> University of California, Los Angeles
> 420 Westwood Plaza, Los Angeles, CA 90095
> phone: (310) 689-6593

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220409/cf2d51b3/attachment.html>

More information about the petsc-users mailing list