<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hi,<br>
<br>
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.<br>
<br>
call DMDAVecGetArrayF90(da_u,u_local,u_array,ierr)<br>
call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0) print
*,"3" <br>
call DMDAVecGetArrayF90(da_v,v_local,v_array,ierr)<br>
call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0) print
*,"4"<br>
call DMDAVecGetArrayF90(da_w,w_local,w_array,ierr)<br>
call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0) print
*,"5"<br>
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)
<br>
call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0) print
*,"6"<br>
call DMDAVecRestoreArrayF90(da_w,w_local,w_array,ierr)
!must be in reverse order<br>
<i><b> call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0)
print *,"7"</b></i><br>
call DMDAVecRestoreArrayF90(da_v,v_local,v_array,ierr) <br>
call MPI_Barrier(MPI_COMM_WORLD,ierr); if (myid==0) print
*,"8"<br>
call DMDAVecRestoreArrayF90(da_u,u_local,u_array,ierr)<br>
<pre class="moz-signature" cols="72">--
Thank you.
Yours sincerely,
TAY wee-beng</pre>
</body>
</html>