Andrew, first: many thanks for tackling this! I've taken your patches and checked them in, with the following small adjustment: I renamed your FEValuesViews::SymmetricSecondOrderTensor<dim> into FEValuesViews::SymmetricTensor<2,dim> since we already have a SymmetricTensor<2,dim> class in namespace dealii to which your new class corresponds.
> The problem is when I run the test file the code simply hangs when I > initialise the FEValues > FEValues<dim> fe_values(fe, quadrature_formula, update_values | > update_gradients | update_quadrature_points | update_JxW_values); That is interesting because it doesn't appear to happen for me. If you update the library (which should give you the files you sent to me, plus the change mentioned above), do make clear ; make all does it still hang? I sometimes find that programs hang or segfault when not all parts of the library have been compiled at the same time. If it still hangs, run the program in gdb, let it run for a while and then interrupt it -- what's the backtrace then in gdb? > I have decided to only add the functionality for get_values and > get_divergences as I don't need anything else and it becomes far more > complicated when one starts asking for gradients of second order > tensors etc. That's fine for the moment. We can add the rest when someone finds a need for it. Best Wolfgang ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
