diff options
Diffstat (limited to 'SD-VBS/common/c/randnWrapper.c')
-rw-r--r-- | SD-VBS/common/c/randnWrapper.c | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/SD-VBS/common/c/randnWrapper.c b/SD-VBS/common/c/randnWrapper.c new file mode 100644 index 0000000..4701b0e --- /dev/null +++ b/SD-VBS/common/c/randnWrapper.c | |||
@@ -0,0 +1,40 @@ | |||
1 | /******************************** | ||
2 | Author: Sravanthi Kota Venkata | ||
3 | ********************************/ | ||
4 | |||
5 | #include "sdvbs_common.h" | ||
6 | #include <math.h> | ||
7 | F2D* randnWrapper(int m, int n) | ||
8 | { | ||
9 | F2D *out; | ||
10 | float seed; | ||
11 | int i,j; | ||
12 | |||
13 | out = fSetArray(m, n, 0); | ||
14 | seed = 0.9; | ||
15 | |||
16 | for(i=0; i<m; i++) | ||
17 | { | ||
18 | for(j=0; j<n; j++) | ||
19 | { | ||
20 | if(i<j) | ||
21 | subsref(out,i,j) = seed * ((i+1.0)/(j+1.0)); | ||
22 | else | ||
23 | subsref(out,i,j) = seed * ((j+1.0)/(i+1.0)); | ||
24 | } | ||
25 | } | ||
26 | |||
27 | for(i=0; i<m ;i++) | ||
28 | { | ||
29 | for(j=0; j<n; j++) | ||
30 | { | ||
31 | float w; | ||
32 | w = subsref(out,i,j); | ||
33 | w = ((-2.0 * log(w))/w); | ||
34 | subsref(out,i,j) = w; | ||
35 | } | ||
36 | } | ||
37 | |||
38 | return out; | ||
39 | } | ||
40 | |||