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>
<div><br></div><div>Thank, </div><div>Mohamad</div><div><br></div><div><br></div>