Le 04/05/2011 19:01, Joris Vankerschaver a écrit :
Do you mean that it is possible to define the RHS as a Cython *callback* function? or is there an other trick ? Can you give me a pointer to that ?On May 3, 11:23 pm, Thierry Dumont<tdum...@math.univ-lyon1.fr> wrote:1) it takes time to do this, 2) we have to solve the callback problem: such program make a lot of callbacks (to the rhs of the system): AFAIK, there are no simple method to make this work fast enough (the ODE solver in Sage is extremely slow, too).Is it really that slow? The ode_solver class, for instance, allows you to define the RHS as a cython function, and behind the scenes a pointer to the RHS is then passed around as it were any other compiled function. I don't think this would be too much slower than a generic C implementation. I haven't done any speed tests, but I always find the ode_solver class to be very fast. J.
Some time ago, I asked here how to define Cython callbacks, and the answer was that it was impossible.
If it is possible, I propose (me) to make a package of ODE/DAE intergrators: ther is a lot of methods which would be nice in Sage like Implicit R-K methods, Symplectic methods and so on...
t.d. -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
<<attachment: tdumont.vcf>>
