On Wed, Nov 12, 2008 at 4:38 PM, Roy Stogner <[EMAIL PROTECTED]> wrote: > > On Wed, 12 Nov 2008, John Peterson wrote: > >> 1.) Just remember that higher-than-normal order quadrature on pyramids >> is required and your answer may be inaccurate. AKA "do nothing" :-) > > Not a good answer, but probably a fair one. > >> 2.) Redefine, within the pyramid quadrature rules, the meaning of >> order. I.e. return a rule several orders higher than what the user >> requests. > > Maybe an "extra_pyramid_quadrature_order()" option? So we can give it > a safe default but still let users override it for efficiency if they > know what they're doing.
Yeah, this sounds like a good idea, at least until something useful comes out of option 3... Maybe a flag in QBase, similar to the way we let knowledgable users disallow the use of rules with negative weights? >> 3.) Research quadrature rules for rational functions. I have a few >> papers on this but haven't looked into it too much yet. Phillipe >> Devloo may be doing something special in his library, so I will check >> there as well... > > This could be ideal. -- John ------------------------------------------------------------------------- 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-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
