clairvoyante

https://software.intel.com/content/www/us/en/develop/articles/the-difference-between-x87-instructions-and-mathematical-functions.html http://notabs.org/fpuaccuracy/index.htm https://github.com/JuliaMath/openlibm/issues/11

exp(x) use Padé approximant. Same order Padé approximant of exp(x)-1 is: > (312*x^11+1201200*x^9+1176215040*x^7+402265543680*x^5+46930980096000*x^3+1295295050649600*x)/(x^12-156*x^11+12012*x^10-600600*x^9+21621600*x^8-588107520*x^7+12350257920*x^6-201132771840*x^5+2514159648000*x^4-23465490048000*x^3+154872234316800*x^2-647647525324800*x+1295295050649600) (cf. https://sagecell.sagemath.org/?z=eJwrSExJ1ShJrMzJL9JIrSjQqNDUNdSp0DHQMTLR1DE0AiJNawDbhwo3&lang=maxima&interacts=eJyLjgUAARUAuQ== ) So... ``` def exp(self): n = int(round(self.x)) x = self - n x2 =...