On Fri, Dec 10, 2010 at 11:46 PM, Mohamad M. Nasr-Azadani <span dir="ltr"><<a href="mailto:mmnasr@gmail.com">mmnasr@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi guys, <div><br></div><div>I was wondering if there is an easy way of accessing the global index of a node which is not within the local and ghost node regions on every processor for DA. </div><div>To be more more specific, I am trying to setup a matrix based on a three-dimensional DA. (Star stencil, width=1). </div>
<div>For some special nodes, I need to insert nonzero values which do not fit in the local plus the ghost regions of the DA. </div><div>I know that I can not use MatSetValuesStencil anymore, but I still can use MatSetValues. But I need to know the global index of those nodes. </div>
<div>I tried to use DAGetGlobalIndices(), but that would again return the global indices of the local plus ghost nodes on current processor. </div><div>I know that I could use some MPI commands to pass those indices among processors, but I was wondering if there a clean and neat way with which each processor can have access to the global index for any given 3-d index. </div>
</blockquote><div><br></div><div>((k*N + j)*M + i)*C + c</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div>Thank, </div>
<div>Mohamad</div><font color="#888888"><div><br></div><div><br></div>
</font></blockquote></div><br><br clear="all"><br>-- <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<br>