Romain Veltz

> In addition I think that all constitutive functions could be like c(x,u,lambda) where lambda is a vector of parameters, and I intend to provide AD based \partial u /\partial...

Hi, Did you progress on this by any chance?

> But I took a couple of hours now to figure out some possibilities. Wow, thank you! > How would you find something like this? You "improved" the specification of...

The continuation steps through it most likely. This solver is well described here: Keller, Herbert B. “The Bordering Algorithm and Path Following Near Singular Points of Higher Nullity.” SIAM Journal...

yes that's the jvp

I thought DE would provide a built in solution like DiffEqOperator used to.

👍

> Let's standardize our notation. If I understand you well, `L(u,p,t)` is a linear operator representing the differential of say `F(u,p,t)` wrt to the first variable at position `u`.

How is the progress on this issue?

Sorry to be pushy... did you make any progress?