mirror of
https://github.com/vale981/arb
synced 2025-03-05 09:21:38 -05:00
commit
ef5ab30492
1 changed files with 3 additions and 3 deletions
|
@ -340,16 +340,16 @@ stieltjes_mag(double n)
|
||||||
return t;
|
return t;
|
||||||
}
|
}
|
||||||
|
|
||||||
static double __hypot(double x, double y) { return sqrt(x * x + y * y); }
|
static double _d_approx_hypot(double x, double y) { return sqrt(x * x + y * y); }
|
||||||
|
|
||||||
/* log2 magnitude of integrand at z = x+yi; alpha = a+bi */
|
/* log2 magnitude of integrand at z = x+yi; alpha = a+bi */
|
||||||
static double
|
static double
|
||||||
integrand_mag(double n, double x, double y, double a, double b)
|
integrand_mag(double n, double x, double y, double a, double b)
|
||||||
{
|
{
|
||||||
double t, u;
|
double t, u;
|
||||||
t = log(__hypot(a - y, b + x));
|
t = log(_d_approx_hypot(a - y, b + x));
|
||||||
u = atan2(b + x, a - y);
|
u = atan2(b + x, a - y);
|
||||||
t = log(__hypot(t,u)) * (n+1) - 2.0 * 3.1415926535897932 * x;
|
t = log(_d_approx_hypot(t,u)) * (n+1) - 2.0 * 3.1415926535897932 * x;
|
||||||
return t * 1.44269504088896341;
|
return t * 1.44269504088896341;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
Loading…
Add table
Reference in a new issue