> > > > 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
