Update: I switched to my old trick of nodal averaging method: take the average of the values at quadrature points, then add this average to the vertex dofs. And if a vertex is shared by n cells, divide the values at the related dofs by n. I am happy with the results because they look quite smooth.
The only concern is that: for a stress component, say Sxx, appears in my resulting vtu files as Sxx_0, Sxx_1 (in 2D) because the variable Sxx is a vector of size n_dofs as required by the add_data_vector function. Is there a way to get rid of the unnecessary duplicates? Essentially I want Sxx to be a vector of size n_vertices. Thank you Jie -- 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.