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 >> -pc

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 works

Re: [petsc-dev] Telescope usage

2022-05-17 Thread 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 > -pc_mg_levels 2 -ksp_view -mg_coarse_pc_type telescope

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 -da_gr

Re: [petsc-dev] Telescope usage

2022-05-17 Thread Loic Gouarin
Le 17/05/2022 à 12:03, Dave May a écrit : 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 defi

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 à 01

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 9:5

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: Th

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Dave May
On Tue, 17 May 2022 at 00:12, Matthew Knepley wrote: > On Mon, May 16, 2022 at 9:59 AM Loic Gouarin < > loic.goua...@polytechnique.edu> wrote: > >> Le 16/05/2022 à 21:54, Barry Smith a écrit : >> >> >> >> On May 16, 2022, at 3:50 PM, Loic Gouarin >> wrote: >> >> Thanks Barry for your quick reply

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Matthew Knepley
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: > > Thanks Barry for your quick reply. > Le 16/05/2022 à 21:41, Barry Smith a écrit : > > > Loic, > > From your code it looks like y

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Loic Gouarin
Le 16/05/2022 à 21:54, Barry Smith a écrit : On May 16, 2022, at 3:50 PM, Loic Gouarin wrote: Thanks Barry for your quick reply. Le 16/05/2022 à 21:41, Barry Smith a écrit :   Loic,     From your code it looks like you are using a DM. Is it a DMDA or a shell DM? If it is a DMDA then t

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Barry Smith
> On May 16, 2022, at 3:50 PM, Loic Gouarin > wrote: > > Thanks Barry for your quick reply. > > Le 16/05/2022 à 21:41, Barry Smith a écrit : >> >> Loic, >> >> From your code it looks like you are using a DM. Is it a DMDA or a shell >> DM? If it is a DMDA then the process is intended

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Loic Gouarin
Thanks Barry for your quick reply. Le 16/05/2022 à 21:41, Barry Smith a écrit :   Loic,     From your code it looks like you are using a DM. Is it a DMDA or a shell DM? If it is a DMDA then the process is intended to be pretty straightforward. PCTELESCOPE should create a new DMDA that has th

Re: [petsc-dev] Telescope usage

2022-05-16 Thread Barry Smith
Loic, From your code it looks like you are using a DM. Is it a DMDA or a shell DM? If it is a DMDA then the process is intended to be pretty straightforward. PCTELESCOPE should create a new DMDA that has the same properties as the coarse grid DMDA but lives on a smaller number of ranks.

[petsc-dev] Telescope usage

2022-05-16 Thread Loic Gouarin
Hello, I could have posted my message on the user list but it seems to me that it's more in the petsc development part. Don't hesitate to tell me if I'm wrong. I am developing a code called cafes that studies fluid-particle interactions in a Stokes fluid. For that, we implemented the whole so