[petsc-users] confusion about VecGetArray on a sequential vec

丁老师 ztdepyahoo at 163.com
Thu Aug 1 06:08:25 CDT 2013


but when we use vecveiw to veiw the seqx, it gives value only with 1 process






At 2013-08-01 18:36:09,"Dave May" <dave.mayhem23 at gmail.com> wrote:

When you called
  VecCreateSeq(PETSC_COMM_SELF,5,&seqx);

you created a sequential vector on each process - not just on rank 0.



On 1 August 2013 12:11, 丁老师 <ztdepyahoo at 163.com> wrote:

I create a sequential vec seqx, and i get it value with VecGetArrary, since seqx is created  only in process 0, why every process can output the value of the seqx.



 VecCreateSeq(PETSC_COMM_SELF,5,&seqx);

 VecView(seqx,PETSC_VIEWER_STDOUT_WORLD);
 double *seqarry;
 VecGetArray(seqx,&seqarry);


 for(int i=0;i<n;i++)
   cout<<MyRank<<"  seqarry "<<s eqarry[i]<<endl;










-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130801/a3bc73a1/attachment.html>


More information about the petsc-users mailing list