Re: [petsc-users] Accessing Vector's ghost values

On Fri, Feb 24, 2012 at 04:43, Mohammad Mirzadeh <[email protected]> wrote:

Thanks Bojan. I suppose that would work. Problem arises, however, when the actual application numbering matters and you have to preserve it somehow. One example, is where you have a library that is currently sequential and a lot of functionality is directly based on the data structure used for the grid. For example when computing derivatives or performing interpolation and so on.

That serial library probably needs the local numbering anyway, if you give it the global indices, it won't be able to access those locations (at least not without a translation step that they did not write). It's better to build the data structure in terms of local indices.