I just recall what I asked before. Is the approach vertex_index*dofs_per_vertex+local_dof
sustainable? I believe there are plans to abandon 'common cell' concept, are there? Then `dofs_per_vertex` could (I'm not sure) become variable across the mesh. Check also the problem of ghost dofs. These are covered by `dof_to_vertex_map` by negative and greater-than-local-range values. As a result this is not invertible and ghosts can't be retained in an inverse `vertex_to_dof_map`. Wouldn't it be better to use some more robust map types? Jan _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
