> then your changesets will be less invasive
> and more easily accepted, do you think that
> would be possible?
>
> c.
Hi Carlo!
I wrote the class without inheriting from DAE and avoiding global pointers for residual function and jacobian evaluation.
You can see the new design of the dld-function [1] and the class [2][3]. Your suggestion really helped.
I'm going to complete the options of the solver. I've already written some test here [4]. In the m-interface I use a private function to check odefun, tspan and init [5]. I think we could adapt it also for ode23 and ode45.