On Tue, 28 Feb 2012, David Knezevic wrote: > I just checked in eigenvalue_ex3. It's an example I made last semester > when I was teaching about eigenvalue problems, and I thought it'd make a > good libMesh example.
It's a great example, when SLEPc's enabled. It's a compilation failure otherwise. ;-) You need #ifdef LIBMESH_ENABLE_SLEPC wrappers around every use of EigenSystem, not just the ones in main(). > Also, it's a good benchmark problem, because the eigenvalues have been > tabulated to high precision in the literature (as indicated in the > example comments). If a benchmark fails in the forest of examples, and nobody's there to evaluate error, does it make a sound? ;-) We really need some benchmarks that get run and tested for excessive error levels as part of the basic open source libMesh suite. "Secret apps that INL runs" plus "Secret apps that UT runs" plus etcetera is much better than nothing, but in light of the recent system-specific (one compiler-specific!) failures of doxygen and xdr on libmesh-users, it would be nice to have a "make check" target that catches numerical errors rather than just logic errors. --- Roy ------------------------------------------------------------------------------ Virtualization & Cloud Management Using Capacity Planning Cloud computing makes use of virtualization - but cloud computing also focuses on allowing computing to be delivered as a service. http://www.accelacomm.com/jaw/sfnl/114/51521223/ _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
