PDE constrained optimization. I was optimizing for many hundreds of thousands of parameters. I only needed the diagonal of the Jacobian by the end of the calculation, however it was necessary to build the full dense matrix first without rewriting massive parts of the logic.