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. 

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...

-Ben 


-------------------------------------------------------------------------
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

Reply via email to