[petsc-users] Problem with DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

[TAY wee-beng]

Hi,

My code hangs and I added in mpi_barrier and print to catch the bug. I 
found that it hangs after printing "7". Is it because I'm doing 
something wrong? I need to access the u,v,w array so I use 
DMDAVecGetArrayF90. After access, I use DMDAVecRestoreArrayF90.

        call DMDAVecGetArrayF90(da_u,u_local,u_array,ierr)
        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"3"
        call DMDAVecGetArrayF90(da_v,v_local,v_array,ierr)
        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"4"
        call DMDAVecGetArrayF90(da_w,w_local,w_array,ierr)
        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"5"
        call 
I_IIB_uv_initial_1st_dm(I_cell_no_u1,I_cell_no_v1,I_cell_no_w1,I_cell_u1,I_cell_v1,I_cell_w1,u_array,v_array,w_array) 

        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"6"
        call DMDAVecRestoreArrayF90(da_w,w_local,w_array,ierr) !must be 
in reverse order
/*        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"7"*/
        call DMDAVecRestoreArrayF90(da_v,v_local,v_array,ierr)
        call MPI_Barrier(MPI_COMM_WORLD,ierr);  if (myid==0) print *,"8"
        call DMDAVecRestoreArrayF90(da_u,u_local,u_array,ierr)

--
Thank you.

Yours sincerely,

TAY wee-beng