On Wed, 7 Oct 2009, Chris Petrovitch wrote: > I'm about to start coding a discontinuous galerkin code to simulate > seismic wave propagation, and I was wondering if libmesh can handle > discontinuous elements?
Yes. See the XYZ and MONOMIAL FE types. > if it is possible, are there any example codes I could see? or code > snippets on how to access the element faces, face normals.. etc. We don't currently have any example DG codes; IIRC the only primary developer to use it was Ben Kirk, for a compressible flow problem too complex to boil down into a simple example. Look at the boundary conditions on the existing examples, though; the get_normals(), etc. methods work on inter-element boundaries as well as on domain boundaries. If you're planning to do adaptive refinement things get a little trickier. The JumpErrorEstimator class is an example of how to properly integrate inter-element boundary jump terms across non-conforming elements. --- Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
