On Tue, 12 Apr 2011, John Peterson wrote:
On Tue, Apr 12, 2011 at 2:05 PM, Roy Stogner <[email protected]> wrote:I think John did most of the pyramid work, and IIRC it does work - the five corner children are pyramids (and similar to the parent since we only support 1st order pyramids) but I believe there's an inverted sixth pyramid and some tets in the middle. The really weird bit is the rational (non-polynomial!) basis functions, but that won't be an issue for fixing a topology-based method.I don't know if pyramid adaptivity actually works, but I was able to find an ancient Matlab script which draws the Pyramid children--they are basically as Roy described. I've attached the PDFs if you are interested Roy. Not sure they will make it to the list though... if anybody else wants them just let me know.
It looks like the adaptivity for Pyramids never got implemented? Pyramid5::embedding_matrix() just throws an error. Which means I'll drop it for now. If we can't refine pyramids yet then that's just a missing feature that we can add at our leisure. I'm fine with the library occasionally saying "we can't do that yet" as long as it never says "we can do that: here's a wrong answer!" Now to figure out the tets.. those we actually do refine. --- Roy
------------------------------------------------------------------------------ Forrester Wave Report - Recovery time is now measured in hours and minutes not days. Key insights are discussed in the 2010 Forrester Wave Report as part of an in-depth evaluation of disaster recovery service providers. Forrester found the best-in-class provider in terms of services and vision. Read this report now! http://p.sf.net/sfu/ibm-webcastpromo
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
