note to self to try this and see if `phase_diagram` can be sped up further, especially for `IntegralModels`: https://github.com/gdalle/ImplicitDifferentiation.jl