squiggle.c/scratchpad/scratchpad.c

35 lines
830 B
C
Raw Normal View History

2023-09-27 13:10:40 +00:00
#include "../squiggle.h"
2023-11-29 23:08:36 +00:00
#include "../squiggle_more.h"
2023-09-27 13:10:40 +00:00
#include <math.h>
#include <stdint.h>
#include <stdio.h>
#include <stdlib.h>
int main()
{
// set randomness seed
uint64_t* seed = malloc(sizeof(uint64_t));
*seed = 1000; // xorshift can't start with a seed of 0
/*
for (int i = 0; i < 100; i++) {
double draw = sample_unit_uniform(seed);
printf("%f\n", draw);
}*/
// Test division
// printf("\n%d\n", 10 % 3);
//
2023-11-29 23:08:36 +00:00
int n_samples = 100, n_threads = 16;
double* results = malloc(n_samples * sizeof(double));
double sampler_scratchpad(uint64_t* seed){
return 1;
}
parallel_sampler(sampler_scratchpad, results, n_threads, n_samples);
for(int i=0; i<n_samples; i++){
printf("Sample %d: %f\n", i, results[i]);
2023-11-29 23:08:36 +00:00
}
2023-09-27 13:10:40 +00:00
free(seed);
}