Implement Hessian sparsity idea discussed with coffrin

[odow]

Gravity produces models with >2x fewer hessian non-zeros. I think they're collating terms for each element, rather than passing the hessian vector with repeated elements, like [(1, 1), (1, 1)].