Pierre, I see that you've already done a merge request for this. Thanks! I have tested this and it works nicely in my application /Carl-Johan
On Wed, Apr 10, 2024 at 8:24 AM Carl-Johan Thore <carljohanth...@gmail.com> wrote: > > > On Tue, Apr 9, 2024 at 5:31 PM Pierre Jolivet <pie...@joliv.et> wrote: > >> >> On 9 Apr 2024, at 4:19 PM, Carl-Johan Thore <carljohanth...@gmail.com> >> 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!d3xUyKqZovzqGxg2_Ek4FStdTXsvnabp6mvY0kSRCCRZXZ6y6pd8MUnEi9sZR8QjNGiwWLxXy3N08YLxKwIcWBQ9Jbc-Qg$ >> 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 >