On Wed, 1 Sep 2010, Derek Gaston wrote:

> On Sep 1, 2010, at 4:36 PM, Roy Stogner wrote:
>
>> Initialize the vector e.g. as GHOSTED with size 0 and an empty ghost
>> index set?
>
> Really?  You would really prefer that?  Seems like a hack to me.  We
> store other info about vectors this same way... like whether or not
> to project them.  Should we put that information into NumericVector
> as well?

A map is okay, if you prefer.  I just don't like having multiple
places holding different copies the same data.  But I suppose "what's
the current state of the vector" and "what should the state of the
vector become" aren't exactly the same data.

> BTW - is it ok if I add_vector() solution and current_local_solution
> other than in the constructor?

Yes.

> Are there pieces of code that actually depend on solution and
> current_local_solution being built but not initialized in the
> constructor like they are now?

Possibly, but I certainly hope not.
---
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

Reply via email to