Constitutive model tweaks

switched to doing a bit more precomputation for calculating R
differentials, so I don't have to do matrix inverses over and over
again. forceDifferentialDensity is now called dEdFDifferential, because
that's what it actually IS, and that and dEnergyDensitydF() are now
compared with their finite difference equivalents in the snow
constitutive model test