On Thu, 27 Nov 2014 20:38:56 +0100 Johan Hake <[email protected]> wrote:
> Hello! > > In some code I have I uses the indices interface to set local dofs in > a vector. It turns out that v[indices] = some_values uses the > GenericVector::set function instead of GenericVector::set_local. This > means that one need to pass global indices. > > I typically use the slicing together with some combination of indices > I got from the vertex_to_dofs functionality. However, now that > returns local dofs and it then makes more sense to switch the > behavior of v[indices] to use local dofs. > > Any objections against switching to local indices in v[indices]? Maybe take into account an issue with GenericVector.__setitem__ https://bitbucket.org/fenics-project/dolfin/issue/54 when fixing this. Jan > > Johan _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
