Re: [petsc-dev] Right-preconditioned GMRES

2019-11-06 Thread Jed Brown via petsc-dev
"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 > barrysmith@barry-smiths-computer.l

Re: [petsc-dev] Right-preconditioned GMRES

2019-11-06 Thread Smith, Barry F. via petsc-dev
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

2019-11-06 Thread Pierre Jolivet via petsc-dev
> 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 preconditi

Re: [petsc-dev] Right-preconditioned GMRES

2019-11-06 Thread Pierre Jolivet via petsc-dev
> 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

2019-11-06 Thread Pierre Jolivet via petsc-dev
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 al

Re: [petsc-dev] GPU counters

2019-11-06 Thread Mark Adams via petsc-dev
Yea, that is what I thought. Oh, I am probably seeing the flops from KSP. The PC is a monolithic code (AMGx). On Wed, Nov 6, 2019 at 11:18 AM Zhang, Junchao wrote: > No. For each vector/matrix operation, PETSc can get its flop count based > on number of nonzeros, for example. > > --Junchao Zhang

Re: [petsc-dev] GPU counters

2019-11-06 Thread Zhang, Junchao via petsc-dev
No. For each vector/matrix operation, PETSc can get its flop count based on number of nonzeros, for example. --Junchao Zhang On Wed, Nov 6, 2019 at 8:44 AM Mark Adams via petsc-dev mailto:petsc-dev@mcs.anl.gov>> wrote: I am puzzled. I am running AMGx now, and I am getting flop counts/rates. H

[petsc-dev] GPU counters

2019-11-06 Thread Mark Adams via petsc-dev
I am puzzled. I am running AMGx now, and I am getting flop counts/rates. How does that happen? Does PETSc use hardware counters to get flops?