On Tue, Jun 17, 2008 at 2:11 PM, Benjamin Kirk <[EMAIL PROTECTED]> wrote: > Have we settled in on CPPUNIT as our unit-testing framework? I'm ready to > reorganize the unit test directory if so. > > On a totally related issue, I was just expanding the quadrature unit tests > and found a bug in QGauss for Tets at 5th-order. The issue is that the > points/weights vector is resized to hold 17 entries, but only 15 entries are > initialized. This led to weirdness which can allow the last two entries of > the weights vector to be nonzero depending on what you had previously > computed on.
That's a nasty one. I would never have seen that -- I'd be to focused on making sure the weights and points were correct to worry if the vector was too long! > Do any of you have access to "The finite element method" vol. 1 by > Zienkiewicz & Taylor? Pg. 222 is referenced as the source for orders 1-5, > and I want to make sure that the 5th-order rule has only 15 points... I found one of Zienkiewicz's books on google books The Finite Element Method: Its Basis and Fundamentals Redacted from the book preview is pg. 167, which should contain Table 5.4 -- tetrahedral quadrature formulae... this always seems to happen to me with Google Books. Maybe Roy can send Vikram to get one of the many copies from the UT Engineering Library: TA 640.2 Z5 -- John ------------------------------------------------------------------------- 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
