I cannot understand your question.<br><br>  Matt<br><br><div class="gmail_quote">On Fri, Feb 27, 2009 at 6:45 PM, Yixun Liu <span dir="ltr">&lt;<a href="mailto:enjoywm@cs.wm.edu">enjoywm@cs.wm.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
<br>
  VecCreate(PETSC_COMM_WORLD,&amp;solutionDisplacementVector);<br>
<br>
  VecSetSizes(solutionDisplacementVector, 3*numOfVerticesOfOneProcessor,<br>
systemSize); //in this case 3*numOfVerticesOfOneProcessor==systemSize<br>
  VecSetFromOptions(solutionDisplacementVector);<br>
<br>
<br>
 VecDuplicate(solutionDisplacementVector, &amp;weightedDisplacementVector);<br>
<br>
<br>
for(i = 0; i &lt; systemSize; i++)<br>
{<br>
output weightedRealDisplacement[i] to file//with this code the values in<br>
weightedDisplacementVector are correct.<br>
<br>
...<br>
VecSetValue(weightedDisplacementVector, i, weightedRealDisplacement[i],<br>
ADD_VALUES);<br>
...<br>
<br>
}<br>
<br>
VecAssemblyBegin(weightedDisplacementVector);<br>
 VecAssemblyEnd(weightedDisplacementVector);<br>
<br>
output weightedDisplacementVector to file to see its values<br>
<br>
<br>
<br>
For simplicity I only use one processor. As I examine the values in<br>
weightedDisplacementVector and I found if  I add the output code in the<br>
for loop the values in weightedDisplacementVector are correct. Without<br>
it there is error.<br>
<br>
<br>
Thanks.<br>
<font color="#888888"><br>
Yixun<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>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<br>