On Tue, 31 Aug 2010, Derek Gaston wrote:
> On Aug 31, 2010, at 2:24 PM, Roy Stogner wrote: > >> However when we're just building a parallel vector, we only call >> NumericVector::build() with the global and local sizes, not with the >> entire context that exists in the caller. The best solution is to add >> a special case alternative that uses the ghosted version of build() >> for those half dozen calls, not to try and figure out how to break our >> own encapsulation within the non-ghosted build() call. > > I'm open for whatever solution. If you tell me what you want I'll take a > crack at an implementation. The effect you want should be to add a get_send_list() argument and change PARALLEL to GHOSTED in system.C lines 230, 246, 262, 310, 593. Try that first and see if it works. If it turns out to be what you need, then we could add it as an option to the library; maybe create a new System::init_parallel_vector() method that gets called on those lines and either inits a PARALLEL by default or a GHOSTED if some System::all_ghosted_vectors boolean is true? --- Roy ------------------------------------------------------------------------------ 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
