<html><head></head><body><span>Hello,<br><br>I use the following routine to transfer the petsc result back to my application which works fine using one core as the local parts of the linear system == global parts: <br><br> VecGetArray(petsc_x, &array);<br> for (i=rstart; i < rend; ++i) { <br> application_x[i] = array[i];<br> }<br> VecRestoreArray(petsc_x, &array);<br><br>As this updates only the local portion of x on each core but nowhere the entire/global application_x vector, multicore simulations don't work correctly.<br><br>How can I update the entire/global application_x <span><span>vector</span></span> with all local parts on all cores?<br><br></span><div><span>Klaus</span></div><div class="ydpc2dc3413yahoo-style-wrap" style="font-family: Helvetica Neue, Helvetica, Arial, sans-serif; font-size: 10px;"></div></body></html>