Add a working nondiagonalized neo Hookean constitutive model
Nonlinear elasticity half works

There are problems with divergent line searches, which I plan to fix
by switching to TAO (since we're doing optimization).
SConstruct: Sync with geode/SConstruct
Use TAO for elastic optimization