[[Imported from SVN: r2017]]

make sure the fd hessian is symmetric; hardwire Dexp(0) for more speed; use symmetry when computing the differential of the geodesic interpolation

[[Imported from SVN: r2016]]

[[Imported from SVN: r2015]]

add a new method which averages out a dg function in order to obtain a p1 function.  This is used to pass Neumann data coming from a rod to the hierarchic error estimator

[[Imported from SVN: r1866]]

[[Imported from SVN: r1865]]

[[Imported from SVN: r1792]]

[[Imported from SVN: r1791]]

[[Imported from SVN: r1790]]

[[Imported from SVN: r1775]]

[[Imported from SVN: r1769]]

[[Imported from SVN: r1768]]

[[Imported from SVN: r1767]]

[[Imported from SVN: r1766]]

[[Imported from SVN: r1719]]

[[Imported from SVN: r1718]]

[[Imported from SVN: r1716]]

[[Imported from SVN: r1715]]

moved the quadrature orders for bending and shear stuff into the header, as a first step towards reduced integration

[[Imported from SVN: r1714]]

[[Imported from SVN: r1713]]

[[Imported from SVN: r1712]]

[[Imported from SVN: r1711]]

[[Imported from SVN: r1710]]

[[Imported from SVN: r1709]]

[[Imported from SVN: r1708]]

[[Imported from SVN: r1703]]

[[Imported from SVN: r1702]]

[[Imported from SVN: r1701]]

[[Imported from SVN: r1700]]

[[Imported from SVN: r1699]]

[[Imported from SVN: r1698]]

[[Imported from SVN: r1696]]

[[Imported from SVN: r1680]]

print max translation correction and max rotation correction separately.  Slightly fiddling with the acceptance criterion

[[Imported from SVN: r1679]]

[[Imported from SVN: r1678]]

[[Imported from SVN: r1677]]

[[Imported from SVN: r1658]]

[[Imported from SVN: r1652]]

[[Imported from SVN: r1651]]

[[Imported from SVN: r1650]]

[[Imported from SVN: r1638]]