Benjamin:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
I am trying to solve a Poisson equation several times with random right-hand side vectors in order to do parallel scalability tests.<br>
Here is part of the code that I use to generate a random vector:<br></blockquote><div> </div><div>This is intended. You can use option ' -random_seed <seed>' to change it.</div><div>See <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PetscRandomGetSeed.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PetscRandomGetSeed.html</a></div>
<div><br></div><div>Hong</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
PetscRandom    :: rctx<br>
<br>
...<br>
<br>
do I = 1,n<br>
<br>
<br>
call PetscRandomCreate(PETSC_COMM_WORLD,rctx,ierr); CHKERRQ(ierr)<br>
call VecSetRandom(f,rctx,ierr); CHKERRQ(ierr)<br>
call PetscRandomDestroy(rctx,ierr); CHKERRQ(ierr)<br>
<br>
call VecView(f,PETSC_VIEWER_STDOUT_WORLD,ierr); CHKERRQ(ierr)<br>
<br>
call Poisson<br>
<br>
end do<br>
<br>
<br>
It appears that f does not change during the execution of the do-loop. In fact its value is even always the same for I=1 when I run the code several times. Apparently I am missing something. Can anybody help?<br>
<br>
Regards,<br>
<br>
<br>
Benjamin<br>
</blockquote></div><br>