Skip to content

Commit

Permalink
Fix precision in sum_radius
Browse files Browse the repository at this point in the history
  • Loading branch information
j-kieffer committed Feb 4, 2025
1 parent ecee23c commit 9ccc9a3
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions src/acb_theta/sum_radius.c
Original file line number Diff line number Diff line change
Expand Up @@ -89,9 +89,9 @@ acb_theta_sum_radius(arf_t R2, arf_t eps, const arb_mat_t cho, slong ord, slong
arb_init(sqrt2pi);
arf_init(cmp);

arb_const_pi(sqrt2pi, prec);
arb_const_pi(sqrt2pi, lp);
arb_mul_2exp_si(sqrt2pi, sqrt2pi, 1);
arb_sqrt(sqrt2pi, sqrt2pi, prec);
arb_sqrt(sqrt2pi, sqrt2pi, lp);

/* Set b such that
(1 + 8/sqrt(pi)) * prod_j (1 + sqrt(2pi)/c_j) * b \leq 2^(-prec) */
Expand Down

0 comments on commit 9ccc9a3

Please sign in to comment.