Re: [petsc-users] ex12 poisson solver

2017-10-10 Thread Matthew Knepley
On Tue, Oct 10, 2017 at 5:09 AM, Barry Smith wrote: > > > On Oct 9, 2017, at 6:56 PM, Matthew Knepley wrote: > > > > On Mon, Oct 9, 2017 at 12:40 PM, Justin Chang > wrote: > > You need this additional command line argument: -petscspace_order 1 > > > > Justin is correct. By default, it is P0, wh

Re: [petsc-users] ex12 poisson solver

2017-10-10 Thread Barry Smith
> On Oct 9, 2017, at 6:56 PM, Matthew Knepley wrote: > > On Mon, Oct 9, 2017 at 12:40 PM, Justin Chang wrote: > You need this additional command line argument: -petscspace_order 1 > > Justin is correct. By default, it is P0, which would be a purely > discontinuous solution to Laplace, without

Re: [petsc-users] ex12 poisson solver

2017-10-09 Thread Matthew Knepley
On Mon, Oct 9, 2017 at 12:56 PM, Matthew Knepley wrote: > On Mon, Oct 9, 2017 at 12:40 PM, Justin Chang wrote: > >> You need this additional command line argument: -petscspace_order 1 >> > > Justin is correct. By default, it is P0, which would be a purely > discontinuous solution to Laplace, wit

Re: [petsc-users] ex12 poisson solver

2017-10-09 Thread Matthew Knepley
On Mon, Oct 9, 2017 at 12:40 PM, Justin Chang wrote: > You need this additional command line argument: -petscspace_order 1 > Justin is correct. By default, it is P0, which would be a purely discontinuous solution to Laplace, without jump terms. Thanks, Matt > On Mon, Oct 9, 2017 at 9:0

Re: [petsc-users] ex12 poisson solver

2017-10-09 Thread Justin Chang
You need this additional command line argument: -petscspace_order 1 On Mon, Oct 9, 2017 at 9:02 AM, David Fuentes wrote: > > Hi, > > I'm trying to use petsc 3.8.0 with ex12.c example to setup a poisson > solver: http://www.mcs.anl.gov/petsc/petsc-current/src/snes/ > examples/tutorials/ex12.c.h

[petsc-users] ex12 poisson solver

2017-10-09 Thread David Fuentes
Hi, I'm trying to use petsc 3.8.0 with ex12.c example to setup a poisson solver: http://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex12.c.html I seem to be getting zeros in my jacobian with this example? I attached a debugger and the assembly routines seems ok... but am getti