1 /* general functions */
18 double absolute_value(double l) {
25 #define __M ((2<<31)-1)
28 /* global variable -- I_{k+1} dedpends on I_{k} */
29 static int idum=123456;
31 int my_rand(int max) {
34 idum=__A*(idum-h*__Q)-h*__R;
36 return(idum*(max/__M));
39 int my_srand(int seed) {
40 printf("debug: seed = %d\n",seed);
44 double my_nrand(void) {
47 idum=__A*(idum-h*__Q)-h*__R;