On Wed, 2 Dec 2009, Rahul Sampath wrote: > Thanks for the response. I too use the PARALLEL type for the most > part. However, I also need access to the ghosts and I'm afraid I'm > interested in more than just a localize to a GHOSTED vector.
Rather than localize a primary PARALLEL vector to a local GHOSTED version, would it be sufficient to make your primary vector GHOSTED in the first place? Even in the library, we do the solution->current_local_solution localizations just for compatibility reasons at this point: the GHOSTED vector interface support with PETSc is pretty recent, and with Trilinos/Epetra it's nonexistent. > I guess I'm trying to understand the parallel data layout. I looked > at how the system.rhs and system.solution vectors were created. It > looks like they are created using Petsc's VecCreateMPI function. So, > can I assume that the following holds: > > 1. Each processor owns a consecutive set of global dofs. > 2. The first dof owned by any processor is equal to 1 more than the > last dof owned by the previous processor. Yes, and yes. --- Roy ------------------------------------------------------------------------------ Join us December 9, 2009 for the Red Hat Virtual Experience, a free event focused on virtualization and cloud computing. Attend in-depth sessions from your desk. Your couch. Anywhere. http://p.sf.net/sfu/redhat-sfdev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
