<div dir="ltr">On Tue, Jan 22, 2013 at 3:34 AM, Hui Zhang <span dir="ltr"><<a href="mailto:mike.hui.zhang@hotmail.com" target="_blank">mike.hui.zhang@hotmail.com</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">Hi petsc-group,<br>
<br>
if I want to map a structured mesh onto processors by partition of mesh **elements**,<br>
can I use DMDACreate3d -> DMDACreateSection, then DMDACreateLocal(Global)Vector to<br>
give me the vector corresponding to my Q1-element d.o.f.?<br></blockquote><div><br></div><div style>No, this is not finished yet. I have only implemented it in serial, mostly because no one is</div><div style>asking for it.</div>
<div style><br></div><div style>Note that the data layout has nothing to do with partitioning. We would still have to rewrite</div><div style>those routines for an element-wise partition. This is on the list of things to do.</div>
<div style><br></div><div style>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I ask because according to the manual, the usual use of DMDA is only for partition of<br>
**nodes** to processors and the local vector is not corresponding to partition of elements.<br>
For example, consider the Q1 element on the 1D mesh, (numbers for nodes)<br>
<br>
1--2--3--4--5<br>
<br>
Given two processors, I want the local vector after partition to reside on the partitioned mesh<br>
<br>
1--2--3   3--4--5<br>
proc 0    proc 1<br>
<br>
The global value of d.o.f at the node 3 may be managed by proc 1, and proc 0 only ghosts it.<br>
If it is the case, proc1 does not ghost any value.<br>
<br>
This is different from the usual DMDA for finite difference stencils, in which the ghosting<br>
is symmetric, i.e. both proc0 and proc1 will ghost some values.<br>
<br>
Thanks in advance!<br>
<br>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>