Re: [petsc-dev] Right-preconditioned GMRES

https://gitlab.com/petsc/petsc/merge_requests/2290 > On Nov 7, 2019, at 4:24 AM, Pierre Jolivet wrote: > > > >> On 7 Nov 2019, at 5:32 AM, Smith, Barry F. wrote: >> >> >> Some idiot logged what they did, but not why they did it. >> >> commit bf108f309acab50613e150419c680842cf4b8a05

Re: [petsc-dev] Right-preconditioned GMRES

> On 7 Nov 2019, at 5:32 AM, Smith, Barry F. wrote: > > > Some idiot logged what they did, but not why they did it. > > commit bf108f309acab50613e150419c680842cf4b8a05 (HEAD) > Author: Barry Smith > Date: Thu Mar 18 20:40:53 2004 -0600 > >bk-changeset-1.2063.1.1 >

Re: [petsc-dev] Right-preconditioned GMRES

"Smith, Barry F. via petsc-dev" writes: > Some idiot logged what they did, but not why they did it. > > commit bf108f309acab50613e150419c680842cf4b8a05 (HEAD) > Author: Barry Smith > Date: Thu Mar 18 20:40:53 2004 -0600 > > bk-changeset-1.2063.1.1 >

Re: [petsc-dev] Right-preconditioned GMRES

Some idiot logged what they did, but not why they did it. commit bf108f309acab50613e150419c680842cf4b8a05 (HEAD) Author: Barry Smith Date: Thu Mar 18 20:40:53 2004 -0600 bk-changeset-1.2063.1.1 barrysmith@barry-smiths-computer.local|ChangeSet|20040319024053|12244 ChangeSet

Re: [petsc-dev] Right-preconditioned GMRES

> On 6 Nov 2019, at 7:44 PM, Matthew Knepley wrote: > > On Wed, Nov 6, 2019 at 1:32 PM Pierre Jolivet via petsc-dev > mailto:petsc-dev@mcs.anl.gov>> wrote: >> On 6 Nov 2019, at 7:30 PM, Fande Kong > > wrote: >> >> Even if both the linear operator and the

Re: [petsc-dev] Right-preconditioned GMRES

> On 6 Nov 2019, at 7:30 PM, Fande Kong wrote: > > Even if both the linear operator and the preconditioner operator are > symmetric, the preconditioned operator by RAS will NOT be symmetric so that > CG is not applicable. > > RAS keeps overlapping elements during restriction but discards

Re: [petsc-dev] Right-preconditioned GMRES

I need to figure this out myself first. In the meantime, here is another (PCASM related) question for you: why is PCASMType switched to BASIC when using a symmetric Pmat? (I guess I don’t have to tell you about the performance of RAS vs. ASM) To me, that would make sense if the default KSP was

Re: [petsc-dev] Right-preconditioned GMRES

Send the code and exact instructions to run a "good" and a "bad" ASM Barry > On Oct 14, 2019, at 10:44 AM, Pierre Jolivet > wrote: > > Here are the three logs. > FGMRES also gives a wrong first iterate. > I think Mark was right in the sense that the problem is _most likely_ in my >

Re: [petsc-dev] Right-preconditioned GMRES

Here are the three logs. FGMRES also gives a wrong first iterate. I think Mark was right in the sense that the problem is _most likely_ in my RHS. But I need to figure out why I only get this problem with right-preconditioned KSPs with restrict or none. Thanks, Pierre dump-basic Description:

Re: [petsc-dev] Right-preconditioned GMRES

