Re: [Libmesh-devel] DenseSubVector

[Roy S. <roy...@ic...>]

On Mon, 9 Mar 2009, Tim Kroeger wrote:

> What do you guys think about having assignment operators in DenseSubVector 
> and DenseSubMatrix?

I'll wait to hear what Ben and John think, but operator= in the
STL container sense would be too misleading, in my opinion.  If I set
up a DenseSubVector with parent a and an identically sized
DenseSubVector with parent b, I would expect "asub = bsub" to change a
slice of a to have values equal to a slice of b.  Your patch looks
like it just resets asub's parent to be b and asub's indexing to match
bsub?  That could lead to silent numerical errors way too easily.

> Actually, I would like to have them since that would 
> enable me to use a std::vector<DenseSubVector>.

Could you use a std::vector<DenseSubVector*> instead?  Or maybe
std::vector<boost::shared_ptr<DenseSubVector> >, if you don't want to
worry about memory management.
---
Roy