<div class="gmail_quote">On Thu, Feb 23, 2012 at 15:05, Bojan Niceno <span dir="ltr"><<a href="mailto:bojan.niceno@psi.ch">bojan.niceno@psi.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">No, I use global.<br></div></blockquote><div><br></div><div>The local form is just a local vector. It doesn't even know that a global problem exists. You can't index into it using global indices. (In general, there is no efficient way to look up information in the local vector (includes ghost points) using global indices.)</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000">
<br>
<tt> for(Int i=0; i<n; i++) {<br>
Int gi = mesh.nodes[i].global_number;<br>
VecGetValues(x, 1, &gi, &unk[i]);<br>
}<br>
</tt><br>
"n" is defined as the number of cells inside, i.e. without buffers.
"unk" is my external array. If I try to access buffer values, I
use:<br>
<br>
<tt> for(Int i=0; i<N; i++) {<br>
Int gi = mesh.nodes[i].global_number;<br>
VecGetValues(x, 1, &gi, &unk[i]);<br>
}<br>
</tt><br>
But then I end up with tons of warnings, presumably because I am
going beyond "n". Vector x was created with VecCreateGhost.</div></blockquote></div><br>