Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

2022-07-20 Thread Bora Jeong
Thank you for the comments. But If the ordering before DMPlexDistribute is what you're looking for, the SF directly tells you the mapping. The leaves of the natural SF gives you the indices from your file. You just need to shift it by +1. *-> I do not agree. Below is my Natural SF with 2 proc,

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

2022-07-20 Thread Bora Jeong
Thank you for the comments. If you define a scalar P1 field of one component, then the dof index will be the vertex index. *-> When you say P1 field of one component, what does exactly mean? I put 1dof scalar field on vertex when the Petsc section is defined. Is this P0 or P1 in Petsc's

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

2022-07-20 Thread Alexis Marboeuf
Hi Bora, I agree with Matt.  I don't know what DMPlexCreateGmshFromFile() is doing in term of distribution then. But If the ordering before DMPlexDistribute is what you're looking for, the SF directly tells you the mapping. The leaves of the natural SF gives you the indices from your file.

Re: [petsc-users] Reusing A00 block with Schur complement

2022-07-20 Thread Barry Smith
It is suppose to not rebuild the preconditioner for the A00 block if the matrix does not change; and automatically rebuild if it does change. > On Jul 20, 2022, at 4:14 PM, Mike Welland wrote: > > My A00 block in an SNES solve with a Schur complement preconditioner is > constant, and I'd

[petsc-users] Reusing A00 block with Schur complement

2022-07-20 Thread Mike Welland
My A00 block in an SNES solve with a Schur complement preconditioner is constant, and I'd like to exploit this for computational speedup. I'm running PETSc through a FEM package (Firedrake) so only have ready access to the options database (I believe that to be the command line arguments) Is this

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

2022-07-20 Thread Matthew Knepley
On Wed, Jul 20, 2022 at 2:27 PM Bora Jeong wrote: > I am not seeing DMPlexDistributeSetDefault(dm,PETSC_FALSE) after you > create your DM from the GMSH file. Your DM is then automatically > distributed and the natural ordering won't be a priori the one from your > file. I don't know how Petsc

Re: [petsc-users] DMPlex: a Mapping Array between Natural and Distributed Cell Index

2022-07-20 Thread Bora Jeong
I am not seeing DMPlexDistributeSetDefault(dm,PETSC_FALSE) after you create your DM from the GMSH file. Your DM is then automatically distributed and the natural ordering won't be a priori the one from your file. I don't know how Petsc distributes a DM during the load step from a file. *->

Re: [petsc-users] Possible to turn off unused option warning in debug mode?

2022-07-20 Thread Barry Smith
You can put these in a .petscrc file in your home directory, the current directory or in the environmental variable PETSC_OPTIONS (I do this) > On Jul 20, 2022, at 2:43 PM, Jacob Faibussowitsch wrote: > > Run with "-options_left 0" (also accepts no, and false) > > Best regards, > >

Re: [petsc-users] Possible to turn off unused option warning in debug mode?

2022-07-20 Thread Jacob Faibussowitsch
Run with "-options_left 0" (also accepts no, and false) Best regards, Jacob Faibussowitsch (Jacob Fai - booss - oh - vitch) > On Jul 20, 2022, at 14:31, Greg Kahanamoku-Meyer > wrote: > > Hi, > > The unused option warning >

[petsc-users] Possible to turn off unused option warning in debug mode?

2022-07-20 Thread Greg Kahanamoku-Meyer
Hi, The unused option warning ( https://petsc.org/release/docs/manualpages/Sys/PetscOptionsLeft.html) is automatically turned on when running in debug mode. Is there a way to turn it off (other than turning off all debugging)? For a bit of background, I am using PETSc through petsc4py, and