> > Using the postprocessor approach according to step-33 I have computed the > strain tensor successfully by means of the > Postprocessor<dim>::compute_derived_quantities_vector function and its > displacement gradient duh. > > Unfortunately, within an arbitrary function called right after I solve the > system, exactly before postprocessing, I try the same, but obtain different > results: > > const QGauss<dim> quadrature_formula(2); > > const unsigned int n_q_points = quadrature_formula.size(); > > FEValues<dim> fe_values(fe, quadrature_formula, update_values | > update_gradients | update_quadrature_points | update_JxW_values); > [...] > Using this FEValues object, you indeed compute the gradients at the Gauss quadrature points. However, DataPostprocessor provides you with the gradients at the patch vertices. In particular, you should be able to see that fe_values.get_quadrature_points() is different from evaluation_points in DataPostprocessor::compute_derived_quantities_*.
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.
