Hi guys, 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. To be more more specific, I am trying to setup a matrix based on a three-dimensional DA. (Star stencil, width=1). For some special nodes, I need to insert nonzero values which do not fit in the local plus the ghost regions of the DA. 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. I tried to use DAGetGlobalIndices(), but that would again return the global indices of the local plus ghost nodes on current processor. 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.
Thank, Mohamad -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101210/4e73476b/attachment.htm>
