src/basesuif/F77/stanford.c

   1 /* DEM */
   2 /* functions sf2c creates that we don't want macro-expanded */
   3
   4 #include <stdarg.h>
   5 #include "sf2c.h"
   6
   7 /* MIN functions */
   8 integer i_min(long num_args, ...)
   9 {
  10         va_list ap;
  11         integer best,temp,i;
  12
  13         va_start(ap, num_args);
  14
  15         best = va_arg(ap,integer);
  16         for (i=2; i<=num_args; i++) {
  17                 temp = va_arg(ap,integer);
  18                 if (temp < best) best = temp;
  19         }
  20
  21         va_end(ap);
  22         return best;
  23 }
  24
  25 shortint h_min(long num_args, ...)
  26 {
  27         va_list ap;
  28         shortint best,temp;
  29         integer i;
  30
  31         va_start(ap, num_args);
  32
  33         best = va_arg(ap,shortint);
  34         for (i=2; i<=num_args; i++) {
  35                 temp = va_arg(ap,shortint);
  36                 if (temp < best) best = temp;
  37         }
  38
  39         va_end(ap);
  40         return best;
  41 }
  42
  43 real r_min(long num_args, ...)
  44 {
  45         va_list ap;
  46         real best,temp;
  47         integer i;
  48
  49         va_start(ap, num_args);
  50
  51         best = va_arg(ap,real);
  52         for (i=2; i<=num_args; i++) {
  53                 temp = va_arg(ap,real);
  54                 if (temp < best) best = temp;
  55         }
  56
  57         va_end(ap);
  58         return best;
  59 }
  60
  61 doublereal d_min(long num_args, ...)
  62 {
  63         va_list ap;
  64         doublereal best,temp;
  65         integer i;
  66
  67         va_start(ap, num_args);
  68
  69         best = va_arg(ap,doublereal);
  70         for (i=2; i<=num_args; i++) {
  71                 temp = va_arg(ap,doublereal);
  72                 if (temp < best) best = temp;
  73         }
  74
  75         va_end(ap);
  76         return best;
  77 }
  78
  79 real ri_min(long num_args, ...)
  80 {
  81         va_list ap;
  82         real best;
  83         integer temp;
  84         integer i;
  85
  86         va_start(ap, num_args);
  87
  88         best = va_arg(ap,integer);
  89         for (i=2; i<=num_args; i++) {
  90                 temp = va_arg(ap,integer);
  91                 if (temp < best) best = temp;
  92         }
  93
  94         va_end(ap);
  95         return best;
  96 }
  97
  98 doublereal di_min(long num_args, ...)
  99 {
 100         va_list ap;
 101         doublereal best;
 102         integer temp;
 103         integer i;
 104
 105         va_start(ap, num_args);
 106
 107         best = va_arg(ap,integer);
 108         for (i=2; i<=num_args; i++) {
 109                 temp = va_arg(ap,integer);
 110                 if (temp < best) best = temp;
 111         }
 112
 113         va_end(ap);
 114         return best;
 115 }
 116
 117 integer ir_min(long num_args, ...)
 118 {
 119         va_list ap;
 120         real best,temp;
 121         integer i;
 122
 123         va_start(ap, num_args);
 124
 125         best = va_arg(ap,real);
 126         for (i=2; i<=num_args; i++) {
 127                 temp = va_arg(ap,real);
 128                 if (temp < best) best = temp;
 129         }
 130
 131         va_end(ap);
 132         return best;
 133 }
 134
 135 integer id_min(long num_args, ...)
 136 {
 137         va_list ap;
 138         doublereal best,temp;
 139         integer i;
 140
 141         va_start(ap, num_args);
 142
 143         best = va_arg(ap,doublereal);
 144         for (i=2; i<=num_args; i++) {
 145                 temp = va_arg(ap,doublereal);
 146                 if (temp < best) best = temp;
 147         }
 148
 149         va_end(ap);
 150         return best;
 151 }
 152
 153 /* MAX functions */
 154 integer i_max(long num_args, ...)
 155 {
 156         va_list ap;
 157         integer best,temp,i;
 158
 159         va_start(ap, num_args);
 160
 161         best = va_arg(ap,integer);
 162         for (i=2; i<=num_args; i++) {
 163                 temp = va_arg(ap,integer);
 164                 if (temp > best) best = temp;
 165         }
 166
 167         va_end(ap);
 168         return best;
 169 }
 170
 171 shortint h_max(long num_args, ...)
 172 {
 173         va_list ap;
 174         shortint best,temp;
 175         integer i;
 176
 177         va_start(ap, num_args);
 178
 179         best = va_arg(ap,shortint);
 180         for (i=2; i<=num_args; i++) {
 181                 temp = va_arg(ap,shortint);
 182                 if (temp > best) best = temp;
 183         }
 184
 185         va_end(ap);
 186         return best;
 187 }
 188
 189 real r_max(long num_args, ...)
 190 {
 191         va_list ap;
 192         real best,temp;
 193         integer i;
 194
 195         va_start(ap, num_args);
 196
 197         best = va_arg(ap,real);
 198         for (i=2; i<=num_args; i++) {
 199                 temp = va_arg(ap,real);
 200                 if (temp > best) best = temp;
 201         }
 202
 203         va_end(ap);
 204         return best;
 205 }
 206
 207 doublereal d_max(long num_args, ...)
 208 {
 209         va_list ap;
 210         doublereal best,temp;
 211         integer i;
 212
 213         va_start(ap, num_args);
 214
 215         best = va_arg(ap,doublereal);
 216         for (i=2; i<=num_args; i++) {
 217                 temp = va_arg(ap,doublereal);
 218                 if (temp > best) best = temp;
 219         }
 220
 221         va_end(ap);
 222         return best;
 223 }
 224
 225 real ri_max(long num_args, ...)
 226 {
 227         va_list ap;
 228         real best;
 229         integer temp;
 230         integer i;
 231
 232         va_start(ap, num_args);
 233
 234         best = va_arg(ap,integer);
 235         for (i=2; i<=num_args; i++) {
 236                 temp = va_arg(ap,integer);
 237                 if (temp > best) best = temp;
 238         }
 239
 240         va_end(ap);
 241         return best;
 242 }
 243
 244 doublereal di_max(long num_args, ...)
 245 {
 246         va_list ap;
 247         doublereal best;
 248         integer temp;
 249         integer i;
 250
 251         va_start(ap, num_args);
 252
 253         best = va_arg(ap,integer);
 254         for (i=2; i<=num_args; i++) {
 255                 temp = va_arg(ap,integer);
 256                 if (temp > best) best = temp;
 257         }
 258
 259         va_end(ap);
 260         return best;
 261 }
 262
 263 integer ir_max(long num_args, ...)
 264 {
 265         va_list ap;
 266         real best,temp;
 267         integer i;
 268
 269         va_start(ap, num_args);
 270
 271         best = va_arg(ap,real);
 272         for (i=2; i<=num_args; i++) {
 273                 temp = va_arg(ap,real);
 274                 if (temp > best) best = temp;
 275         }
 276
 277         va_end(ap);
 278         return best;
 279 }
 280
 281 integer id_max(long num_args, ...)
 282 {
 283         va_list ap;
 284         doublereal best,temp;
 285         integer i;
 286
 287         va_start(ap, num_args);
 288
 289         best = va_arg(ap,doublereal);
 290         for (i=2; i<=num_args; i++) {
 291                 temp = va_arg(ap,doublereal);
 292                 if (temp > best) best = temp;
 293         }
 294
 295         va_end(ap);
 296         return best;
 297 }
 298
 299
 300 /* ABS functions */
 301
 302 double sd_abs(x)
 303 doublereal x;
 304 {
 305 if(x >= 0)
 306         return(x);
 307 return(- x);
 308 }
 309
 310 double sc_abs(z)
 311 complex *z;
 312 {
 313 return sqrt((z->_r * z->_r) + (z->_i * z->_i));
 314 }
 315
 316 double sr_abs(x)
 317 real x;
 318 {
 319 if(x >= 0)
 320         return(x);
 321 return(- x);
 322 }
 323
 324 double sz_abs(z)
 325 doublecomplex *z;
 326 {
 327 return sqrt((z->_r * z->_r) + (z->_i * z->_i));
 328 }
 329
 330 integer si_abs(x)
 331 integer x;
 332 {
 333 if(x >= 0)
 334         return(x);
 335 return(- x);
 336 }
 337
 338 shortint sh_abs(x)
 339 shortint x;
 340 {
 341 if(x >= 0)
 342         return(x);
 343 return(- x);
 344 }