Dear Jean-Paul Pelteret, Thank you very much for your help. It is very helpful.
> However, you'll note that this function only makes sense (and its > therefore only defined) for a vector field. Your solution field is scalar > though, which presents a bit of an issue. On the face of it you'd have to > first project the gradient of the solution (scaled by the function "f") > onto an axillary vector-valued FE space and then compute the divergence of > that projected field. > I'm sorry to tell you that actually I don't know how to project grad_solution into the vector-valued FE space.... Do you happen to know some example tutorial code that includes it....? Thank you. Kyusik. -- 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.