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

[Matthew Knepley]

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

> 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?

  Thanks,

     Matt


> 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>