Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Carl-Johan Thore
t;>> Kind regards, >>>> Carl-Johan >>>> >>>> >>>> On Sat, Aug 26, 2023 at 5:21 PM Pierre Jolivet via petsc-users < >>>> petsc-users@mcs.anl.gov> wrote: >>>> >>>>> >>>>> >>>>>

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Pierre Jolivet
ue.The branch should be usable in its current state.Thanks,Pierre From: Pierre Jolivet <pierre.joli...@lip6.fr> Sent: Saturday, August 26, 2023 4:36 PMTo: Carl-Johan Thore <carl-johan.th...@liu.se>Cc: Carl-Johan Thore <carljohanth...@gmail.com>; petsc-users@mcs.anl.govSubject: R

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Pierre Jolivet
> On Sat, Aug 26, 2023 at 5:21 PM Pierre Jolivet via petsc-users >>>>>>> wrote: >>>>>>> >>>>>>> >>>>>>>>> On 27 Aug 2023, at 12:14 AM, Carl-Johan Thore >>>>>>>>> wrote: >>&

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Carl-Johan Thore
en if symmetry in A00 and A11 cannot be exploited; using SBAIJ for >>> them would just be a >>> (pretty big) bonus. >>> >>> “I’ll rebase on top of main and try to get it integrated if it could be >>> useful to you (but I’m traveling >>> right n

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Pierre Jolivet
t;>>>> from these blocks, or is there >>>>>>> somewhere else in the code that things will go wrong? >>>>>> >>>>>> Ah, no, you are right, in that case it should work. >>>>>> >>>>>> For the MATNEST I was thinking

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-06 Thread Carl-Johan Thore
tsc/-/merge_requests/6841. >> I need to add a new code path in MatCreateRedundantMatrix() to make sure >> the resulting Mat is indeed SBAIJ, but that is orthogonal to the >> PCFIELDSPLIT issue. >> The branch should be usable in its current state. >> >> Thanks, >

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-09-04 Thread Pierre Jolivet
mcs.anl.govSubject: Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ  On 26 Aug 2023, at 11:16 PM, Carl-Johan Thore <carl-johan.th...@liu.se> wrote: "(Sadly) MATSBAIJ is extremely broken, in particular, it cannot be used to retrieve rectangular blocks in MatCreateSubMatrices, t

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-28 Thread Zhang, Hong via petsc-users
gt; Sent: Saturday, August 26, 2023 10:20 AM To: Carl-Johan Thore mailto:carl-johan.th...@liu.se>> Cc: petsc-users mailto:petsc-users@mcs.anl.gov>> Subject: Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ On 27 Aug 2023, at 12:14 AM, Carl-Johan Thore mailto:carl-johan.th...@liu.se>

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-28 Thread Pierre Jolivet
hogonal to the PCFIELDSPLIT >> issue. >> The branch should be usable in its current state. >> >> Thanks, >> Pierre >> >>> >>> From: Pierre Jolivet >> <mailto:pierre.joli...@lip6.fr>> >>> Sent: Saturday, August 26, 20

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-28 Thread Carl-Johan Thore
:20 AM > *To:* Carl-Johan Thore > *Cc:* petsc-users > *Subject:* Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ > > > > On 27 Aug 2023, at 12:14 AM, Carl-Johan Thore > wrote: > > “Well, your A00 and A11 will possibly be SBAIJ also, so you’ll end up with > the same issu

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-28 Thread Carl-Johan Thore
l to the > PCFIELDSPLIT issue. > The branch should be usable in its current state. > > Thanks, > Pierre > > > *From:* Pierre Jolivet > *Sent:* Saturday, August 26, 2023 4:36 PM > *To:* Carl-Johan Thore > *Cc:* Carl-Johan Thore ; petsc-users@mcs.anl.gov > *Sub

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-26 Thread Carl-Johan Thore via petsc-users
"(Sadly) MATSBAIJ is extremely broken, in particular, it cannot be used to retrieve rectangular blocks in MatCreateSubMatrices, thus you cannot get the A01 and A10 blocks in PCFIELDSPLIT. I have a branch that fixes this, but I haven’t rebased in a while (and I’m AFK right now), would you want

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-26 Thread Zhang, Hong via petsc-users
of smaller storage of SBAIJ (could at cost of communication overhead). Hong From: petsc-users on behalf of Pierre Jolivet via petsc-users Sent: Saturday, August 26, 2023 10:20 AM To: Carl-Johan Thore Cc: petsc-users Subject: Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-26 Thread Pierre Jolivet via petsc-users
2023 4:36 PM > To: Carl-Johan Thore > Cc: Carl-Johan Thore ; petsc-users@mcs.anl.gov > Subject: Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ > > > > > On 26 Aug 2023, at 11:16 PM, Carl-Johan Thore <mailto:carl-johan.th...@liu.se>> wrote: > > "

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-26 Thread Pierre Jolivet
> On 26 Aug 2023, at 11:16 PM, Carl-Johan Thore wrote: > > "(Sadly) MATSBAIJ is extremely broken, in particular, it cannot be used to > retrieve rectangular blocks in MatCreateSubMatrices, thus you cannot get the > A01 and A10 blocks in PCFIELDSPLIT. > I have a branch that fixes this, but I

Re: [petsc-users] PCFIELDSPLIT with MATSBAIJ

2023-08-26 Thread Pierre Jolivet
(Sadly) MATSBAIJ is extremely broken, in particular, it cannot be used to retrieve rectangular blocks in MatCreateSubMatrices, thus you cannot get the A01 and A10 blocks in PCFIELDSPLIT. I have a branch that fixes this, but I haven’t rebased in a while (and I’m AFK right now), would you want me