[petsc-users] Reference for multi-grid tuning

2016-08-03 Thread Justin Chang
Hi all, I want to play around with algebraic multi-grid parameters to optimize performance for my 3D problems but I don't have a good idea of where to start. Don't really see a detailed reference or tutorial for tuning these parameters in the PETSc manual. Can someone point me to a good paper/ref

Re: [petsc-users] Multi-physics meshes with PETSc DM?

2016-08-03 Thread Matthew Knepley
On Wed, Aug 3, 2016 at 4:57 PM, Andrew Ho wrote: > It will not be periodic since periodic meshes in Plex are topologically >> periodic, meaning that >> there is no separate aliasing array. I could possibly read in the GMsh >> thing and do surgery >> on the mesh. >> > > I'm not sure what you mean

Re: [petsc-users] Multi-physics meshes with PETSc DM?

2016-08-03 Thread Andrew Ho
> > It will not be periodic since periodic meshes in Plex are topologically > periodic, meaning that > there is no separate aliasing array. I could possibly read in the GMsh > thing and do surgery > on the mesh. > I'm not sure what you mean by this, but being able to handle periodic domains is imp

Re: [petsc-users] turn off the snes linesearch

2016-08-03 Thread Matthew Knepley
On Wed, Aug 3, 2016 at 4:39 PM, Xiangdong wrote: > Hello everyone, > > If I want to turn off the line search and always use the full newton step, > does the option "-snes_linesearch_type basic" guarantee this? In my case, I > found that even with that option, it still try to do function evaluatio

[petsc-users] turn off the snes linesearch

2016-08-03 Thread Xiangdong
Hello everyone, If I want to turn off the line search and always use the full newton step, does the option "-snes_linesearch_type basic" guarantee this? In my case, I found that even with that option, it still try to do function evaluation at a point different from the exact newton step. Any sugg

Re: [petsc-users] Multi-physics meshes with PETSc DM?

2016-08-03 Thread Matthew Knepley
On Sat, Jul 30, 2016 at 1:35 PM, Andrew Ho wrote: > Is there a reason the physical groups aren't sufficient for handling this? > As far as I can tell, this is the only way in GMsh to have any kind of > grouping of elements. > > The Gmsh file format can be found here (happens to be the ASCII versi

Re: [petsc-users] false-positive leak report in log_view?

2016-08-03 Thread Mohammad Mirzadeh
OK so I just ran the example under valgrind, and if I use two VecViews, it complains about following leak: ==66838== 24,802 (544 direct, 24,258 indirect) bytes in 1 blocks are definitely lost in loss record 924 of 926 ==66838==at 0x19EBB: malloc (in /usr/local/Cellar/valgrind/3.11.0/lib/va

Re: [petsc-users] false-positive leak report in log_view?

2016-08-03 Thread Mohammad Mirzadeh
On Wed, Aug 3, 2016 at 10:59 AM, Matthew Knepley wrote: > On Tue, Aug 2, 2016 at 12:40 PM, Mohammad Mirzadeh > wrote: > >> I often use the memory usage information in log_view as a way to check on >> memory leaks and so far it has worked perfect. However, I had long noticed >> a false-positive r

Re: [petsc-users] Meaning of error message (gamg & fieldsplit related)

2016-08-03 Thread Barry Smith
This should be fixed in the master branch > On Aug 1, 2016, at 11:13 PM, Safin, Artur wrote: > > Hello, > > I am running some code that employs gamg preconditioning within a fieldsplit, > and for sufficiently large/refined meshes, I am getting the following error: > > ---

Re: [petsc-users] Meaning of error message (gamg & fieldsplit related)

2016-08-03 Thread Barry Smith
This is ixed in the master branch of PETSc. Barry > On Aug 1, 2016, at 11:13 PM, Safin, Artur wrote: > > Hello, > > I am running some code that employs gamg preconditioning within a fieldsplit, > and for sufficiently large/refined meshes, I am getting the following error: > >

Re: [petsc-users] Implementing discontinuous Galerkin FEM?

2016-08-03 Thread Matthew Knepley
On Thu, Jul 28, 2016 at 1:31 PM, Andrew Ho wrote: > I am trying to implement a discontinuous Galerkin discretization using the > PETSc DM features to handle most of the topology/geometry specific > functions. However, I'm not really sure which direction to approach this > from since DG is kind of

Re: [petsc-users] false-positive leak report in log_view?

2016-08-03 Thread Matthew Knepley
On Tue, Aug 2, 2016 at 12:40 PM, Mohammad Mirzadeh wrote: > I often use the memory usage information in log_view as a way to check on > memory leaks and so far it has worked perfect. However, I had long noticed > a false-positive report in memory leak for Viewers, i.e. destruction count > is alwa

Re: [petsc-users] Meaning of error message (gamg & fieldsplit related)

2016-08-03 Thread Mark Adams
Are you saying it works for a while but fails when the problem is large, or that it never works with fieldsplit_1? And how many processors are you using? On Mon, Aug 1, 2016 at 9:13 PM, Safin, Artur wrote: > Hello, > > > I am running some code that employs gamg preconditioning within a > fieldsp

Re: [petsc-users] TS and petscFE

2016-08-03 Thread Matthew Knepley
On Tue, Aug 2, 2016 at 8:22 AM, Maximilian Hartig wrote: > Hello all, > > I would like to run a transient problem with PetscFE. Example ex11.c seems > relevant since it uses the PestcFV context to create boundary conditions > and RHS Functions for the TS. > Is there an easy way to do transient an