<div dir="ltr">On Wed, Dec 26, 2012 at 12:05 PM, Zhenglun (Alan) Wei <span dir="ltr"><<a href="mailto:z240w014@ku.edu" target="_blank">z240w014@ku.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  

    
  
  <div text="#000000" bgcolor="#FFFFFF">
    Dear folks,<br>
         I have a quick question on the DMDACreate3d.<br>
         In the manual, it says that the input format of this function
    is:<br>
    <pre style="line-height:normal;text-indent:0px;letter-spacing:normal;text-align:-webkit-auto;font-variant:normal;text-transform:none;font-style:normal;font-weight:normal;word-spacing:0px">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)
</pre>
    <br>
         Now, I'm trying to manually define the <span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;display:inline!important;font-weight:normal;float:none;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:'Times New Roman';word-spacing:0px">"arrays containing the number of nodes in each cell along the x,
      y, and z coordinates</span>". Therefore, my focus turns to 'lx[]',
    'ly[]' and 'lz[]'. I suppose that they're not simply just three
    integers; they may be three integer type arrays, as I guess.
    However, I checked all examples listed for this function. None of
    them teaches me how to implement this three parameters except
    'PETSC_NULL'. Could you please provide me an extra example to
    demonstrate how to use DMDACreate3d or DMDACreate2d with non-null
    'lx[]', 'ly[]' and 'lz[]'.<br></div></blockquote><div><br></div><div style>It is used by snes/examples/tutorials/ex28.c in the 1D case to ensure that the staggered grid has a specific compatible layout. As the docs say, these are arrays of length m,n,p and must sum to M, N, and P.</div>
<div style><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF">
         Or, a demonstration in 1D would be a good example. Say, I have
    a 1D uniform mesh; the number of grid in x-direction is 300. I want
    to use 4 processes to evenly divide this mesh. What should I input
    for 'lx[]' for each process?<br>
    <br></div></blockquote><div><br></div><div style>The defaults do this.</div></div></div></div>