<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Jul 21, 2015 at 1:35 AM, TAY wee-beng <span dir="ltr"><<a href="mailto:zonexo@gmail.com" target="_blank">zonexo@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I need to check the contents of the array which was declared using:<br>
<br>
PetscScalar,pointer :: u_array(:,:,:),v_array(:,:,:),w_array(:,:,:),p_array(:,:,:)<br>
<br>
I tried to use :<br>
<br>
call PetscViewerASCIIOpen(MPI_COMM_WORLD,"pres.txt",viewer,ierr)<br>
<br>
call VecView(p_array,viewer,ierr)<br>
<br>
or<br>
<br>
call MatView(p_array,viewer,ierr)<br>
<br>
call PetscViewerDestroy(viewer,ierr)<br>
<br>
but I got segmentation error. So is there a PETSc routine I can use?</blockquote><div><br></div><div>No. Those routines work only for Vec objects. You could</div><div><br></div><div> a) Declare a DMDA of the same size</div><div><br></div><div> b) Use DMDAVecGetArrayF90() to get out the multidimensional array</div><div><br></div><div> c) Use that in your code</div><div><br></div><div> d) Use VecView() on the original vector</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
-- <br>
Thank you<br>
<br>
Yours sincerely,<br>
<br>
TAY wee-beng<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>