line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
1
|
|
|
|
|
|
|
#define PERL_NO_GET_CONTEXT |
2
|
|
|
|
|
|
|
#include "EXTERN.h" |
3
|
|
|
|
|
|
|
#include "perl.h" |
4
|
|
|
|
|
|
|
#include "XSUB.h" |
5
|
|
|
|
|
|
|
|
6
|
|
|
|
|
|
|
#include "ppport.h" |
7
|
|
|
|
|
|
|
|
8
|
|
|
|
|
|
|
#include "sandy_rng.h" |
9
|
|
|
|
|
|
|
#include "RNG.h" |
10
|
|
|
|
|
|
|
|
11
|
|
|
|
|
|
|
MODULE = App::Sandy::RNG PACKAGE = App::Sandy::RNG |
12
|
|
|
|
|
|
|
|
13
|
|
|
|
|
|
|
PROTOTYPES: DISABLED |
14
|
|
|
|
|
|
|
|
15
|
|
|
|
|
|
|
RNGContext * |
16
|
|
|
|
|
|
|
App::Sandy::RNG::new (unsigned long seed) |
17
|
|
|
|
|
|
|
PREINIT: |
18
|
|
|
|
|
|
|
RNGContext *ctx; |
19
|
|
|
|
|
|
|
CODE: |
20
|
114
|
|
|
|
|
|
Newx (ctx, 1, RNGContext); |
21
|
114
|
|
|
|
|
|
ctx->rng = gsl_rng_alloc (gsl_rng_ranlxd2); |
22
|
114
|
|
|
|
|
|
gsl_rng_set (ctx->rng, seed); |
23
|
114
|
|
|
|
|
|
RETVAL = ctx; |
24
|
|
|
|
|
|
|
OUTPUT: |
25
|
|
|
|
|
|
|
RETVAL |
26
|
|
|
|
|
|
|
|
27
|
|
|
|
|
|
|
void |
28
|
|
|
|
|
|
|
DESTROY (RNGContext *ctx) |
29
|
|
|
|
|
|
|
CODE: |
30
|
114
|
|
|
|
|
|
Safefree (ctx->rng); |
31
|
114
|
|
|
|
|
|
Safefree (ctx); |
32
|
|
|
|
|
|
|
|
33
|
|
|
|
|
|
|
void |
34
|
|
|
|
|
|
|
set (RNGContext *ctx, unsigned long seed) |
35
|
|
|
|
|
|
|
CODE: |
36
|
0
|
|
|
|
|
|
gsl_rng_set (ctx->rng, seed); |
37
|
|
|
|
|
|
|
|
38
|
|
|
|
|
|
|
unsigned long |
39
|
|
|
|
|
|
|
max (RNGContext *ctx) |
40
|
|
|
|
|
|
|
CODE: |
41
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_max (ctx->rng); |
42
|
|
|
|
|
|
|
OUTPUT: |
43
|
|
|
|
|
|
|
RETVAL |
44
|
|
|
|
|
|
|
|
45
|
|
|
|
|
|
|
unsigned long |
46
|
|
|
|
|
|
|
min (RNGContext *ctx) |
47
|
|
|
|
|
|
|
CODE: |
48
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_min (ctx->rng); |
49
|
|
|
|
|
|
|
OUTPUT: |
50
|
|
|
|
|
|
|
RETVAL |
51
|
|
|
|
|
|
|
|
52
|
|
|
|
|
|
|
size_t |
53
|
|
|
|
|
|
|
size (RNGContext *ctx) |
54
|
|
|
|
|
|
|
CODE: |
55
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_size (ctx->rng); |
56
|
|
|
|
|
|
|
OUTPUT: |
57
|
|
|
|
|
|
|
RETVAL |
58
|
|
|
|
|
|
|
|
59
|
|
|
|
|
|
|
char * |
60
|
|
|
|
|
|
|
name (RNGContext *ctx) |
61
|
|
|
|
|
|
|
CODE: |
62
|
0
|
|
|
|
|
|
RETVAL = (char *) gsl_rng_name (ctx->rng); |
63
|
|
|
|
|
|
|
OUTPUT: |
64
|
|
|
|
|
|
|
RETVAL |
65
|
|
|
|
|
|
|
|
66
|
|
|
|
|
|
|
unsigned long |
67
|
|
|
|
|
|
|
get (RNGContext *ctx) |
68
|
|
|
|
|
|
|
CODE: |
69
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_get (ctx->rng); |
70
|
|
|
|
|
|
|
OUTPUT: |
71
|
|
|
|
|
|
|
RETVAL |
72
|
|
|
|
|
|
|
|
73
|
|
|
|
|
|
|
double |
74
|
|
|
|
|
|
|
uniform (RNGContext *ctx) |
75
|
|
|
|
|
|
|
CODE: |
76
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_uniform (ctx->rng); |
77
|
|
|
|
|
|
|
OUTPUT: |
78
|
|
|
|
|
|
|
RETVAL |
79
|
|
|
|
|
|
|
|
80
|
|
|
|
|
|
|
double |
81
|
|
|
|
|
|
|
uniform_pos (RNGContext *ctx) |
82
|
|
|
|
|
|
|
CODE: |
83
|
0
|
|
|
|
|
|
RETVAL = gsl_rng_uniform_pos (ctx->rng); |
84
|
|
|
|
|
|
|
OUTPUT: |
85
|
|
|
|
|
|
|
RETVAL |
86
|
|
|
|
|
|
|
|
87
|
|
|
|
|
|
|
double |
88
|
|
|
|
|
|
|
ran_gaussian (RNGContext *ctx, double sigma) |
89
|
|
|
|
|
|
|
CODE: |
90
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_gaussian (ctx->rng, sigma); |
91
|
|
|
|
|
|
|
OUTPUT: |
92
|
|
|
|
|
|
|
RETVAL |
93
|
|
|
|
|
|
|
|
94
|
|
|
|
|
|
|
double |
95
|
|
|
|
|
|
|
ran_gaussian_ratio_method (RNGContext *ctx, double sigma) |
96
|
|
|
|
|
|
|
CODE: |
97
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_gaussian_ratio_method (ctx->rng, sigma); |
98
|
|
|
|
|
|
|
OUTPUT: |
99
|
|
|
|
|
|
|
RETVAL |
100
|
|
|
|
|
|
|
|
101
|
|
|
|
|
|
|
double |
102
|
|
|
|
|
|
|
ran_gaussian_pdf (double x, double sigma) |
103
|
|
|
|
|
|
|
CODE: |
104
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_gaussian_pdf (x, sigma); |
105
|
|
|
|
|
|
|
OUTPUT: |
106
|
|
|
|
|
|
|
RETVAL |
107
|
|
|
|
|
|
|
|
108
|
|
|
|
|
|
|
double |
109
|
|
|
|
|
|
|
ran_ugaussian (RNGContext *ctx) |
110
|
|
|
|
|
|
|
CODE: |
111
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_ugaussian (ctx->rng); |
112
|
|
|
|
|
|
|
OUTPUT: |
113
|
|
|
|
|
|
|
RETVAL |
114
|
|
|
|
|
|
|
|
115
|
|
|
|
|
|
|
double |
116
|
|
|
|
|
|
|
ran_ugaussian_ratio_method (RNGContext *ctx) |
117
|
|
|
|
|
|
|
CODE: |
118
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_ugaussian_ratio_method (ctx->rng); |
119
|
|
|
|
|
|
|
OUTPUT: |
120
|
|
|
|
|
|
|
RETVAL |
121
|
|
|
|
|
|
|
|
122
|
|
|
|
|
|
|
double |
123
|
|
|
|
|
|
|
ran_ugaussian_pdf (double x) |
124
|
|
|
|
|
|
|
CODE: |
125
|
0
|
|
|
|
|
|
RETVAL = gsl_ran_ugaussian_pdf (x); |
126
|
|
|
|
|
|
|
OUTPUT: |
127
|
|
|
|
|
|
|
RETVAL |
128
|
|
|
|
|
|
|
|
129
|
|
|
|
|
|
|
unsigned long |
130
|
|
|
|
|
|
|
get_n (RNGContext *ctx, unsigned long n) |
131
|
|
|
|
|
|
|
INIT: |
132
|
123840
|
50
|
|
|
|
|
if (n <= 0) |
133
|
0
|
|
|
|
|
|
croak ("n must be greater than zero"); |
134
|
|
|
|
|
|
|
CODE: |
135
|
123840
|
|
|
|
|
|
RETVAL = gsl_rng_get (ctx->rng) % n; |
136
|
|
|
|
|
|
|
OUTPUT: |
137
|
|
|
|
|
|
|
RETVAL |
138
|
|
|
|
|
|
|
|
139
|
|
|
|
|
|
|
double |
140
|
|
|
|
|
|
|
get_norm (RNGContext *ctx, double mean, double stdd) |
141
|
|
|
|
|
|
|
CODE: |
142
|
3332
|
|
|
|
|
|
RETVAL = mean + (int) (gsl_ran_gaussian (ctx->rng, stdd) + 0.5); |
143
|
|
|
|
|
|
|
OUTPUT: |
144
|
|
|
|
|
|
|
RETVAL |