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
>
> On Tue, 17 Nov
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 s
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 VecGetA
On Tue, Nov 17, 2015 at 10:11 AM, Morten Nobel-Jørgensen
wrote:
> After distributing a DMPlex it seems like my cells are appearing twice (or
> rather multiple cells maps onto the same vertices).
>
Your code is creating the serial mesh on each process. You only want
nonzero sizes on one proc.
Yo
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 file
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, Ue(x