Linear system partitioning

[Matthew Knepley]

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 at numericable.fr> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20060504/440c815f/attachment.htm>