[petsc-users] global index distributed arrays

[Mohamad M. Nasr-Azadani]

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>