@@ -433,7 +433,7 @@ irk_pareto_vec(irk_state *state, npy_intp len, double *res, const double alp)
433
433
void
434
434
irk_weibull_vec (irk_state *state, npy_intp len, double *res, const double alp)
435
435
{
436
- int i, err;
436
+ int err;
437
437
const double d_zero = 0.0 , d_one = 1.0 ;
438
438
double rec_alp = 1.0 /alp;
439
439
@@ -458,7 +458,7 @@ irk_weibull_vec(irk_state *state, npy_intp len, double *res, const double alp)
458
458
void
459
459
irk_power_vec (irk_state *state, npy_intp len, double *res, const double alp)
460
460
{
461
- int i, err;
461
+ int err;
462
462
const double d_zero = 0.0 , d_one = 1.0 ;
463
463
double rec_alp = 1.0 /alp;
464
464
534
534
irk_f_vec (irk_state *state, npy_intp len, double *res, const double df_num, const double df_den)
535
535
{
536
536
int err;
537
- const double d_zero = 0.0 , d_one = 1.0 ;
537
+ const double d_zero = 0.0 ;
538
538
double shape = 0.5 *df_num, scale = 2.0 /df_num;
539
539
double *den = NULL ;
540
540
@@ -642,7 +642,6 @@ irk_noncentral_chisquare_vec(irk_state *state, npy_intp len, double *res, const
642
642
for (i = 0 ; i < len; ) {
643
643
int k, j, cv = pvec[idx[i]];
644
644
645
- DIST_PRAGMA_VECTOR
646
645
for (j=i+1 ; (j < len) && (pvec[idx[j]] == cv); j++) {}
647
646
648
647
assert (j > i);
728
727
irk_logistic_vec (irk_state *state, npy_intp len, double *res, const double loc, const double scale)
729
728
{
730
729
int i, err;
731
- const double d_one = 1.0 , d_mone = - 1.0 , d_zero = 0.0 ;
730
+ const double d_one = 1.0 , d_zero = 0.0 ;
732
731
733
732
if (len < 1 )
734
733
return ;
@@ -862,7 +861,7 @@ irk_wald_vec(irk_state *state, npy_intp len, double *res, const double mean, con
862
861
static void
863
862
irk_vonmises_vec_small_kappa (irk_state *state, npy_intp len, double *res, const double mu, const double kappa)
864
863
{
865
- int i, err, n, size, blen = 1 << 20 ;
864
+ int i, err, n, size;
866
865
double rho_over_kappa, rho, r, s_kappa, Z, W, Y, V;
867
866
double *Uvec = NULL , *Vvec = NULL ;
868
867
float *VFvec = NULL ;
@@ -923,9 +922,9 @@ irk_vonmises_vec_small_kappa(irk_state *state, npy_intp len, double *res, const
923
922
static void
924
923
irk_vonmises_vec_large_kappa (irk_state *state, npy_intp len, double *res, const double mu, const double kappa)
925
924
{
926
- int i, err, n, size, blen = 1 << 20 ;
925
+ int i, err, n, size;
927
926
double r_over_two_kappa, recip_two_kappa;
928
- double s_minus_one, hpt, r_over_two_kappa_minus_one, rho_minus_one, neg_W_minus_one ;
927
+ double s_minus_one, hpt, r_over_two_kappa_minus_one, rho_minus_one;
929
928
double *Uvec = NULL , *Vvec = NULL ;
930
929
float *VFvec = NULL ;
931
930
const double d_zero = 0.0 , d_one = 1.0 ;
@@ -1000,8 +999,6 @@ irk_vonmises_vec_large_kappa(irk_state *state, npy_intp len, double *res, const
1000
999
void
1001
1000
irk_vonmises_vec (irk_state *state, npy_intp len, double *res, const double mu, const double kappa)
1002
1001
{
1003
- int blen = 1 << 20 ;
1004
-
1005
1002
if (len < 1 )
1006
1003
return ;
1007
1004
@@ -1022,7 +1019,7 @@ irk_vonmises_vec(irk_state *state, npy_intp len, double *res, const double mu, c
1022
1019
void
1023
1020
irk_noncentral_f_vec (irk_state *state, npy_intp len, double *res, const double df_num, const double df_den, const double nonc)
1024
1021
{
1025
- int i, blen = 1 << 20 ;
1022
+ int i;
1026
1023
double *den = NULL , fctr;
1027
1024
1028
1025
if (len < 1 )
@@ -1281,7 +1278,7 @@ irk_poisson_vec_POISNORM(irk_state *state, npy_intp len, int *res, const double
1281
1278
void
1282
1279
irk_poisson_vec_V (irk_state *state, npy_intp len, int *res, double *lambdas)
1283
1280
{
1284
- int err, blen ;
1281
+ int err;
1285
1282
1286
1283
if (len < 1 )
1287
1284
return ;
@@ -1304,7 +1301,7 @@ irk_poisson_vec_V(irk_state *state, npy_intp len, int *res, double *lambdas)
1304
1301
void
1305
1302
irk_zipf_long_vec (irk_state *state, npy_intp len, long *res, const double a)
1306
1303
{
1307
- int i, err, n_accepted, batch_size, blen = 1 << 20 ;
1304
+ int i, err, n_accepted, batch_size;
1308
1305
double T, U, V, am1, b;
1309
1306
double *Uvec = NULL , *Vvec = NULL ;
1310
1307
long X;
@@ -1897,7 +1894,6 @@ void
1897
1894
irk_rand_int64_vec (irk_state *state, npy_intp len, npy_int64 *res, const npy_int64 lo, const npy_int64 hi)
1898
1895
{
1899
1896
npy_uint64 rng;
1900
- int err;
1901
1897
npy_intp i;
1902
1898
1903
1899
if (len < 1 )
0 commit comments