Re: [petsc-users] non-manifold DMPLEX

2021-12-12 Thread Patrick Sanan
This would be a particularly useful example for testing the PETSc library, as most applications are set on manifolds. There’s a recently-reformatted chapter in the developers docs on how the testing works - this explains the /* TEST */ blocks you see. If parts of these docs are confusing, that’s ve

Re: [petsc-users] SNES always ends at iteration 0

2021-12-12 Thread Mark Adams
On Sun, Dec 12, 2021 at 7:47 PM celestechev...@gmail.com < celestechev...@gmail.com> wrote: > Thank you so much for your reply ! > > In fact I didn't know how to set tolerances, so I proceeded without > specifying the tolerances, hoping that this could lead to the > implementation of default value

Re: [petsc-users] SNES always ends at iteration 0

2021-12-12 Thread Matthew Knepley
On Sun, Dec 12, 2021 at 7:48 PM celestechev...@gmail.com < celestechev...@gmail.com> wrote: > Thank you so much for your reply ! > > In fact I didn't know how to set tolerances, so I proceeded without > specifying the tolerances, hoping that this could lead to the > implementation of default value

Re: [petsc-users] non-manifold DMPLEX

2021-12-12 Thread Matthew Knepley
On Sun, Dec 12, 2021 at 4:36 PM TARDIEU Nicolas wrote: > Dear Patrick and Matthew, > > Thank you very much for your answers. I am gonna try to set up such a test > by assigning cell types. > Shall I open a MR in order to contribute this example ? > Yes, that would be great. Thanks, Matt

Re: [petsc-users] SNES always ends at iteration 0

2021-12-12 Thread celestechev...@gmail.com
Thank you so much for your reply ! In fact I didn't know how to set tolerances, so I proceeded without specifying the tolerances, hoping that this could lead to the implementation of default values... I just added *SNESSetFromOptions(snes); *However, it doesn't make any difference... But it's tru

Re: [petsc-users] SNES always ends at iteration 0

2021-12-12 Thread Mark Adams
The three "tol" values should be finite. It sounds like you set them to 0. Don't do that and use the defaults to start. The behavior with zero tolerances is not defined. You can use -snes_monitor to print out the iterations. On Sun, Dec 12, 2021 at 6:22 PM celestechev...@gmail.com < celestechev...

[petsc-users] SNES always ends at iteration 0

2021-12-12 Thread celestechev...@gmail.com
Hello, I encountered a strange problem concerning the convergence of SNES. In my recent test runs I found that SNES always stops at iteration 0... At first I thought there may be an error with the tolerance setting, so I output the tolerances : *atol = 0.00, rtol = 0.00, stol = 0.

Re: [petsc-users] Unstructured mesh

2021-12-12 Thread Mark Adams
On Sun, Dec 12, 2021 at 3:19 PM Karthikeyan Chockalingam - STFC UKRI < karthikeyan.chockalin...@stfc.ac.uk> wrote: > Thank for your response that was helpful. I have a couple of questions: > > > > (i) How can I control the level of refinement? I tried > to pass the flag “-ex56_dm_

Re: [petsc-users] non-manifold DMPLEX

2021-12-12 Thread TARDIEU Nicolas via petsc-users
Dear Patrick and Matthew, Thank you very much for your answers. I am gonna try to set up such a test by assigning cell types. Shall I open a MR in order to contribute this example ? Regards, Nicolas De : knep...@gmail.com Envoyé : dimanche 12 décembre 2021 12:1

Re: [petsc-users] Unstructured mesh

2021-12-12 Thread Karthikeyan Chockalingam - STFC UKRI
Thank for your response that was helpful. I have a couple of questions: (i) How can I control the level of refinement? I tried to pass the flag “-ex56_dm_refine 0” but that didn’t stop the refinement from 8 giving 32 cubes. (ii)What does -cell 2,2,1 correspond

Re: [petsc-users] non-manifold DMPLEX

2021-12-12 Thread Matthew Knepley
On Sun, Dec 12, 2021 at 6:11 AM Patrick Sanan wrote: > Here you have the following "points": > > - 1 3-cell (the cube volume) > - 7 2-cells (the 6 faces of the cube plus the extra one) > - 16 1-cells (the 12 edges of the cube, plus 3 extra ones from the extra > face, plus the extra edge) > - 11

Re: [petsc-users] non-manifold DMPLEX

2021-12-12 Thread Patrick Sanan
Here you have the following "points": - 1 3-cell (the cube volume) - 7 2-cells (the 6 faces of the cube plus the extra one) - 16 1-cells (the 12 edges of the cube, plus 3 extra ones from the extra face, plus the extra edge) - 11 0-cells (the 8 vertices of the cube, pus 2 extra ones from the extra