Output 2nd derivatives of static model when doing identification

This is necessary for global identification.

The code for 2nd static derivatives already exists (used for planner_objective).