arb
arb copied to clipboard
flintlib
Arb has been merged into FLINT -- use https://github.com/flintlib/flint/ instead
Would it make sense to add functions like conj, real_part, imag_part, and abs for acb_poly? My use case is that I have an acb->acb function and I am using acb_poly...
An `acb_mat` type matrix (say 100x100 cells) is the output of a summation done in a long (n=1..10^12) loop. To speed up the processing, the loop is split up in...
I was wondering if it was possible to perform rigorous integration using Arb with integrands which do not have a closed-form expression (i.e. that are obtained numerically for only a...
Hi Fredrik, I am one of the developers of the [Anima](https://github.com/Inria-Visages/Anima-Public) library at INRIA and maintained by @ocommowi. This is a C++11 library for medical image processing. It occurs relatively...
The [Riemann theta function](https://en.wikipedia.org/wiki/Theta_function#Riemann_theta_function) is a natural generalization of the Jacobi theta functions. In this more general situation, the invariance involves the integral symplectic group acting on the Siegel half...
``` AF = ArbField(123) posinf = AF("Inf"); neginf = -posinf; posinf > neginf, posinf, neginf (false, true, [+/- inf], [+/- inf]) altinf = AF(1)//AF(0) [+/- inf] altinf == posinf, altinf...
Similarly to [flint issue #447](https://github.com/wbhart/flint2/issues/447) arb fails to build in Sage when make is compiled with libguile support (that happens to be the case on archlinux). As a consequence, modifying...
It's not trivial to parallelize the main work queue, but the function evaluations in the Gaussian quadrature could be parallelized, including the computation of quadrature nodes. It would definitely be...
The integration code could be extended to allow computing f_0(x), ..., f_n(x) simultaneously. It's often useful to compute a sequence of integrals where most of the integrand evaluation can be...
Consider the program: ```cpp #include "arb.h" #include int main() { arb_t x; arb_init(x); arb_add_error_2exp_si(x, -30); // Add 2^(-30) to radius of x. arb_t y; arb_init(y); // y = 0, exactly....
