Fix bug in gasdev

milestone_5_without_improvements-logs
chortas 3 years ago
parent f43eec7034
commit ec130dddd5

@ -422,7 +422,7 @@ float gammp(float a, float x);
/*and unit variance, using ran1(idum) as the source */ /*and unit variance, using ran1(idum) as the source */
/*of uniform deviates */ /*of uniform deviates */
/*idum: seed */ /*idum: seed */
double gasdev(long* idum, long* idum2, long* iy, long* iv, int* iset); double gasdev(long* idum, long* idum2, long* iy, long* iv);
/*gaussian covariance value for lag h*/ /*gaussian covariance value for lag h*/
double gaussian(double h); double gaussian(double h);

@ -49,7 +49,7 @@ void generate(long* seed, int n, struct realization_mod* realization, int cores)
/*Gaussian white noise generation*/ /*Gaussian white noise generation*/
for (i = 0; i < n; i++) for (i = 0; i < n; i++)
(*realization).vector[i] = gasdev(seed, &idum2, &iy, iv, &iset); (*realization).vector[i] = gasdev(seed, &idum2, &iy, iv);
t = clock() - t; t = clock() - t;
double time_taken = ((double)t)/CLOCKS_PER_SEC; // calculate the elapsed time double time_taken = ((double)t)/CLOCKS_PER_SEC; // calculate the elapsed time

@ -423,7 +423,7 @@ float gammp(float a, float x);
/*and unit variance, using ran1(idum) as the source */ /*and unit variance, using ran1(idum) as the source */
/*of uniform deviates */ /*of uniform deviates */
/*idum: seed */ /*idum: seed */
double gasdev(long* idum, long* idum2, long* iy, long* iv, int* iset); double gasdev(long* idum, long* idum2, long* iy, long* iv);
/*gaussian covariance value for lag h*/ /*gaussian covariance value for lag h*/
double gaussian(double h); double gaussian(double h);

@ -150,8 +150,9 @@ void log_log(int level, const char *file, int line, const char *fmt, ...) {
char* env_var = getenv("ENV"); char* env_var = getenv("ENV");
if (env_var != NULL && strcmp("false", env_var) == 0) return; if (env_var != NULL && strcmp("false", env_var) == 0) return;
char* substr = strstr(fmt, "MEM"); char* substr_mem = strstr(fmt, "MEM");
if (env_var != NULL && strcmp("analysis", env_var) == 0 && substr == NULL) return; char* substr_cpu = strstr(fmt, "CPU");
if (env_var != NULL && strcmp("analysis", env_var) == 0 && substr_mem == NULL && substr_cpu == NULL) return;
lock(); lock();

Loading…
Cancel
Save