[sage-devel] Re: Automatic Differentiation - Taylor Arithmetics

2019-05-17 Thread Nisoli Isaia
: > > Nisoli Isaia wrote: > > 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 qu

[sage-devel] Re: Automatic Differentiation - Taylor Arithmetics

2019-04-23 Thread Nisoli Isaia
Thank you! I think I could start by exposing these ARB functions to SAGE, so that the interface is more complete, and then try to do the same for Real Interval. Best wishes Isaia On Monday, April 22, 2019 at 4:59:36 AM UTC-3, Marc Mezzarobba wrote: > > Nisoli Isaia wrote: > > P.=

[sage-devel] Re: Automatic Differentiation - Taylor Arithmetics

2019-04-20 Thread Nisoli Isaia
On Saturday, April 20, 2019 at 5:36:46 AM UTC-3, Marc Mezzarobba wrote: > > Nisoli Isaia wrote: > > I was planning in doing a Cython implementation of Forward automatic > > differentiation and > > Taylor arithmetics as in > > https://press.princeton.ed

[sage-devel] Automatic Differentiation - Taylor Arithmetics

2019-04-19 Thread Nisoli Isaia
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

[sage-devel] cocalc image, cython and cypari.h

2017-11-06 Thread Nisoli Isaia
Dear all, after upgrading to the last docker image of cocalc, cython complains with a missing "cypari.h" header. _projects_38d5792a_e2ea_4974_a00f_8fe2468b2337_dynamic_cython_spyx_1.c:441:20: fatal error: cypari.h: No such file or directory #include "cypari.h" Before cython worked with load()

[sage-devel] Interval Sparse/Dense Matrices

2016-05-25 Thread Nisoli Isaia
I'm thinking in implementing cython classes for mpfr_t and mpfi_t dense/sparse matrices, imitating the dense/sparse matrices for the mpq_t type. Would it be of interest for the community? Is there any previous work on the issue, ongoing discussion on the topic? Best wishes Isaia -- You

[sage-devel] matrix vector product for generic sparse matrices

2016-05-25 Thread Nisoli Isaia
I opened a ticket due to the following issue: It seems like this simple cython code is faster 30% faster than the implemented matrix vector product for generic sparse matrices. I run some (many) tests with matrices with RIF coefficients and all seem to point in this direction. Is it possible