dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

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).

Garth

> --
> Anders
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev