Hi all, I have a doubt regarding the mesh movement in step42.cc Below is the code taken from the example.
Lets say cell 1 and cell 2 are two adjacent cells of a fem mesh (share 2 nodes / vertex) and they live on processor 1 and 2 respectively. Processor 1 only keeps track of vertices that it has moved via the vertex touched vector. It has no record of whether that vertex is moved by any other processor or not. In this case, when the processor 1 moves all the vertices of cell 1 (which includes the shared vertices), I am unable to see why the processor 2 wont move that shared vertex. template <int dim> void PlasticityContactProblem<dim>:: move_mesh (const TrilinosWrappers::MPI::Vector <http://dealii.org/8.4.0/doxygen/deal.II/classTrilinosWrappers_1_1MPI_1_1Vector.html> &displacement) const { std::vector<bool> vertex_touched(triangulation.n_vertices(), false); for (typename DoFHandler<dim>::active_cell_iterator <http://dealii.org/8.4.0/doxygen/deal.II/classDoFHandler.html> cell = dof_handler.begin_active <http://dealii.org/8.4.0/doxygen/deal.II/classDoFHandler.html#ad2df51a32781906218910731e063ac62> (); cell != dof_handler.end <http://dealii.org/8.4.0/doxygen/deal.II/classDoFHandler.html#a05c70e1862a2ad145f91f9da1f44cc28> (); ++cell) if (cell->is_locally_owned()) for (unsigned int v = 0; v < GeometryInfo<dim>::vertices_per_cell; ++v) if (vertex_touched[cell->vertex_index(v)] == false) { vertex_touched[cell->vertex_index(v)] = true; Point<dim> <http://dealii.org/8.4.0/doxygen/deal.II/classPoint.html> vertex_displacement; for (unsigned int d = 0; d < dim; ++d) vertex_displacement[d] = displacement(cell->vertex_dof_index(v, d)); cell->vertex(v) += vertex_displacement; } } -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
