testrandomu.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054 #ifdef HAVE_CONFIG_H
00055 #include <config.h>
00056 #endif
00057
00059 #define _GNU_SOURCE
00060
00061
00062 #ifdef HAVE_SYS_TYPES_H
00063 #include <sys/types.h>
00064 #endif
00065 #ifdef HAVE_SYS_STAT_H
00066 #include <sys/stat.h>
00067 #endif
00068 #ifdef HAVE_FCNTL_H
00069 #include <fcntl.h>
00070 #endif
00071 #ifdef HAVE_UNSTD_H
00072 #include <unistd.h>
00073 #endif
00074 #ifdef HAVE_STDIO_H
00075 #include <stdio.h>
00076 #endif
00077 #ifdef HAVE_STRING_H
00078 #include <string.h>
00079 #endif
00080 #ifdef HAVE_STDLIB_H
00081 #include <stdlib.h>
00082 #endif
00083 #ifdef HAVE_MATH_H
00084 #include <math.h>
00085 #endif
00086 #ifdef HAVE_TIME_H
00087 #include <time.h>
00088 #endif
00089 #ifdef HAVE_LIBGEN_H
00090 # include <libgen.h>
00091 #endif
00092
00093 #include <gsl/gsl_rng.h>
00094 #include <gsl/gsl_randist.h>
00095
00096 #include <utils.h>
00097
00099 void show_usage(char *pgm);
00100
00102 int main(int argc, char **argv)
00103 {
00111 unsigned long int random_num, i;
00112
00113 const gsl_rng_type *T;
00114 gsl_rng *rng;
00115
00116
00117 (void) banner(basename(argv[0]), "1.0", "BEGIN");
00118
00119
00120 for (i=1; i<argc; i++) {
00121 if ( !strcmp(argv[i], "-h") ) {
00122 (void) show_usage(basename(argv[0]));
00123 (void) banner(basename(argv[0]), "OK", "END");
00124 return 0;
00125 }
00126 else {
00127 (void) fprintf(stderr, "%s:: Wrong arg %s.\n\n", basename(argv[0]), argv[i]);
00128 (void) show_usage(basename(argv[0]));
00129 (void) banner(basename(argv[0]), "ABORT", "END");
00130 (void) abort();
00131 }
00132 }
00133
00134
00135 T = gsl_rng_default;
00136 rng = gsl_rng_alloc(T);
00137 (void) gsl_rng_set(rng, time(NULL));
00138
00139 for (i=0; i<1000000; i++) {
00140 random_num = gsl_rng_uniform_int(rng, 100);
00141 (void) fprintf(stdout,"%ld\n",random_num);
00142 }
00143
00144 (void) gsl_rng_free(rng);
00145
00146
00147 (void) banner(basename(argv[0]), "OK", "END");
00148
00149 return 0;
00150 }
00151
00152
00156 void show_usage(char *pgm) {
00161 (void) fprintf(stderr, "%s: usage:\n", pgm);
00162 (void) fprintf(stderr, "-h: help\n");
00163
00164 }