at the moment the defaults for the jacobian_inplace functions are to not use sparsity, which works well for dense matricies, but for sparse you might want to know if you are pre-calculating the sparsity or not. Perhaps the default should fail with notimplemented?