Hi Unless you are using elements with higher-order continuity the divergence of your solution will not be continuous at support points. You can export the average value of the divergence per cell perhaps?
Andrew On 06.04.2011, at 11:43, Timo Koellner wrote: > Hi everyone, > > after you helped me to get rid of my last problem, here's a new one... No big > thing. > > After having solved the system I'm interested in the divergence of my > solution, ideally at the same points as the actual solution is computed and > saved to e.g. a gnuplot file by DataOut (I assumed them to be the DoFs or are > they just the vertices?) > > I found the get_function_gradients member of FEValuesBase and this actually > works. But using this I get the divergence on each quadrature point, but not > on the same points as DataOut prints the solution. > > Is there a way to obtain this? > > Thanks. > > Timo > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
