<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On May 27, 2012, at 8:58 PM, Mohammad Mirzadeh wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><pre style=""><div>> Of course, this new function sounds unnecessary because reordering v is equivalent to reordering </div><div>> idxm[] and idxn[]. But the second method seems not an easy task to me, is it?</div>
</pre><br><div class="gmail_quote">Why not use a single AO object to hold the mapping between PETSc mapping and your maping? </div></div></blockquote><div><br></div><div>Yes, thanks! I'm doing so.</div><br><blockquote type="cite"><div dir="ltr"><div class="gmail_quote">Seems to me that the indecies used in v[] are completely arbitrary; they should always go in [0,1,...,m*n], don't they?</div>
</div>
</blockquote><br></div><div>No, I only take some non-contiguous entries of v[] to the mat. So, I need to first take out</div><div>those values to a new array, what I worried about the cost. But as Barry Smith told us, </div><div>the cost is small.</div><div><br></div><div><br></div><br></body></html>