Ted Sariyski <[email protected]> writes: > On 10/31/2013 04:53 PM, Jed Brown wrote: >> Ted Sariyski <[email protected]> writes: >>> Q1: I have the stiffness matrix and the rhs vector assembled (outside) >>> and in PETSc order. Do I still need to create a scatter:"ierr = >>> VecScatterCreate(bb,isglobal,localX,islocal,&scatter)"? >> What do you want to do with the Vec? > Use KSP for solving a linear system.
It looks like your vectors are already distributed, so you don't need a VecScatter. If you were computing nonlinear residuals, you would probably need to communicate with neighbors; VecScatter is a flexible way to do that.
pgp1SSEmuHHm3.pgp
Description: PGP signature
