[petsc-users] Problem Iterating DMPlex

2015-11-17 Thread Morten Nobel-Jørgensen
After distributing a DMPlex it seems like my cells are appearing twice (or rather multiple cells maps onto the same vertices). I’m assuming the way I’m iterating the DMPlex is wrong. Essentially I iterate the DMPlex the following way after distribution (see code snippet below – or attached

[petsc-users] Understanding PETSC: boundary layer flow with SNES

2015-11-17 Thread Pavel Schor
Dear PETSC users, I am a newbie to PETSC and I try to understand it. My first attempt is to solve 2D laminar boundary layer equations: u*du/dx +v*du/dy = Ue*(dUe/dx). -nu*d2u/dy2 dudx + dvdy = 0 B.C.: y=0:u=0,v=0 y=H:u(x)=Ue(x) Where H is height of the domain,

Re: [petsc-users] problem using 64-bit-indices

2015-11-17 Thread Satish Balay
ierr is of type 'PetscErrorCode' - not 'PetscInt' Perhaps the dmdavecgetarrayf90 issue is also due to the difference in types used for parameters Satish On Tue, 17 Nov 2015, Randall Mackie wrote: > I ran into a problem yesterday where a call to DMDACreate3d gave an error > message about the

[petsc-users] problem using 64-bit-indices

2015-11-17 Thread Randall Mackie
I ran into a problem yesterday where a call to DMDACreate3d gave an error message about the size being too big and that I should use —with-64-bit-indices. So I recompiled PETSc (latest version, 3.6.2) with that option, but when I recompiled my code, I found the following errors: call

Re: [petsc-users] problem using 64-bit-indices

2015-11-17 Thread Randall Mackie
Thanks Satish, This fixed all error messages. Randy > On Nov 17, 2015, at 9:12 AM, Satish Balay wrote: > > ierr is of type 'PetscErrorCode' - not 'PetscInt' > > Perhaps the dmdavecgetarrayf90 issue is also due to the difference in types > used for parameters > > Satish