On Mon, Jul 2, 2012 at 3:58 AM, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I have used DMDACreate2d for my code and then use:<br>
<br>
call DMLocalToGlobalBegin(da,b_rhs_<u></u>semi_local,INSERT_VALUES,b_<u></u>rhs_semi_global,ierr)<br>
<br>
call DMLocalToGlobalEnd(da,b_rhs_<u></u>semi_local,INSERT_VALUES,b_<u></u>rhs_semi_global,ierr)<br>
<br>
to construct the global DM vector b_rhs_semi_global<br>
<br>
Now I want to get the values with ghost values in a 2d array locally which is declared as:<br>
<br>
real(8), allocatable :: array2d(:,:)<br>
<br>
I guess I should use DMDAGetGhostCorners to get the corressponding indices and allocate it.  But what should I do next? How can I use something like VecGetArrayF90 to get to the pointer to access the local vector?<br>
<br>
I can't use DMDAVecGetArrayF90/<u></u>DMDAVecRestoreArrayF90 since I'm using intel fortran and they can't work. I can't use gfortran at the moment since I've problems with HYPRE with gfortran in 3D.<br>
</blockquote><div><br></div><div>Are you certain of this? That used to be true, but the current version should work for any F90.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Yours sincerely,<br>
<br>
TAY wee-beng<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>