Romain Veltz

You can already have this feature `br.X` when `X` is in `record_from_solution`

As for NonlinearSolve, I am thinking about this with @Datseris My main blocking point for now is that I want to allow a vector inferface which is not an AbstractArray...

No, define ``` prob = BifurcationProblem(..., record_from_solution = (x,p) -> (X = x[1], Y = x[2])) ```

The biggest advantage would be to tackle parameter derivatives not with finite differences

Or really?? you mean this https://jutho.github.io/KrylovKit.jl/stable/ and [VectorInterface.jl](https://github.com/Jutho/VectorInterface.jl) ?

Should we close this? It is seems that BK is a now an external dep of MTK

Sorry I misread. > Would it make sense to actually store the NonlinearProblem all the way through BifurcationKit's internal process? You could save the `NonlinearProblem` in your own `BifurcationProblem`, it...

I am not sure, I have not designed it specifically. Althoough I dont think there are any globals so it should be fine

Thats what I had in mind!

> Hi, what is the intended geometry ? Rectangle > Can you use a grid aligned with the coordinate directions ? Yes I can