On Sat, Oct 03, 2009 at 12:12:34PM +0100, Garth N. Wells wrote:
Anders Logg wrote:
Why are there two different assertions on the vector size on the
Function constructors:
assert(V->dofmap().global_dimension() == x.size());
assert(V->dofmap().global_dimension() <= x->size());
Is this a typo? I hope so. Otherwise it's strange that the behavior
should be different depending on whether x is a pointer or reference.
They should actually have been '<='. This is to handle sub-function
'views' where the sub-function points to the 'mother' vector but has
it's own function space (the sub-space).
There's still one == there. Should that one also be changed for
consistency?