"Smith, Barry F. via petsc-dev" writes: > Is this one process with one subdomain? (And hence no meaningful overlap > since there is nothing to overlap?) And you expect to get the "exact" answer > on one iteration? > > Please run the right preconditioned GMRES with -pc_asm_type [restrict

Re: [petsc-dev] Right-preconditioned GMRES

Is this one process with one subdomain? (And hence no meaningful overlap since there is nothing to overlap?) And you expect to get the "exact" answer on one iteration? Please run the right preconditioned GMRES with -pc_asm_type [restrict and basic and none] -ksp_monitor_true_solution

Re: [petsc-dev] Right-preconditioned GMRES

> On 13 Oct 2019, at 3:48 PM, Mark Adams wrote: > > Oh I see you have established that it is the overlap and got something that > looks like the load on a homogenous BC and the processor boundary ... very > odd. > > Does the bad version converge to the correct solution if you iterate? It >

Re: [petsc-dev] Right-preconditioned GMRES

Thanks for reporting the issue and your patience. I have made a fix for this problem but need to have Jed review it: https://bitbucket.org/petsc/petsc/pull-requests/634/allow-users-to-change-the-pcside-and/diff If Jed is happy with it we'll get it into master within a couple of days.

Re: [petsc-dev] Right-preconditioned GMRES

I forgot to mention this, but you just need to run: mpirun -np 1 ./ex1 Can you reproduce the issue? Thanks. > On Feb 9, 2017, at 10:52 PM, Pierre Jolivet > wrote: > > There, attached. > > Thanks. > > PS: here is the patch also, just in case > diff --git

Re: [petsc-dev] Right-preconditioned GMRES

There, attached.Thanks.PS: here is the patch also, just in casediff --git a/src/ksp/ksp/examples/tutorials/ex1.c b/src/ksp/ksp/examples/tutorials/ex1.cindex 6382c09..d9dc69a 100644--- a/src/ksp/ksp/examples/tutorials/ex1.c+++ b/src/ksp/ksp/examples/tutorials/ex1.c@@ -123,6 +123,11 @@ int main(int

Re: [petsc-dev] Right-preconditioned GMRES

On Thu, Feb 9, 2017 at 7:48 AM, Pierre Jolivet wrote: > On Feb 9, 2017, at 2:31 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 7:29 AM, Pierre Jolivet > wrote: > >> >> On Feb 9, 2017, at 2:17 PM, Matthew Knepley

Re: [petsc-dev] Right-preconditioned GMRES

> On Feb 9, 2017, at 2:31 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 7:29 AM, Pierre Jolivet > wrote: > >> On Feb 9, 2017, at 2:17 PM, Matthew Knepley >

Re: [petsc-dev] Right-preconditioned GMRES

On Thu, Feb 9, 2017 at 7:29 AM, Pierre Jolivet wrote: > > On Feb 9, 2017, at 2:17 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 7:06 AM, Pierre Jolivet > wrote: > >> >> On Feb 9, 2017, at 1:37 PM, Matthew

Re: [petsc-dev] Right-preconditioned GMRES

> On Feb 9, 2017, at 2:17 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 7:06 AM, Pierre Jolivet > wrote: > >> On Feb 9, 2017, at 1:37 PM, Matthew Knepley >

Re: [petsc-dev] Right-preconditioned GMRES

On Thu, Feb 9, 2017 at 7:06 AM, Pierre Jolivet wrote: > > On Feb 9, 2017, at 1:37 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 4:56 AM, Pierre Jolivet > wrote: > >> Hello, >> I _might_ be doing something wrong

Re: [petsc-dev] Right-preconditioned GMRES

> On Feb 9, 2017, at 1:37 PM, Matthew Knepley wrote: > > On Thu, Feb 9, 2017 at 4:56 AM, Pierre Jolivet > wrote: > Hello, > I _might_ be doing something wrong when calling KSPSetFromOptions/KSPSetUp, > but

Re: [petsc-dev] Right-preconditioned GMRES

On Thu, Feb 9, 2017 at 4:56 AM, Pierre Jolivet wrote: > Hello, > I _might_ be doing something wrong when calling > KSPSetFromOptions/KSPSetUp, but still, am I supposed to get this kind of > error? > [0]PETSC ERROR: No support for this operation for this object type >

[petsc-dev] Right-preconditioned GMRES

Hello, I _might_ be doing something wrong when calling KSPSetFromOptions/KSPSetUp, but still, am I supposed to get this kind of error? [0]PETSC ERROR: No support for this operation for this object type [0]PETSC ERROR: KSP gmres does not support PRECONDITIONED with RIGHT If not, I’ll try to send