Re: [petsc-users] strong-scaling vs weak-scaling

2016-08-30 Thread Jed Brown
Mark Adams writes: >> >> >> Anyway, what I really wanted to say is, it's good to know that these >> "dynamic range/performance spectrum/static scaling" plots are designed to >> go past the sweet spots. I also agree that it would be interesting to see a >> time vs

Re: [petsc-users] strong-scaling vs weak-scaling

2016-08-30 Thread Mark Adams
> > > Anyway, what I really wanted to say is, it's good to know that these > "dynamic range/performance spectrum/static scaling" plots are designed to > go past the sweet spots. I also agree that it would be interesting to see a > time vs dofs*iterations/time plot. Would it then also be useful to

Re: [petsc-users] strong-scaling vs weak-scaling

2016-08-30 Thread Justin Chang
Thanks everyone. I still think there is a even better phrase for this, like, static scaling? Because unlike strong/weak scaling, concurrency is fixed (hence "static") and we only scale the problem, so this is a mix between strong and weak scaling. ¯\_(ツ)_/¯ Anyway, what I really wanted to say

Re: [petsc-users] Distribution of DMPlex for FEM

2016-08-30 Thread Morten Nobel-Jørgensen
We have now hit a related problem. If we change the dof from 1 to 3 we the following error message. I'm using the next branch (pulled from git today). mpiexec -np 2 ./ex18k [1]PETSC ERROR: - Error Message --

Re: [petsc-users] strong-scaling vs weak-scaling

2016-08-30 Thread Jed Brown
Mark Adams writes: > I would guess it is the latter. In this case, definitely. > It is hard to get "rollover" to the right. You could get it on KNL > (cache configuration of HBM) when you spill out of HBM. Yes, but the same occurs if you start repeatedly spilling from some

Re: [petsc-users] strong-scaling vs weak-scaling

2016-08-30 Thread Mark Adams
On Tue, Aug 23, 2016 at 4:33 AM, Justin Chang wrote: > Redid some of those experiments for 8 and 20 cores and scaled it up to > even larger problems. Attached is the plot. > > Looking at this "dynamic plot" (if you ask me, I honestly think there > could be a better word for

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Florian Lindner
Thanks everybody, just to be sure, it's max/min of the last iteration? Florian Am 30.08.2016 um 16:10 schrieb Matthew Knepley: > On Tue, Aug 30, 2016 at 9:03 AM, Florian Lindner > wrote: > > Hi, > > Am 30.08.2016 um 14:05 schrieb

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Hong
max/min 2.829815526858e+06 is an estimate to your condition number. Hong On Tue, Aug 30, 2016 at 9:03 AM, Florian Lindner wrote: > Hi, > > Am 30.08.2016 um 14:05 schrieb Barry Smith: > > > >The format of .petscrc requires each option to be on its own line > > > >>

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Song Gao
condition number is max_eigenvalue/min_eigenvalue. so I guess it is max/min 2.829815526858e+06? 2016-08-30 10:03 GMT-04:00 Florian Lindner : > Hi, > > Am 30.08.2016 um 14:05 schrieb Barry Smith: > > > >The format of .petscrc requires each option to be on its own line >

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Matthew Knepley
On Tue, Aug 30, 2016 at 9:03 AM, Florian Lindner wrote: > Hi, > > Am 30.08.2016 um 14:05 schrieb Barry Smith: > > > >The format of .petscrc requires each option to be on its own line > > > >> -ksp_view > >> -pc_type none > >> -ksp_type gmres > >>

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Florian Lindner
Hi, Am 30.08.2016 um 14:05 schrieb Barry Smith: > >The format of .petscrc requires each option to be on its own line > >> -ksp_view >> -pc_type none >> -ksp_type gmres >> -ksp_monitor_singular_value >> -ksp_gmres_restart 1000 Oh man, didn't know that. Sorry! Is using a hash # ok for

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Matthew Knepley
On Tue, Aug 30, 2016 at 7:01 AM, Florian Lindner wrote: > Hello, > > there is a FAQ and a Stackoverflow article about getting the condition > number of a petsc matrix: > > http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber >

Re: [petsc-users] Condition number of matrix

2016-08-30 Thread Barry Smith
The format of .petscrc requires each option to be on its own line > -ksp_view > -pc_type none > -ksp_type gmres > -ksp_monitor_singular_value > -ksp_gmres_restart 1000 > On Aug 30, 2016, at 7:01 AM, Florian Lindner wrote: > > Hello, > > there is a FAQ and a

[petsc-users] Condition number of matrix

2016-08-30 Thread Florian Lindner
Hello, there is a FAQ and a Stackoverflow article about getting the condition number of a petsc matrix: http://www.mcs.anl.gov/petsc/documentation/faq.html#conditionnumber http://scicomp.stackexchange.com/questions/34/how-can-i-estimate-the-condition-number-of-a-large-sparse-matrix-using-petsc