Mark ``no return'' functions for static analysis

This enables the clang static analysis to recognize conditions in
BFAM_ASSERT, BFAM_ABORT_IF, and friends.
Add more floating point type information

This adds floating point type information that will be used to compute
Jacobi polynomials.
Add a function to evaluate Jacobi polynomials
Ignore cast-align warnings from p4est
Make zeroing gopt structure explicit

This avoids a warning from the clang static analyzer.
Make Jacobi comments compile in pandoc

Here we make the comments compile with pandoc.  In the future I would
like to write a script to extract the comments and code and compile
them for documentation.  This could be done by modifying
[shooco](http://rtomayko.github.com/shocco/).
Add approximate equality macros for floats
Add argument checks to Jacobi functions
Fix indent
Fix typo in the comments of the Jacobi header
Add gradient of the orthonormal Jacobi polynomials

This adds a function that evaluates the gradient of the orthonormal
Jacobi polynomials.  Currently this function is not tested to its full
precision (we are using MATLAB to compute double precision test values).
Cleanup Jacobi tests
Add more long real type information
Add tan macro for long real
Add function to compute Jacobi-Gauss quadrature

This initial attempt works for Legendre-Gauss quadrature but does not
work general Jacobi-Gauss quadratures.  A better initial guess is
needed.
Add Jacobi Gauss Lobatto Quadrature function

The values for the test program were generated from Walter Gautschi's
MATLAB code from the paper ``High-order generalized Gauss-Radau and
Gauss-Lobatto formulae for Jacobi and Laguerre weight functions''.
Gautschi's codes can be found at:

  http://www.cs.purdue.edu/archives/2002/wxg/codes/globatto_jacobi.m
  http://www.cs.purdue.edu/archives/2002/wxg/codes/r_jacobi.m
  http://www.cs.purdue.edu/archives/2002/wxg/codes/gauss.m