We use the numbering and MatGetSubmatrix() to permute the<br>
matrix in those examples. You can permute the vector by<br>
creating a VecScatter between two global vectors. The first<br>
IS is "is" from the example, and the second is 0...n-1 (use<br>
ISCreateStride).<br>
<br>
Thanks,<br>
<br>
Matt<br><br><div><span class="gmail_quote">On 5/4/06, <b class="gmail_sendername">Pierre Gay</b> <<a href="mailto:gay-pierre@numericable.fr">gay-pierre@numericable.fr</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>Hello,<br><br>I'm trying to use PETSc to partition and renumber a linear system. I<br>found some examples in the distribution:<br><br>$(PETSC_DIR)/ksp/ksp/examples/tutorials/ex10.c<br>$(PETSC_DIR)/mat/examples/tests/ex73.c
<br><br>But these seem to only renumber the matrix. I also need to renumber<br>associated vectors (right hand side, solution) and I did not manage to<br>get some working code.<br><br>Could someone lead me to an example?<br>
<br>Many thanks,<br><br>Pierre<br><br><br><br></blockquote></div><br><br clear="all"><br>-- <br>"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness