[petsc-users] global index distributed arrays

[Matthew Knepley]

On Fri, Dec 10, 2010 at 11:46 PM, Mohamad M. Nasr-Azadani
<mmnasr at gmail.com>wrote:

> 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.
>

((k*N + j)*M + i)*C + c

   Matt


> Thank,
> Mohamad
>
>
>


-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101211/7392f742/attachment.htm>