On Thu, 2011-03-10 at 23:53 -0300, Lisandro Dalcin wrote: > While updating petsc4py for the latest DMDA chages, and alto taking > into account some new in NumPy C API related to neighborhood > iterators, I've got new idea. What about the API below?: > > enum DMDABoundaryType: > DMDA_BOUNDARY_NONE = 0 > DMDA_BOUNDARY_GHOSTED = 1 > DMDA_BOUNDARY_MIRROR = 2 > DMDA_BOUNDARY_PERIODIC = 4 > > PetscErrorCode DMDASetBoundaryType(DM dm, DMDABoundaryType bx, > DMDABoundaryType by, DMDABoundaryType bz) > > Changes for GHOSTED and PERIODIC should be easy > (http://lkml.org/lkml/2000/8/25/132) ;-). > > The new MIRROR boundary type could be implemented later. >
How does this fit into DMDACreate*D()? Are you changing the create API as well? PetscErrorCode DMDACreate3d(MPI_Comm comm,DMDABoundaryType bx, DMDABoundaryType by, DMDABoundaryType bz, DMDAStencilType stencil_type,PetscInt M,PetscInt N,PetscInt P,PetscInt m,PetscInt n,PetscInt p,PetscInt dof,PetscInt s,const PetscInt lx[],const PetscInt ly[],const PetscInt lz[],DM *da) Ethan > > -- ------------------------------------ Ethan Coon Post-Doctoral Researcher Applied Mathematics - T-5 Los Alamos National Laboratory 505-665-8289 http://www.ldeo.columbia.edu/~ecoon/ ------------------------------------
