On 2/7/2012 2:49 PM, Matthew Knepley wrote:

On Mon, Jul 2, 2012 at 3:58 AM, TAY wee-beng <zonexo@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

I just tested 3.3-p1 and it still doesn't work (example ex11f90 in dm). Is there a chance petsc-dev can work?

 

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