squiggle.c

Self-contained Monte Carlo estimation in C99
Log | Files | Refs | README

example.c (876B)

      1 #include "../../../squiggle.h"
      2 #include <stdio.h>
      3 #include <stdlib.h>
      4 
      5 int main()
      6 {
      7     // set randomness seed
      8     uint64_t* seed = malloc(sizeof(uint64_t));
      9     *seed = 1000; // xorshift can't start with 0
     10 
     11     int n = 1000 * 1000;
     12     double* gamma_array = malloc(sizeof(double) * (size_t)n);
     13     for (int i = 0; i < n; i++) {
     14         gamma_array[i] = sample_gamma(1.0, seed);
     15     }
     16     printf("gamma(1) summary statistics = mean: %f, std: %f\n", array_mean(gamma_array, n), array_std(gamma_array, n));
     17     printf("\n");
     18 
     19     double* beta_array = malloc(sizeof(double) * (size_t)n);
     20     for (int i = 0; i < n; i++) {
     21         beta_array[i] = sample_beta(1, 2.0, seed);
     22     }
     23     printf("beta(1,2) summary statistics: mean: %f, std: %f\n", array_mean(beta_array, n), array_std(beta_array, n));
     24     printf("\n");
     25 
     26     free(gamma_array);
     27     free(beta_array);
     28     free(seed);
     29 }