You might be interested in taking a look at tides [1]. It is a library for integration of ODE's in high precission using precisely automatic differentiation and Taylor arithmetic. It doesn't use interval arithemtic (so, no certified computation), but maybe some ideas could be useful. The way it works is by translating the function that defines your ODE to c-code and then compiling it.
There is a tides optional package for sage, and an interface to it. [1] https://iuma.unizar.es/es/investigacion/software/TIDES El viernes, 19 de abril de 2019, 14:48:09 (UTC+2), Nisoli Isaia escribió: > > Dear all, > I was planning in doing a Cython implementation of Forward automatic > differentiation and > Taylor arithmetics as in > https://press.princeton.edu/titles/9488.html > to use to implement a library for Sage with rigorous quadrature and > integration of ODE. > > I'm trying to understand which parent class could be the best for these > objects for them > to be compatible with Sage coercion model. > > What is the usual procedure, I read through > > http://doc.sagemath.org/html/en/reference/structure/sage/structure/element.html > but I'm still confused. > > Best wishes > Isaia > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-devel@googlegroups.com. Visit this group at https://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.