1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
/* ------------------------------------------------------------- */
/* File: example1.c */
/* ------------------------------------------------------------- */
/* Include UNURAN header files. */
#include <unuran.h>
#include <unuran_urng_rngstreams.h>
/* ------------------------------------------------------------- */
int main(void)
{
int i; /* loop variable */
double x; int k; /* will hold the random number */
/* Declare UNU.RAN objects. */
UNUR_GEN *gen1, *gen2, *gen3; /* generator objects */
/* Declare objects for uniform random number generators. */
UNUR_URNG *urng2; /* uniform RN generator object */
/* -- Optional: Set seed for RNGSTREAMS library -------------- */
/* The RNGSTREAMS library sets a package seed. */
unsigned long seed[] = {111u, 222u, 333u, 444u, 555u, 666u};
RngStream_SetPackageSeed(seed);
/* -- Example 1 ---------------------------------------------- */
/* Beta distribution with shape parameters 2 and 3. */
/* Use method 'AUTO' (AUTOmatic). */
/* Create generator object. */
gen1 = unur_str2gen("beta(2,3)");
/* It is important to check if the creation of the generator */
/* object was successful. Otherwise `gen1' is the NULL pointer */
/* and would cause a segmentation fault if used for sampling. */
if (gen1 == NULL) {
fprintf(stderr, "ERROR: cannot create generator object\n");
exit (EXIT_FAILURE);
}
/* Now you can use the generator object `gen1' to sample from */
/* the target distribution. Eg.: */
for (i=0; i<10; i++) {
x = unur_sample_cont(gen1);
printf("%f\n",x);
}
/* -- Example 2 ---------------------------------------------- */
/* Student's t distribution with 3 degrees of freedom. */
/* Use method 'TDR' (Transformed Density Rejection) with */
/* "immediate acception" */
gen2 = unur_str2gen("student(3) & method=TDR; variant_ia");
if (gen2 == NULL) exit (EXIT_FAILURE);
/* However, this time we use a (new) independent stream of */
/* uniformrandom numbers. */
urng2 = unur_urng_rngstream_new("urng2");
unur_chg_urng( gen2, urng2 );
/* Draw a sample. */
for (i=0; i<10; i++) {
x = unur_sample_cont(gen2); printf("%f\n",x);
}
/* -- Example 3 ---------------------------------------------- */
/* Discrete distribution with given probability vector. */
/* Use method 'DGT' (Discrete Guide Table method). */
gen3 = unur_str2gen("discr; pv=(0.5,1.5,1.0,0.3) & method=DGT");
if (gen3 == NULL) exit (EXIT_FAILURE);
/* we use the default URNG again. So there is nothing to do. */
/* Draw a sample. Notice that we get integers! */
for (i=0; i<10; i++) {
k = unur_sample_discr(gen3); printf("%d\n",k);
}
/* -- Call destructor ---------------------------------------- */
/* When generators are not needed any they can be destroyed. */
unur_free(gen1);
unur_free(gen2); unur_urng_free(urng2);
unur_free(gen3);
exit (EXIT_SUCCESS);
} /* end of main() */
/* ------------------------------------------------------------- */
|
