Hey Matt, Yes I have tried all possible variants: [(8,8,4),] ([8,8,4],) [[8,8,4],] ((8,8,4),) and they all fail with the same TypeError: Expected tuple, got list However, I think I see what is happening now. In line 156 cdef object ranges = list(ownership_ranges) it looks like ranges is being defined as a list, rather than a tuple. So maybe there should be a cast from list to tuple in the return statement? -Chris On Thu, Apr 7, 2016 at 9:19 PM, Matthew Knepley <[email protected]> wrote:

On Thu, Apr 7, 2016 at 2:14 PM, Chris Eldred <[email protected]> wrote:

...

I think what happens on the 1st version is that the input list has len(3), and since the DMDA only has 1 dimension petsc4py throws the first ValueError.

Yes, you are right.

...

As far as I can tell, the code is expecting a list of lists: ie for each dimension it gets ranges[i] (where i is the index for the dimension) and casts it into an array. So I tried the 2nd version, which fails regardless of whether I use lists or tuples, and I don't understand why.

That function returns 'range' and the type is declared as tuple. Did you try

ownership_ranges=((8,8,4),)

Matt

...

-Chris

On Thu, Apr 7, 2016 at 9:00 PM, Matthew Knepley <[email protected]> wrote:

...

On Thu, Apr 7, 2016 at 1:52 PM, Chris Eldred <[email protected]> wrote:

...

Hey All,

(I sent this to petsc-users a few days ago and didn't get any replies, figured I would try here as well...)

Does anyone know the correct syntax for setting ownership ranges for DMDAs in petsc4py? For example, if I have a 1D DMDA with 20 vertices and 3 processes, and I want to distribute the vertices as [8,8,4], should the correct create call be:

da =

PETSc.DMDA().create(dim=1,dof=1,sizes=[20,],proc_sizes=[3,],boundary_type=[PETSc.DM.BoundaryType.PERIODIC,], stencil_type= PETSc.DMDA.StencilType.BOX, stencil_width=1,ownership_ranges=[8,8,4])

or

da =

PETSc.DMDA().create(dim=1,dof=1,sizes=[20,],proc_sizes=[3,],boundary_type=[PETSc.DM.BoundaryType.PERIODIC,], stencil_type= PETSc.DMDA.StencilType.BOX, stencil_width=1,ownership_ranges=[[8,8,4],])

For me, the 1st version fails with ValueError: number of dimensions 1 and number ownership ranges 3 (which I expected) but the 2nd version also fails with TypeError: Expected tuple, got list

I have tried using tuples instead of lists and nothing changes. Creating without the ownership_ranges argument works fine as well. Is there something else I am missing?

Here is the code:

https://bitbucket.org/petsc/petsc4py/src/b234a50085607b793fce93d5dc46307c0d9...

I do not understand why the 1st version fails.

Matt

...

Thanks in advance for any assistance!

Regards, -Chris

-- Chris Eldred Postdoctoral Fellow, LAGA, University of Paris 13 PhD, Atmospheric Science, Colorado State University, 2015 DOE Computational Science Graduate Fellow (Alumni) B.S. Applied Computational Physics, Carnegie Mellon University, 2009 [email protected]

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

-- Chris Eldred Postdoctoral Fellow, LAGA, University of Paris 13 PhD, Atmospheric Science, Colorado State University, 2015 DOE Computational Science Graduate Fellow (Alumni) B.S. Applied Computational Physics, Carnegie Mellon University, 2009 [email protected]

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener

-- Chris Eldred Postdoctoral Fellow, LAGA, University of Paris 13 PhD, Atmospheric Science, Colorado State University, 2015 DOE Computational Science Graduate Fellow (Alumni) B.S. Applied Computational Physics, Carnegie Mellon University, 2009 [email protected]