utils.c

Go to the documentation of this file.

/*
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#include "utils.h"
#include <limits.h>
#include <stdbool.h>
#include <time.h>
 
void
rand_init(void)
{
    time_t now = time(0);
    const unsigned char *p = (unsigned char *) &now;
    uint32_t seed = 0;
    for (size_t i = 0; i < sizeof(now); ++i) {
        seed = (seed * (UCHAR_MAX + 2U)) + p[i];
    }
    dsfmt_gv_init_gen_rand(seed);
}
 
void
rand_init_seed(const uint32_t seed)
{
    dsfmt_gv_init_gen_rand(seed);
}
 
double
rand_uniform(const double min, const double max)
{
    return min + (dsfmt_gv_genrand_open_open() * (max - min));
}
 
int
rand_uniform_int(const int min, const int max)
{
    return (int) floor(rand_uniform(min, max));
}
 
double
rand_normal(const double mu, const double sigma)
{
    static const double two_pi = 2 * M_PI;
    static double z1;
    static bool generate;
    generate = !generate;
    if (!generate) {
        return z1 * sigma + mu;
    }
    const double u1 = dsfmt_gv_genrand_open_open();
    const double u2 = dsfmt_gv_genrand_open_open();
    const double z0 = sqrt(-2 * log(u1)) * cos(two_pi * u2);
    z1 = sqrt(-2 * log(u1)) * sin(two_pi * u2);
    return z0 * sigma + mu;
}
 
void
utils_json_parse_check(const cJSON *json)
{
    if (json == NULL) {
        printf("Error reading JSON\n");
        const char *error_ptr = cJSON_GetErrorPtr();
        if (error_ptr != NULL) {
            printf("Error before: %s\n", error_ptr);
        }
        exit(EXIT_FAILURE);
    }
}