fixed one indexing bug in the hessian computation; added various helper methods for the derivatives of the strains and fd approximations for them

[[Imported from SVN: r1501]]