Re: [petsc-dev] Telescope usage

2022-05-17 Thread Loïc Gouarin
17 mai 2022 17:53:46 Matthew Knepley : > On Tue, May 17, 2022 at 4:53 AM Dave May wrote: >> Dear Loic, >> >> I can confirm that PCTELESCOPE works nicely when using >> KSPSetComputeOperators. >> Here is an example >> >> petsc-3.17.1/src/ksp/ksp/tutorials $ mpiexec -n 4 ./ex29 -pc_type mg >>

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Loïc Gouarin
Thanks Dave for your time and this detailed reply. My implementation is old and should probably be updated. I will read your email carefully and I come back to you with a better implementation. Thanks, Loic 17 mai 2022 16:53:37 Dave May : > Dear Loic, > > I can confirm that PCTELESCOPE

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Dave May
Dear Loic, I can confirm that PCTELESCOPE works nicely when using KSPSetComputeOperators. Here is an example petsc-3.17.1/src/ksp/ksp/tutorials $ mpiexec -n 4 ./ex29 -pc_type mg -pc_mg_levels 2 -ksp_view -mg_coarse_pc_type telescope -mg_coarse_pc_telescope_reduction_factor 4 -da_grid_x 65

Re: [petsc-dev] odd log behavior

2022-05-17 Thread Zhang, Hong via petsc-dev
Python users including myself would love NaN since NaN is the default missing value marker for reasons of computational speed and convenience. For example, if you load these values into pandas, no extra code is needed to handle them. Other choices such as N/A would require some extra work for

Re: [petsc-dev] odd log behavior

2022-05-17 Thread Mark Adams
Oh, this is not hard: /* /\* Put NaN into the time for all events that may not be time accurately since they may happen asynchronously on the GPU *\/ */ /* #if defined(PETSC_HAVE_DEVICE) */ /* if (!PetscLogGpuTimeFlag && petsc_gflops > 0) { */ /*

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Dave May
On Tue 17. May 2022 at 11:56, Loic Gouarin wrote: > Hi Dave, > > could you explain what you mean by state ? > Ah - by state Matt and I mean any auxiliary data stored in the user context passed to KSPSetComputeOperators which is required to define your operator and is distributed. For example, a

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Loic Gouarin
Hi Dave, could you explain what you mean by state ? Thanks, Loic Le 17/05/2022 à 11:50, Dave May a écrit : Hi Loic, Can you confirm if your problem has stored state which needs to be repartitioned? Thanks, Dave On Tue 17. May 2022 at 09:07, Loic Gouarin wrote: Le 17/05/2022 à

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Dave May
Hi Loic, Can you confirm if your problem has stored state which needs to be repartitioned? Thanks, Dave On Tue 17. May 2022 at 09:07, Loic Gouarin wrote: > > Le 17/05/2022 à 01:43, Dave May a écrit : > > > > On Tue, 17 May 2022 at 00:12, Matthew Knepley wrote: > >> On Mon, May 16, 2022 at

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Loic Gouarin
Le 17/05/2022 à 01:43, Dave May a écrit : On Tue, 17 May 2022 at 00:12, Matthew Knepley wrote: On Mon, May 16, 2022 at 9:59 AM Loic Gouarin wrote: Le 16/05/2022 à 21:54, Barry Smith a écrit : On May 16, 2022, at 3:50 PM, Loic Gouarin wrote: