On Wed, Jun 18, 2008 at 5:18 PM, John Peterson <[EMAIL PROTECTED]> wrote: > On Wed, Jun 18, 2008 at 1:52 PM, David Knezevic <[EMAIL PROTECTED]> wrote: >> >> hehe, well the way the paper is written doesn't inspire me with confidence >> in their results, so I'd be very interested to hear the results of the >> comparisons. > > So far, so good. This 14-point rule described by Walkington returns > the same results for polynomials of fifth order as the 15 point rule > we already have in LibMesh. Can anyone suggest some good test > functions to integrate so I can test it further? Preferably something > that has an analytical solution over the reference tet but that's not > crucial.
I've done a bit more extensive testing of the 14-point rule. The attached figure shows the results of integrating the function w^3 sin(wx) sin(wy) sin(wz) over the unit tet with the existing 5th degree rule and Walkington's 14-point rule. In case the attachment doesn't make it through: both rules perform about the same as the frequency w increases, with Walkington's rule generally being slightly better for this particular example. I think it's worth implementing in the library; I'll leave the original 15 point rule there, commented out with a note about why. -- John
comparison.pdf
Description: Adobe PDF document
------------------------------------------------------------------------- Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://sourceforge.net/services/buy/index.php
_______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
