<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Aug 17, 2015 at 10:49 AM, Xujun Zhao <span dir="ltr"><<a href="mailto:xzhao99@gmail.com" target="_blank">xzhao99@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div>I want PETSc to generate random vector using VecSetRandom() following given examples, but failed and showed some "out of memory" error. The following is the code, which goes well until it reaches VecSetRandom(). Can anyone help me figure out the reason? Thanks a lot.</div></div></blockquote><div><br></div><div>Does src/vec/vec/examples/tests/ex43.c run for you?</div><div><br></div><div> Thanks,</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"><div dir="ltr"><div>XZ<br><br><br>--------------------------------------------------------------------------------------------<br>  Vec             u;<br>  PetscRandom     rand_ctx;     /* random number generator context */<br>  PetscMPIInt     size, rank;<br>  PetscInt        n, dn;<br>  <br><br>  MPI_Comm_rank(PETSC_COMM_WORLD,&rank);//CHKERRQ(ierr);<br>  MPI_Comm_size(PETSC_COMM_WORLD,&size);//CHKERRQ(ierr);<br>  n  = N/size + 1;<br>  dn = n*size - N;<br>  if ( dn>0 && rank<dn ) n -= 1;<br>  printf("--->test in petsc_random_vector(): rank = %d, n = %d\n",rank,n);<br><br><br>  VecCreate(PETSC_COMM_WORLD,&u);<br>  VecSetSizes(u,n,N);<br>  PetscRandomCreate(PETSC_COMM_WORLD, &rand_ctx);<br>#if defined(PETSC_HAVE_DRAND48)<br>  PetscRandomSetType(rand_ctx,PETSCRAND48);<br>#elif defined(PETSC_HAVE_RAND)<br>  PetscRandomSetType(rand_ctx,PETSCRAND);<br>#endif<br>  PetscRandomSetFromOptions(rand_ctx);<br><br><br>  VecSetRandom(u,rand_ctx);<br>  PetscRandomDestroy(&rand_ctx);<br></div></div>
</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>