Hi Barry,
PetscSplitOwnership works indeed even though it needs some further work to
determine the local ownership istart and iend (only gives the ownership size)
I found out you can "cheat" and call
MatMPIAIJSetallocation(Mat,0,PETSC_NULL_INTEGER,0,PETSC_NULL_INTEGER), then
MatGetOwnershipRa
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
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
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
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