I think the class MeshFunction is supposed to have this kind of functionality, but it seems it quite outdated. You should write some code that suit your need. The zz error patch recovery is easy to implement using the DenseMatrix class. See the article:
Superconvergence and the superconvergent patch recovery O. C. Zienkiewicza, Corresponding Author Contact Information and J. Z. Zhub http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6V36-3YB56X8-8&_user=2459694&_rdoc=1&_fmt=&_orig=search&_sort=d&view=c&_acct=C000057393&_version=1&_urlVersion=0&_userid=2459694&md5=00bdb8afa675d4930896d95dd508330f Here is other related article: A Posteriori Error Estimation for Finite Element Methods, an Overview http://www.math.ntnu.no/preprint/numerics/2001/N6-2001.ps Paulo >Message: 5 >Date: Thu, 27 Nov 2008 17:56:04 +0200 >From: "Vasilis Vavourakis" <[EMAIL PROTECTED]> >Subject: [Libmesh-users] regarding post-process... >To: [email protected] > >Hello all :) > >I attached a source file that implements a displacement-based formulation >for steady-state 2D elasticity (very simple example for someone to >understand it) and this piece of code works FINE. For my practice i also >implemented a mixed formulation (displacements-stresses) and i did manage to >evaluate them correctly. If there is anyone that may need it i can post it >as well!!! > >My problem in the displacement formulation (see source code file) is that >after the numerical evaluation of displacements i couldn't find a way to >calculate nodal stresses in the post-process. > >Although the same question it was posted both from me and others in the past >( see http://sourceforge.net/mailarchive/forum.php?thread_name=.139.91.194.170.12269 >27623.squirrel%40www.mead.upatras.gr&forum_name=libmesh-users), >i still cannot implement the ExplicitSystem class so as to do my >post-process! I red over & over again the documentation for this class but >it was really hard for me to figure out a way to do it! :( Unfortunately i >got very confused due to the fact that i'm bery new to LibMesh as well as >its "vast" size of documentation for classes... > >I'd really appreciate if someone can help me out with this... >Apart from useful advices, any handy source code is really really welcomed >:) > >many thanks in advance ~:^) ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
