DOK matrix is not adquate: does not store accidental zeros, which leads
to problems when crosschecking Jacobian