Hello,
I use the following routine to transfer the petsc result back to my application
which works fine using one core as the local parts of the linear system ==
global parts:
VecGetArray(petsc_x, &array);
for (i=rstart; i < rend; ++i) {
application_x[i] = array[i];
}
VecRestoreArray(petsc_x, &array);
As this updates only the local portion of x on each core but nowhere the
entire/global application_x vector, multicore simulations don't work correctly.
How can I update the entire/global application_x vector with all local parts on
all cores?
Klaus
