arb
arb copied to clipboard
flintlib
Arb has been merged into FLINT -- use https://github.com/flintlib/flint/ instead
With the new algorithm to evaluate R_J, there are some regressions in elliptic_pi_inc (see https://trac.sagemath.org/ticket/29719): * Slightly lower accuracy * NaN result when evaluated at inexact phi = pi/2 (this...
Calling acb_elliptic_pi_inc(a_res, a_n, a_phi, a_m, quasiperiod_1, 52) like this (Mathematica syntax) ``` EllipticPi[+0.000000 -3.000000 Sqrt[-1], +0.769972 +0.680250 Sqrt[-1], -9.000000 -0.000000 Sqrt[-1]] = +0.575274 -0.173587 i ``` but Wolfram Mathematica states...
Doing `./configure && make && sudo make install` on a Linux system results in `librarb.so*` installed in `/usr/local/lib64`, but not known to the loader, as `ldconfig` is not run. IMHO...
"make" is okay but not for "make check". flint has successfully passed "make" and "make check". Here is the error message ``` partitions_fmpz_ui_threaded....warning: newton_refine_root: improvement failed make[1]: *** [../build/partitions/test/t-partitions_fmpz_ui_threaded_RUN] Segmentation...
Would there be any interest in switching the build system over to a standard autoconf setup? I might have time to work on this…
I'd like to integrate a function over cubical region in CC^n (n-fold integral). At the moment I define integrand recursively, but of course i'd be much faster to compute multi-dimensional...
I believe the result should be complex infinity.
We want to have good support for Arb in CSymPy, essentially allow users to use Arb to evaluate expressions, and it will give the result with a given accuracy. What...
The below (and attached) code arb_one(a); arb_set_str(b,"[-0.236000000000000000000000000000000 +/- 5.12e-34]",prec); arb_set_str(z,"[69.314718055994530941723212145817656807550013436025525 +/- 4.13e-52]",prec); arb_indeterminate(res); arb_hypgeom_u(res,a,b,z,prec); arb_printn(res,53,0); printf("\n"); arb_set_str(b,"-0.236000000000000000000000000000000",prec); arb_set_str(z,"69.314718055994530941723212145817656807550013436025525",prec); arb_indeterminate(res); arb_hypgeom_u(res,a,b,z,prec); arb_printn(res,53,0); printf("\n"); gives: [+/- 30.4] ## with [0.0139820478288060060036870204487 +/- 7.41e-32]...
