[petsc-users] Getting 2d array with updated ghost values from DM global vector

[Matthew Knepley]

On Mon, Jul 2, 2012 at 3:58 AM, TAY wee-beng <zonexo at gmail.com> wrote:

> Hi,
>
> I have used DMDACreate2d for my code and then use:
>
> call DMLocalToGlobalBegin(da,b_rhs_**semi_local,INSERT_VALUES,b_**
> rhs_semi_global,ierr)
>
> call DMLocalToGlobalEnd(da,b_rhs_**semi_local,INSERT_VALUES,b_**
> rhs_semi_global,ierr)
>
> to construct the global DM vector b_rhs_semi_global
>
> Now I want to get the values with ghost values in a 2d array locally which
> is declared as:
>
> real(8), allocatable :: array2d(:,:)
>
> 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?
>
> I can't use DMDAVecGetArrayF90/**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.
>

Are you certain of this? That used to be true, but the current version
should work for any F90.

   Matt


> Thanks
>
> --
> Yours sincerely,
>
> TAY wee-beng
>
>


-- 
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/20120702/91769539/attachment.html>