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>
