<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:΢ÈíÑźÚ
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Dear Petsc developers,<div><br></div><div>I have a vector M which consists of a series of 3d vectors, and I want to reset M by normalizing each 3d vector. Here is my code:<br><div><br></div><div>/**********************************************************************</div><div><div>VecGetArray(M, &aM); </div><div>DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0); </div><div><span class="Apple-tab-span" style="white-space:pre"> </span></div><div>for(node=xs; node<xs+xm; node++)</div><div><span style="font-size: 12pt;">{</span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>mag=0.0;</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>for(l=0; l<3; l++)</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>{</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>val[l]=aM[node*3+l];</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>mag=mag+val[l]*val[l];<span class="Apple-tab-span" style="font-size: 12pt; white-space: pre;"> </span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>}</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>mag=sqrt(mag);</div><div><span class="Apple-tab-span" style="white-space:pre"> </span><span style="font-size: 12pt;">for(l=0; l<user->NDOF; l++)</span></div><div><span class="Apple-tab-span" style="white-space:pre"> </span>{</div><div><span class="Apple-tab-span" style="white-space:pre"> </span> aM[node*3+l]=val[l]/mag;</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>}</div><div>}</div><div>VecRestoreArray(M, &aM); </div><div>VecAssemblyBegin(M);</div><div>VecAssemblyEnd(M); </div><div>VecView(M, PETSC_VIEWER_STDOUT_WORLD);</div></div></div><div>**********************************************************************/</div><div><br></div><div>but I got the error at the last step:</div><div><div>--------------------------------------------------------------------------</div><div>mpiexec noticed that process rank 3 with PID 17156 on node compute-21-8.local exited on signal 6 (Aborted).</div><div>--------------------------------------------------------------------------</div></div><div>and if I commented out VecView, and used vector M for other operations, e.g.</div><div>KSPSolve(ksp, M, b), I got "memory corruption" message. Your comment on this issue is well appreciated.</div><div><br></div><div><br></div><div><br></div> </div></body>
</html>