projects
/
physik
/
nlsop.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
support: meanvalue/sum of multiple fts
[physik/nlsop.git]
/
random.c
diff --git
a/random.c
b/random.c
index
cb3e900
..
a664fa7
100644
(file)
--- a/
random.c
+++ b/
random.c
@@
-12,6
+12,10
@@
#include <stdlib.h>
#include <math.h>
#include <stdlib.h>
#include <math.h>
+#ifdef USE_LIBC_RAND
+#include <time.h>
+#endif
+
#include "random.h"
static int rand_fd;
#include "random.h"
static int rand_fd;
@@
-21,11
+25,13
@@
int rand_init(char *rf)
{
if(rf==NULL)
{
{
if(rf==NULL)
{
+#ifndef USE_LIBC_RAND
if((rand_fd=open("/dev/urandom",O_RDONLY))<0)
{
puts("cannot open /dev/urandom");
return -1;
}
if((rand_fd=open("/dev/urandom",O_RDONLY))<0)
{
puts("cannot open /dev/urandom");
return -1;
}
+#endif
} else
{
if((rand_fd=open(rf,O_RDONLY))<0)
} else
{
if((rand_fd=open(rf,O_RDONLY))<0)
@@
-46,7
+52,9
@@
int rand_init(char *rf)
int rand_close(void)
{
int rand_close(void)
{
+#ifndef USE_LIBC_RAND
close(rand_fd);
close(rand_fd);
+#endif
return 1;
}
return 1;
}
@@
-56,8
+64,14
@@
u32 get_rand(u32 max)
if(c_ptr>=b_ptr+BUFSIZE)
{
#ifdef MORE_PRINTF
if(c_ptr>=b_ptr+BUFSIZE)
{
#ifdef MORE_PRINTF
- printf("getting another %d bytes of random data ...",BUFSIZE);
+ printf("getting another %d bytes of random data ...
\n
",BUFSIZE);
#endif
#endif
+#ifdef USE_LIBC_RAND
+ c_ptr=b_ptr;
+ srand((int)time(NULL));
+ while(c_ptr<b_ptr+BUFSIZE) *(c_ptr++)=(u32)rand();
+ c_ptr=b_ptr;
+#else
if(read(rand_fd,b_ptr,BUFSIZE*sizeof(u32))<BUFSIZE*sizeof(u32))
{
/* -> assume random file, end reached */
if(read(rand_fd,b_ptr,BUFSIZE*sizeof(u32))<BUFSIZE*sizeof(u32))
{
/* -> assume random file, end reached */
@@
-66,12
+80,17
@@
u32 get_rand(u32 max)
read(rand_fd,b_ptr,BUFSIZE*sizeof(u32));
}
c_ptr=b_ptr;
read(rand_fd,b_ptr,BUFSIZE*sizeof(u32));
}
c_ptr=b_ptr;
+#endif
#ifdef MORE_PRINTF
#ifdef MORE_PRINTF
- printf("
done
\n");
+ printf("
got it!
\n");
#endif
}
#endif
}
+#ifdef USE_LIBC_RAND
+ return((u32)(*(c_ptr++)*(max*1.0/((long long unsigned int)RAND_MAX+1))));
+#else
return((u32)(*(c_ptr++)*(max*1.0/((long long unsigned int)URAND_MAX+1))));
return((u32)(*(c_ptr++)*(max*1.0/((long long unsigned int)URAND_MAX+1))));
+#endif
}
u32 get_rand_lgp(u32 max,double a,double b)
}
u32 get_rand_lgp(u32 max,double a,double b)