(Daniel Wirtz, 2013-10-24) Providing a default sparse identity matrix as dydxp-jacobian (i.e. matrix M and M-sparsity pattern). This yields a major increase in performance when providing f-Jacobian information. This is due to a ode15i sparse Miter matrix, which is dense if odenumjac is used to compute the dydxp value/pattern).