hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-04 17:01:40 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

heuristic parameter values for greater heights

Browse source
This commit is contained in:
p15-git-acc 2020-09-20 11:00:00 -05:00
parent dd50190726
commit fafce0f7b3
1 changed files with 9 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
acb_dirichlet/platt_local_hardy_z_zeros.c
View file

@ -1373,7 +1373,7 @@ _create_heuristic_context(const fmpz_t n, slong prec)
{
goto finish;
}
else if (_fmpz_cmp_a_10exp_b(n, 1, 7) < 0)
else if (_fmpz_cmp_a_10exp_b(n, 1, 7) <= 0)
{
A = 8;
B = 4096;
@ -1385,7 +1385,7 @@ _create_heuristic_context(const fmpz_t n, slong prec)
dh = 189.300555621 + -18.9197620910*x + 0.830573279739*x*x;
dH = 0.996344708432 + 0.00726343756987*x + 0.000214553523479*x*x;
}
else if (_fmpz_cmp_a_10exp_b(n, 1, 17) < 0)
else if (_fmpz_cmp_a_10exp_b(n, 1, 17) <= 0)
{
A = 8;
B = 4096;
@ -1397,17 +1397,17 @@ _create_heuristic_context(const fmpz_t n, slong prec)
dh = 94.8311227685 + -0.883336123977*x + 0.0556261501085*x*x;
dH = 1.18099035487 + -0.0119992923413*x + 0.00073741022749*x*x;
}
else if (_fmpz_cmp_a_10exp_b(n, 1, 22) < 0)
else if (_fmpz_cmp_a_10exp_b(n, 1, 37) <= 0)
{
A = 16;
B = 8192;
Ns_max = 300;
sigma_interp = 19;
dlogJ = -0.506132895307 + 0.497565124541*x + 2.25705259466e-06*x*x;
dK = 103.630048145 + -0.867460371390*x + 0.00106534037344*x*x;
dgrid = 4085.55616952 + -1.36695690825*x + -0.0198932963863*x*x;
dh = 142.523364391 + -0.057841628146*x + 0.0157997314814*x*x;
dH = 0.727682145349 + -0.00353076385259*x + 0.000111471602416*x*x;
sigma_interp = _fmpz_cmp_a_10exp_b(n, 1, 27) < 0 ? 19 : 17;
dlogJ = -0.510731356439 + 0.497749405271*x + 4.7480897930e-07*x*x;
dK = 98.3296668968 + -0.629042293257*x + -0.00158584251152*x*x;
dgrid = 3810.92686930 + 9.86688705118*x + -0.132241869233*x*x;
dh = 189.640536974 + -1.94782986530*x + 0.0342403597641*x*x;
dH = 0.897872657566 + -0.0100363606454*x + 0.000170869048778*x*x;
}
else
{