Modifications have been made to correct behavior of h-adaptivity with higher order SZABAB and BERNSTEIN:

--  defined SZABAB to be hierarchic

-- defined compute_constraints function for SZABAB and BERNSTEIN

-- defined extra_hanging_dofs for SZABAB and BERNSTEIN to be true

-- FEInterface::compute_constraints for SZABAB and BERNSTEIN now call the compute constraint functions for the respective elements.
forward declare Elem
use openmp when enabled
just pass vectors through to the PETSc interface
update download link
Avoid narrowing warnings when initializing arrays.
Bugfix for System::point_*
Merge pull request #183 from friedmud/petsc_vector_optimize

just pass vectors through to the PETSc interface
Make libmesh-config --ldflags return whatever libmesh_LDFLAGS was set
to during configure.
when using openmp set omp_num_threads
use static scheduling and remove unnecessary omp_get_thread_num
Merge pull request #182 from friedmud/openmp_merge

use openmp when enabled
guard the new add_vector implementation properly
rbOOmit change: Fixed bugs in RBEIMEvaluation::read_in_interpolation_points_elem
rbOOmit change: Removed inverse_map from RBEIMAssembly::evaluate_basis_function. Should be a significant optimization for large systems.
rbOOmit change: Added logging in RBEIMAssembly::evaluate_basis_function
rbOOmit change: Further optimized RBEIMAssembly::evaluate_basis_function by storing FE objects. This means we benefit from caching in the calls to reinit.
Updating citations:
* Fuentes, 2013b
* Rong, 2013
* Fixing Dodd III
* Tonks, 2014
Allow TimeSolver tolerances relative to first step

Specifying a negative tolerance now lets you avoid trying to directly
calculate or guess what a decent L2/H1/whatever error growth rate is,
by using a multiple of the growth rate observed at the initial time
step.
Better AdaptiveTimeSolver documentation