| 1 |
chuckv |
4 |
|
| 2 |
|
|
|
| 3 |
|
|
#ifndef _sprngf_h_ |
| 4 |
|
|
|
| 5 |
|
|
#define SPRNG_LFG 0 |
| 6 |
|
|
#define SPRNG_LCG 1 |
| 7 |
|
|
#define SPRNG_LCG64 2 |
| 8 |
|
|
#define SPRNG_CMRG 3 |
| 9 |
|
|
#define SPRNG_MLFG 4 |
| 10 |
|
|
#define SPRNG_PMLCG 5 |
| 11 |
|
|
#define DEFAULT_RNG_TYPE SPRNG_LFG |
| 12 |
|
|
|
| 13 |
|
|
#define SPRNG_DEFAULT 0 |
| 14 |
|
|
#define CRAYLCG 0 |
| 15 |
|
|
#define DRAND48 1 |
| 16 |
|
|
#define FISH1 2 |
| 17 |
|
|
#define FISH2 3 |
| 18 |
|
|
#define FISH3 4 |
| 19 |
|
|
#define FISH4 5 |
| 20 |
|
|
#define FISH5 6 |
| 21 |
|
|
#define LECU1 0 |
| 22 |
|
|
#define LECU2 1 |
| 23 |
|
|
#define LECU3 2 |
| 24 |
|
|
#define LAG1279 0 |
| 25 |
|
|
#define LAG17 1 |
| 26 |
|
|
#define LAG31 2 |
| 27 |
|
|
#define LAG55 3 |
| 28 |
|
|
#define LAG63 4 |
| 29 |
|
|
#define LAG127 5 |
| 30 |
|
|
#define LAG521 6 |
| 31 |
|
|
#define LAG521B 7 |
| 32 |
|
|
#define LAG607 8 |
| 33 |
|
|
#define LAG607B 9 |
| 34 |
|
|
#define LAG1279B 10 |
| 35 |
|
|
|
| 36 |
|
|
#ifdef CHECK_POINTERS |
| 37 |
|
|
#define CHECK 1 |
| 38 |
|
|
#else |
| 39 |
|
|
#define CHECK 0 |
| 40 |
|
|
#endif /* ifdef CHECK_POINTERS */ |
| 41 |
|
|
|
| 42 |
|
|
#define MAX_PACKED_LENGTH 24000 |
| 43 |
|
|
|
| 44 |
|
|
#ifdef POINTER_SIZE |
| 45 |
|
|
#if POINTER_SIZE == 8 |
| 46 |
|
|
#define SPRNG_POINTER integer*8 |
| 47 |
|
|
#else |
| 48 |
|
|
#define SPRNG_POINTER integer*4 |
| 49 |
|
|
#endif |
| 50 |
|
|
#else |
| 51 |
|
|
#define SPRNG_POINTER integer*4 |
| 52 |
|
|
#endif /* ifdef POINTER_SIZE */ |
| 53 |
|
|
|
| 54 |
|
|
#ifdef USE_MPI |
| 55 |
|
|
#define make_sprng_seed fseed_mpi |
| 56 |
|
|
#else |
| 57 |
|
|
#define make_sprng_seed fmake_new_seed |
| 58 |
|
|
#endif |
| 59 |
|
|
|
| 60 |
|
|
#endif /* ifdef _sprng_h */ |
| 61 |
|
|
|
| 62 |
|
|
#ifdef USE_MPI |
| 63 |
|
|
external fseed_mpi |
| 64 |
|
|
integer fseed_mpi |
| 65 |
|
|
#else |
| 66 |
|
|
external fmake_new_seed |
| 67 |
|
|
integer fmake_new_seed |
| 68 |
|
|
#endif |
| 69 |
|
|
|
| 70 |
|
|
#ifndef DEFAULTINT |
| 71 |
|
|
#define DEFAULTINT |
| 72 |
|
|
#endif |
| 73 |
|
|
#ifndef FLOAT_GEN |
| 74 |
|
|
#define DBLGEN |
| 75 |
|
|
#endif |
| 76 |
|
|
|
| 77 |
|
|
#if defined(SIMPLE_SPRNG) |
| 78 |
|
|
#undef DEFAULTINT |
| 79 |
|
|
|
| 80 |
|
|
#ifndef _sprngf_h_ |
| 81 |
|
|
#define pack_sprng fpack_rng_simple |
| 82 |
|
|
#define unpack_sprng funpack_rng_simple |
| 83 |
|
|
#ifdef USE_MPI |
| 84 |
|
|
#define isprng fget_rn_int_simmpi |
| 85 |
|
|
#define init_sprng finit_rng_simmpi |
| 86 |
|
|
#else |
| 87 |
|
|
#define isprng fget_rn_int_sim |
| 88 |
|
|
#define init_sprng finit_rng_sim |
| 89 |
|
|
#endif /* ifdef USE_MPI */ |
| 90 |
|
|
#define print_sprng fprint_rng_simple |
| 91 |
|
|
|
| 92 |
|
|
#if defined(FLOAT_GEN) && defined(USE_MPI) |
| 93 |
|
|
#define sprng fget_rn_flt_simmpi |
| 94 |
|
|
#endif |
| 95 |
|
|
#if defined(FLOAT_GEN) && !defined(USE_MPI) |
| 96 |
|
|
#define sprng fget_rn_flt_sim |
| 97 |
|
|
#endif |
| 98 |
|
|
#if defined(DBLGEN) && defined(USE_MPI) |
| 99 |
|
|
#define sprng fget_rn_dbl_simmpi |
| 100 |
|
|
#endif |
| 101 |
|
|
#if defined(DBLGEN) && !defined(USE_MPI) |
| 102 |
|
|
#define sprng fget_rn_dbl_sim |
| 103 |
|
|
#endif |
| 104 |
|
|
|
| 105 |
|
|
#endif /* ifdef _sprng_h */ |
| 106 |
|
|
external isprng |
| 107 |
|
|
external fget_rn_dbl_sim, fget_rn_flt_sim |
| 108 |
|
|
external init_sprng, fpack_rng_simple |
| 109 |
|
|
external funpack_rng_simple, fprint_rng_simple |
| 110 |
|
|
#ifdef USE_MPI |
| 111 |
|
|
external fget_rn_flt_simmpi, fget_rn_dbl_simmpi |
| 112 |
|
|
real*4 fget_rn_flt_simmpi |
| 113 |
|
|
real*8 fget_rn_dbl_simmpi |
| 114 |
|
|
#endif |
| 115 |
|
|
integer isprng,fpack_rng_simple,fprint_rng_simple |
| 116 |
|
|
SPRNG_POINTER init_sprng, funpack_rng_simple |
| 117 |
|
|
real*4 fget_rn_flt_sim |
| 118 |
|
|
real*8 fget_rn_dbl_sim |
| 119 |
|
|
#endif |
| 120 |
|
|
|
| 121 |
|
|
#if defined(CHECK_POINTERS) |
| 122 |
|
|
#undef DEFAULTINT |
| 123 |
|
|
external fget_rn_int_ptr, fget_rn_flt_ptr, fget_rn_dbl_ptr |
| 124 |
|
|
external fspawn_rng_ptr, ffree_rng_ptr, finit_rng_ptr |
| 125 |
|
|
external fpack_rng_ptr, funpack_rng_ptr, fprint_rng_ptr |
| 126 |
|
|
|
| 127 |
|
|
integer fget_rn_int_ptr, ffree_rng_ptr, fpack_rng_ptr |
| 128 |
|
|
SPRNG_POINTER finit_rng_ptr, funpack_rng_ptr |
| 129 |
|
|
integer fspawn_rng_ptr, fprint_rng_ptr |
| 130 |
|
|
real*4 fget_rn_flt_ptr |
| 131 |
|
|
real*8 fget_rn_dbl_ptr |
| 132 |
|
|
|
| 133 |
|
|
#ifndef _sprngf_h_ |
| 134 |
|
|
#define isprng fget_rn_int_ptr |
| 135 |
|
|
#define free_sprng ffree_rng_ptr |
| 136 |
|
|
#define spawn_sprng(A,B,C) fspawn_rng_ptr(A,B,C,CHECK) |
| 137 |
|
|
#define pack_sprng fpack_rng_ptr |
| 138 |
|
|
#define unpack_sprng funpack_rng_ptr |
| 139 |
|
|
#define init_sprng finit_rng_ptr |
| 140 |
|
|
#define print_sprng fprint_rng_ptr |
| 141 |
|
|
#ifdef FLOAT_GEN |
| 142 |
|
|
#define sprng fget_rn_flt_ptr |
| 143 |
|
|
#else |
| 144 |
|
|
#define sprng fget_rn_dbl_ptr |
| 145 |
|
|
#endif |
| 146 |
|
|
#endif |
| 147 |
|
|
#endif |
| 148 |
|
|
|
| 149 |
|
|
#if defined(DEFAULTINT) |
| 150 |
|
|
external fget_rn_int, fget_rn_flt, fget_rn_dbl |
| 151 |
|
|
external fspawn_rng, ffree_rng, finit_rng |
| 152 |
|
|
external fpack_rng, funpack_rng, fprint_rng |
| 153 |
|
|
|
| 154 |
|
|
integer fget_rn_int, ffree_rng, fpack_rng |
| 155 |
|
|
SPRNG_POINTER finit_rng, funpack_rng |
| 156 |
|
|
integer fspawn_rng, fprint_rng |
| 157 |
|
|
real*4 fget_rn_flt |
| 158 |
|
|
real*8 fget_rn_dbl |
| 159 |
|
|
|
| 160 |
|
|
#ifndef _sprngf_h_ |
| 161 |
|
|
#define isprng fget_rn_int |
| 162 |
|
|
#define free_sprng ffree_rng |
| 163 |
|
|
#define spawn_sprng(A,B,C) fspawn_rng(A,B,C,CHECK) |
| 164 |
|
|
#define pack_sprng fpack_rng |
| 165 |
|
|
#define unpack_sprng funpack_rng |
| 166 |
|
|
#define init_sprng finit_rng |
| 167 |
|
|
#define print_sprng fprint_rng |
| 168 |
|
|
#ifdef FLOAT_GEN |
| 169 |
|
|
#define sprng fget_rn_flt |
| 170 |
|
|
#else |
| 171 |
|
|
#define sprng fget_rn_dbl |
| 172 |
|
|
#endif |
| 173 |
|
|
#endif |
| 174 |
|
|
|
| 175 |
|
|
#endif |
| 176 |
|
|
|
| 177 |
|
|
|
| 178 |
|
|
|
| 179 |
|
|
#ifndef _sprngf_h_ |
| 180 |
|
|
#define _sprngf_h_ |
| 181 |
|
|
#endif |
