Gabriel, [...] > u_x.reinit(n_dofs_in_one_direction); > *FEValues<dim> fe_values(fe, > quadrature_formula, update_values);* > *std::vector<Tensor<1, dim>> velocity_values(n_q_points);* > *const FEValuesExtractors::Vector velocities(0);* > > *for (const auto &cell : dof_handler.active_cell_iterators())* > * {* > * fe_values.reinit(cell);* > * fe_values[velocities].get_function_values(solution, > velocity_values);* > > * [...]* > > *}* > > With this procedure I get the vector of Tensors velocity_values, which > contains the velocity-values of "solution" in every direction at the > quadratue points. So now I an easily get the values of velocity_values in > the x-direction at the quadrature points. But how do I get the values of > u_x at the dofs? > > Does somebody have an idea how to fix this, or a better way to do handle > this? >
You can use fe.system_to_component_index(i).first; to obtain the component of the local degree with number i, see https://www.dealii.org/current/doxygen/deal.II/step_8.html#ElasticProblemassemble_system . 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAOYDWbJBw8LSCkjPtcQCJWgA89bdkpekUhD1CA1%2Bm4HxQZHr4Q%40mail.gmail.com.