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 | |||
