[[Imported from SVN: r8973]]

[[Imported from SVN: r8972]]

[[Imported from SVN: r8968]]

[[Imported from SVN: r8967]]

[[Imported from SVN: r8966]]

[[Imported from SVN: r8965]]

[[Imported from SVN: r8878]]

This more general method works even when using parametrized
boundaries.

[[Imported from SVN: r8877]]

[[Imported from SVN: r8876]]

[[Imported from SVN: r8871]]

[[Imported from SVN: r8870]]

properly implement getting the Lagrange nodes.  Works without a priori information for all orders and dimensions

[[Imported from SVN: r8869]]

[[Imported from SVN: r8868]]

[[Imported from SVN: r8867]]

[[Imported from SVN: r8866]]

[[Imported from SVN: r8865]]

[[Imported from SVN: r8864]]

[[Imported from SVN: r8863]]

[[Imported from SVN: r8862]]

[[Imported from SVN: r8861]]

[[Imported from SVN: r8687]]

[[Imported from SVN: r8638]]

[[Imported from SVN: r8637]]

[[Imported from SVN: r8636]]

[[Imported from SVN: r8635]]

[[Imported from SVN: r8634]]

[[Imported from SVN: r8633]]

[[Imported from SVN: r8632]]

[[Imported from SVN: r8631]]

[[Imported from SVN: r8630]]

[[Imported from SVN: r8629]]

[[Imported from SVN: r8628]]

[[Imported from SVN: r8627]]

There were two issues fixed by this patch:

a) The bending energy only involves the third director.
   One must look at it as a function in S^2, not as
   a subset of coefficients of a function to SO(3).
   This means that we have to use a different projection
   in one place.
b) A certain subset of the coefficients (see the code)
   must always be zero, because they are parts of
   3-vectors which are known to be perpendicular to (0,0,1).
   For some strange reason they aren't.  I now hardwire
   them to be 0 anyway.  This makes my test pass.
   Maybe one day I'll find the true reason :-)

[[Imported from SVN: r8626]]

of Neumann boundary forces.  But I am not sure and I want to
postpone the issue.  Therefore this patch adds a cpp warning.

[[Imported from SVN: r8625]]

[[Imported from SVN: r8597]]

Don't use a loop with 4 iteration to clear an array,
if the array has only length 3.

Boy, this did take a while to find...

[[Imported from SVN: r8593]]

Another easy patch which makes the code a lot faster.

Also, this patch removes the quasi-Newton code.  So far that code
hasn't helped at all.  Quasi-Newton methods will need more work
in the future.

[[Imported from SVN: r8588]]

the ones with i<=j.  The rest follows from symmetry of
the Hessian.

Needless to say, this makes the Hessian assembly quite
a bit faster.

[[Imported from SVN: r8587]]

space of orthogonal matrices.

[[Imported from SVN: r8573]]