Hi, right now I'm implementing the function MatSolve for MATMPIAIJ matrices, so I can assume that output vector x is an MPI vector. the problem is that my library returns the solution in an array to all my MPI process, so how do I know which elements of the solution should be safe in the local vector x. Thanks