On Tue, 31 Aug 2010 15:13:02 -0500 (CDT), Roy Stogner <royst...@ices.utexas.edu> wrote: > We don't want to create ghosted instead of serial vectors - sometimes > code is going to ask for a serial vector because it really needs data > that's not on the local processor or its ghost cells. (and because > it's too lazy to build a proper extended send list...)
That wouldn't make sense because ghosted vectors are parallel vectors, not serial vectors. > But it would be sufficient for your case to just create ghosted > instead of parallel vectors. We can't do that in general in > NumericVector::build() because that doesn't have access to DofMap's > send_list. But we could do that as an option (run-time, even - make > it a flag in the System class) in System::add_vector() and > *System::init(). I don't remember how this is done in libmesh, but it would be typical to create a solution/residual vector from the function space (mesh), and the overlap (distribution and ghost indices) is available in that context. Jed ------------------------------------------------------------------------------ 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
