On Wed, 1 Sep 2010, Derek Gaston wrote:
> add a member variable to system called
> "default_parallel_vector_type" that defaults to PARALLEL.
>
> add_vector will get an optional argument for the ParallelType
> (defaulting to PARALLEL).
>
> add_vector will also store off the ParallelType for each vector in a
> map (I _do_ think this is necessary because there is no other way to
> capture the ParallelType within add_vector).
>
> All internal vectors will be added using add_vector and passing in
> "default_parallel_vector_type" (or GHOSTED or PARALLEL if only those
> make sense... like for current_local_solution).
>
> Finally, I'm also going to create a helper function in System called
> init_vector that just takes the name of the vector to be inited...
> this is because there are a few places where vectors get explicitly
> inited.
>
> I think that should cover it. If you _want_ to you will be able to
> override the defaults by either setting default_parallel_vector_type
> OR by doing as you mentioned and calling add_vector yourself with a
> different type.
This all sounds good except for the map - what would
paralleltype_map["vector_name"] get you that
get_vector("vector_name").type() doesn't?
---
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
