On Aug 31, 2010, at 1:11 PM, Jed Brown wrote: > If a VecGhost was passed in to SNES, then all the Krylov vectors will > also be VecGhost. If not, create a new VecGhost (VecCreateGhost() or > VecDuplicate() an existing VecGhost) and VecCopy from the plain Vec to > the ghost Vec.
Well, that is odd then. Because I'm sure that libMesh is passing in ghosted vectors (hmmm... maybe the residual (rhs) isn't ghosted by default?). Maybe libMesh is failing at detecting that they are ghosted? It looks like libMesh is looking for a local to global mapping to decide if it's ghosted or not. Is that the proper way? Derek ------------------------------------------------------------------------------ This SF.net Dev2Dev email is sponsored by: Show off your parallel programming skills. Enter the Intel(R) Threading Challenge 2010. http://p.sf.net/sfu/intel-thread-sfd _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
