Giu/analistica
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

misc: deleted directory

Browse Source
This commit is contained in:
Giù Marcer 2020-04-27 23:23:08 +02:00 committed by rnhmjoj
parent c560f3c57c
commit 1b39f1cf0b
2 changed files with 0 additions and 73 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
misc/pdfs-plot.py
View File

@ -1,15 +0,0 @@
#!/usr/bin/env python
import matplotlib.pyplot as plt
import numpy as np
import sys
lows, upps, values = np.loadtxt(sys.stdin, unpack=True,)
dim = upps[1] - lows[1]
values = values/np.sum((values)*dim)
plt.bar(lows, align='edge', width=dim, height=values,
color='#dbbf0d', edgecolor='#595856')
plt.annotate("N = 100000\nbins=100", (15, 0.16))
plt.ylabel("pdf(x)")
plt.xlabel("x")
plt.show()

58
misc/pdfs.c
View File

@ -1,58 +0,0 @@
// This program generates a collection of random numbers which follow the
// distribution whose name is given in input by the user when calling the
// program and outputs an histogram of the result. Maximum and minimum of
// the range must also be given. The random number distribution functions
// are taken from the GSL library at:
// https://www.gnu.org/software/gsl/doc/html/randist.html
#include <stdio.h>
#include <string.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_histogram.h>
int main (int argc, char** argv)
{
// An histogram with 100 bins which range from 0 to 100 is set:
size_t bins = 100;
int low = atoi(argv[2]);
int upp = atoi(argv[3]);
gsl_histogram * histo = gsl_histogram_alloc(bins);
gsl_histogram_set_ranges_uniform(histo, low, upp);
// The random number generator is set:
gsl_rng* gen = gsl_rng_alloc(gsl_rng_taus);
// 100'000 random numbers are generated and put into the histogram:
size_t N = 100000;
for (size_t i = 0; i < N; i++)
{
if (strcmp(argv[1], "binomial") == 0)
gsl_histogram_increment(histo, gsl_ran_binomial(gen, 0.3, 10));
if (strcmp(argv[1], "poisson") == 0)
gsl_histogram_increment(histo, gsl_ran_poisson(gen, 1));
if (strcmp(argv[1], "uniform") == 0)
gsl_histogram_increment(histo, gsl_rng_get(gen)*100);
if (strcmp(argv[1], "gaussian") == 0)
gsl_histogram_increment(histo, gsl_ran_gaussian(gen, 5)+30);
if (strcmp(argv[1], "wigner") == 0)
gsl_histogram_increment(histo, gsl_ran_cauchy(gen, 10)+20);
if (strcmp(argv[1], "landau") == 0)
gsl_histogram_increment(histo, gsl_ran_landau(gen));
if (strcmp(argv[1], "chi2") == 0)
gsl_histogram_increment(histo, gsl_ran_gamma(gen, 5, 2));
if (strcmp(argv[1], "exponential") == 0)
gsl_histogram_increment(histo, gsl_ran_exponential(gen, 5));
if (strcmp(argv[1], "student") == 0)
gsl_histogram_increment(histo, gsl_ran_tdist(gen, 3));
if (strcmp(argv[1], "fischer") == 0)
gsl_histogram_increment(histo, gsl_ran_fdist(gen, 3, 4));
}
// The histogram is given in output:
gsl_histogram_fprintf (stdout, histo, "%g", "%g");
// The memory is freed.
gsl_rng_free(gen);
gsl_histogram_free(histo);
return 0;
}