[petsc-users] DMDA with dimension of size 1

[Phil Tooley]

Hi all,

Is it valid to have a DMDA with one of the dimensions of size 1.  I was
hoping to avoid having to write explicit logic to handle the case that I
am working on a 2D rather than a 3D image for my current application.  
Unfortunately when I try to construct such a DMDA I get an error:

[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: Overflow in integer operation:
[0]PETSC ERROR: Mesh of 10 by 10 by 1 (dof) is too large for 32 bit indices

Is there a workaround other than "write everything twice"?

Thanks

Phil

-- 
Phil Tooley
Research Software Engineering
University of Sheffield