
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Junchao,<div class=""><br class=""></div><div class="">I tested by commenting out the AOApplicationToPetsc calls as you suggest, but it doesn’t work because it doesn’t maintain the proper order of the elements in the scattered vectors.</div><div class=""><br class=""></div><div class="">I attach a modified version of the test code where I put elements into the global vector, then carry out the scatter, and check on the subcomms that they are correct.</div><div class=""><br class=""></div><div class="">You can see everything is fine with the AOApplicationToPetsc calls, but the comparison fails when those are commented out.</div><div class=""><br class=""></div><div class="">If there is some way I can achieve the right VecScatters without those calls, I would be happy to know how to do that.</div><div class=""><br class=""></div><div class="">Thank you again for your help.</div><div class=""><br class=""></div><div class="">Randy</div><div class=""><br class=""></div><div class="">ps. I suggest you run this test with nx=ny=nz=10 and only a couple subcomms and maybe 4 processes to demonstrate the behavior</div></body></html>