> Of course, this new function sounds unnecessary because reordering v is > equivalent to reordering > idxm[] and idxn[]. But the second method seems not an easy task to me, is it?
Why not use a single AO object to hold the mapping between PETSc mapping and your maping? Seems to me that the indecies used in v[] are completely arbitrary; they should always go in [0,1,...,m*n], don't they? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120527/3c6e429a/attachment.html>