[petsc-users] Reshaping a vector into a matrix

[S V N Vishwanathan]

> > 
> > In my application I need access to (x_i.vec1, x_i.vec2, ...., x_i.veck)
> > on the same processor. If I use the local size provided by the
> > MatGetLocalSize how do I ensure that the local vectors are aligned w.r.t
> 
>    After you do the MatLoad() of the AIJ matrix get its local size, say nlocal. The local size of the MAIJ matrix created from the AIJ matrix will be k*nlocal. 
> You will then create your vector with a local size of k*nlocal. I think that everything lines up perfectly. 
> If this is not the case then please send a demonstration of its failure to petsc-maint at mcs.anl.gov (code) and we'll see what is going on.

This is indeed the case but I was not sure if it was just fluke or by
design. Thanks for the clarification.

vishy