<div class="gmail_quote">On Sat, Dec 11, 2010 at 01:40, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div class="im"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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></div><div>((k*N + j)*M + i)*C + c</div></blockquote></div><div><br></div><div>Matt, you have seriously missed the point.</div><div><br></div><div>Mohamad, there is not an easy way to calculate it for an arbitrary index, it necessarily involves a search. Do you need it for arbitrary indices, or do you know in advance which global indices are needed? DAGetAO() can give you what you are after, but there is usually a way to avoid using that code since it is not scalable.</div>
<div><br></div><div>Jed</div>