Nick

> Why would we exclude the case of roots that touch but do not cross the origin? As far as I know, TOMS 748 imposes no such restriction. It's not...

For those finding this issue, I have managed to figure out how to get this working: Simply add ```julia Zygote.@adjoint function Quantity(x::Number, u) return Quantity(x, u), (∂y -> (∂y *...

The numerical integrators already work this way. For example the tanh-sinh integrator gets all the nodes/weights created in the constructor, and then many integrands can be integrated using the `.integrate`...

Patch file which implements this feature: [unwrap.patch](https://github.com/user-attachments/files/23915712/unwrap.patch)