I need to provide a function to a c++ library and I want the user to pass a Julia function. The c++ function is like
f(dy,t,y) and is called all over the place in the c++ code. I have succeeded in the past to write closures when the function definition is
f(dy,t,y,param) by copying the method used in Sundials.
However I have no clue how to deal with this situation.
Does anyone have an idea?
Thank you very much,