In the mixed-mode Jacobian example in #684, how can we specify the known sparsity pattern in the backend construction? And the related documentation on mixed mode Jacobian is really scarce, maybe we can put the benchmarking example as documentation to illustrate how to exploit MixedMode Jacobian to achieve faster Jacobian computation.
(Open a new issue to not disrupt #684)