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@fenics.org > http://www.fenics.org/mailman/listinfo/dolfin-dev _______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev