<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Dec 20, 2014 at 7:51 AM, Victor Magri @Sinmec <span dir="ltr"><<a href="mailto:victorapm@sinmec.ufsc.br" target="_blank">victorapm@sinmec.ufsc.br</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello!<div><br></div><div>I'm using MatPermute function to create a permuted matrix from two given index sets. According to line 4711 of ${PETSC_DIR}/src/mat/interface/matrix.c they map from row/col of permuted matrix to row/col of the original one. So, the index sets are inverted before the act of permutation. Shouldn't it be better to call this function MatInvertPermute or something like this in oppose to MatPermute? Besides that, I think it would be good to define a new function which does not invert the index sets before permutation. </div></blockquote><div><br></div><div>The old ordering routines of Alan George, which we wrap up in MatOrdering give back</div><div><br></div><div>  p[new number] = old number</div><div><br></div><div>Thus this is what we use for input. We get the inverse, because for parallelism, it makes</div><div>sense to have</div><div><br></div><div>  p[old number] = new number</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Thank you!</div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>