On Tue, Apr 9, 2024 at 5:31 PM Pierre Jolivet <[email protected]> wrote:
> > On 9 Apr 2024, at 4:19 PM, Carl-Johan Thore <[email protected]> > wrote: > > This Message Is From an External Sender > This message came from outside your organization. > Thanks for the suggestion. I don't have a factored matrix (and can't > really use direct linear solvers) so MatSolveTranspose doesn't seem to be > an option. > I should have mentioned that I've also tried KSPSolveTranspose but that > doesn't work with pcredistribute > > > I’m not a frequent PCREDISTRIBUTE user, but it looks like > https://urldefense.us/v3/__https://petsc.org/release/src/ksp/pc/impls/redistribute/redistribute.c.html*line332__;Iw!!G_uCfscf7eWS!cgAsxyEH_bcrB52QRg6oG5BdTF8bLQr9pb71JxdsoP7FOO0BCi8XGqL6D3oXT9GrnsuqWfA8MuLi8xxWvmzAKeSxmS6R2g$ > could > be copy/paste’d into PCApplyTranspose_Redistribute() by just changing a > MatMult() to MatMultTranspose() and KSPSolve() to KSPSolveTranspose(). > Would you be willing to contribute (and test) this? > Then, KSPSolveTranspose() — which should be the function you call — will > work. > > Thanks, > Pierre > Thanks, that sounds promising. Yes, I'll try to make a contribution /Carl-Johan
