Thanks Daniel, I had confused the dof_index with the vertex_number. Now It works as expected.
Best, Joel On Wednesday, September 7, 2016 at 3:47:21 PM UTC+2, Daniel Arndt wrote: > > Joel, > > In GridTools::find_cells_adjacent_to_vertex(dof_handler,center_id) you > need to give the number of the vertex_number as center_id, not the number > of the degree of freedom. > Then, you want to use > > Quadrature<dim> > quadrature_formula(dof_handler.get_fe().get_unit_support_points()); > > to create a Quadrature object that gives you the support points on each > cell. > Finally, you can output the global dof number and their support point on > each cell by: > > for (unsigned int j=0; j<dofs_per_cell; ++j) > std::cout << local_dof_indices[j] << ": "<< > fe_values.quadrature_point(j) << std::endl; > > In this ouput, you can then search for the dofs with support_point nearest > to the vertex you wanted to consider. > > Best, > Daniel > -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
