projects
/
physik
/
nlsop.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
some more bugfixes + testing
[physik/nlsop.git]
/
random.c
diff --git
a/random.c
b/random.c
index
6690531
..
468a3a5
100644
(file)
--- a/
random.c
+++ b/
random.c
@@
-55,6
+55,9
@@
u32 get_rand(u32 max)
{
if(c_ptr>=b_ptr+BUFSIZE)
{
{
if(c_ptr>=b_ptr+BUFSIZE)
{
+#ifdef MORE_PRINTF
+ printf("getting another %d bytes of random data ...\n",BUFSIZE);
+#endif
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 */
@@
-63,6
+66,9
@@
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;
+#ifdef MORE_PRINTF
+ printf("got it!\n");
+#endif
}
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))));
@@
-73,3
+79,18
@@
u32 get_rand_lgp(u32 max,double a,double b)
return((u32)(1.0*max*(-1.0*b+sqrt(b*b+2*a*((b+a/2)*get_rand(URAND_MAX)/((long long unsigned int)URAND_MAX+1))))/a));
}
return((u32)(1.0*max*(-1.0*b+sqrt(b*b+2*a*((b+a/2)*get_rand(URAND_MAX)/((long long unsigned int)URAND_MAX+1))))/a));
}
+u32 get_rand_reject(u32 max_x,u32 max_y,u32 *graph)
+{
+ u32 x,y;
+ unsigned char ok;
+
+ ok=0;
+ while(!ok)
+ {
+ x=get_rand(max_x);
+ y=get_rand(max_y);
+ if(y<=graph[x]) ok=1;
+ }
+
+ return x;
+}