I implemented an exact solution for the mean curvature for ellipsoids
(it's a slightly ugly expression, I just took it from
here: http://math.stackexchange.com/a/540820), and I computed the L2 and
Linfty errors between my computed values and that expression. The
result: Convergence wrt mesh refinement - about an order of magnitude
decrease in each norm per increment of refine_global(). That's not a
tidy way to put it, but it was quick and easy.
That's a fairly safe check. Well done!
I think this would actually make for a really nice addition to the set
of tutorials, or the code gallery at
http://dealii.org/code-gallery.html
<http://dealii.org/code-gallery.html>
Would you be interested in getting it there? (Maybe after rewriting it
in such a way that you solve for all three components at once?)
Yes, I am very interested in doing this ( but it will have to wait until
September), thank you for asking.
That would be fantastic. Send me an off-list email when you are willing
to go forward with this!
Cheers
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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.
