diff options
author | imushir | 2016-01-18 11:47:44 +0530 |
---|---|---|
committer | imushir | 2016-01-18 11:47:44 +0530 |
commit | bcb2bf1c014a1078f4306cd0111f4e596fe934ce (patch) | |
tree | 2aaae4926f8354f89b2061599d62b5d650e1535e /2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c | |
parent | 03a27e6884942182e984433fd2ebba7e5f55a4d9 (diff) | |
download | Scilab2C-bcb2bf1c014a1078f4306cd0111f4e596fe934ce.tar.gz Scilab2C-bcb2bf1c014a1078f4306cd0111f4e596fe934ce.tar.bz2 Scilab2C-bcb2bf1c014a1078f4306cd0111f4e596fe934ce.zip |
added support for uin8,uint16,int8,int16 and not in GeneratedSetup.sci,u8cmd_digital_outs.c,getAllInterfaces.sci,getAllHeaders.sci
Diffstat (limited to '2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c')
-rw-r--r-- | 2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c b/2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c new file mode 100644 index 00000000..55df9093 --- /dev/null +++ b/2.3-1/src/c/auxiliaryFunctions/rand/i8rands.c @@ -0,0 +1,42 @@ +/* Scilab2C FOSSEE IIT Bombay */ +#include "rand.h" +#include <stdio.h> + +double i8rands(void) +{ + +int m=1; +const int itwo=2; +static int m2=0,halfm=0,ia=0,ic=0,mic=0,iy=0; +static double s=0.0; +if(m2==0) +{ + /*if first entry,compute machine integer word length*/ + while(m>m2) + { + m2=m; + m=itwo*m2; + } + halfm = m2; + /* compute multiplier and increment for linear congruential methos */ + ia = 8*(int)(halfm*atan(1.0)/8.0) + 5; + ic = 2*(int)(halfm*(0.5-sqrt(3.0)/6.0)) +1; + mic = (m2 - ic) + m2; + /* s is the scale factor for converting to floating point*/ + s = 0.5/halfm; + +} + +/* compute next random number */ + iy = iy*ia; +/* the following statement is for computers which do not allow interger overflow on addition*/ +if(iy > mic) iy = (iy - m2) - m2; + iy = iy + ic; + + /* the following statement is for computers where the word length for addition is greater than for multiplication */ + if(iy/2 > m2) iy = (iy - m2) - m2; + /* the following statement is for computers where integer overflow affects the sign bit */ + if(iy < 0) iy = (iy + m2) + m2; + return (double)iy*s; + +} |