Re: Linear system partitioning

["Matthew Knepley" <knepley@xxxxxxxxx>]

We use the numbering and MatGetSubmatrix() to permute the
matrix in those examples. You can permute the vector by
creating a VecScatter between two global vectors. The first
IS is "is" from the example, and the second is 0...n-1 (use
ISCreateStride).

  Thanks,

     Matt

On 5/4/06, Pierre Gay <gay-pierre@xxxxxxxxxxxxxx> wrote:

Hello,

I'm trying to use PETSc to partition and renumber a linear system. I
found some examples in the distribution:

$(PETSC_DIR)/ksp/ksp/examples/tutorials/ex10.c
$(PETSC_DIR)/mat/examples/tests/ex73.c

But these seem to only renumber the matrix. I also need to renumber
associated vectors (right hand side, solution) and I did not manage to
get some working code.

Could someone lead me to an example?

Many thanks,

Pierre

-- 
"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness