diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/graph_genfunc.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/graph_genfunc.c b/src/graph_genfunc.c index 39a8480..54d9daa 100644 --- a/src/graph_genfunc.c +++ b/src/graph_genfunc.c @@ -2,7 +2,7 @@ #include "fracture.h" double *genfunc_uniform(unsigned int L, bound_t boundary, gsl_rng *r, unsigned int *num) { - *num = pow(L / 2 + 1, 2) + pow((L + 1) / 2, 2); + *num = 2 * pow(L / 2, 2); double *lattice = (double *)malloc(2 * (*num) * sizeof(double)); for (unsigned int i = 0; i < (*num); i++) { @@ -18,7 +18,7 @@ double g(double rho, double dist) { } double *genfunc_hyperuniform(unsigned int L, bound_t boundary, gsl_rng *r, unsigned int *num) { - *num = pow(L / 2 + 1, 2) + pow((L + 1) / 2, 2); + *num = 2 * pow(L / 2, 2); // necessary to prevent crashing when underflow occurs gsl_set_error_handler_off(); |