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]? Johan
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
