Move the DIS benchmarks up a directory and update hardcoded paths

Note that this repo does not attempt to keep a copy of the original
DIS benchmark distributions. UNC real-time has another repo for that.

2 files changed, 343 insertions, 0 deletions

diff --git a/dis/Field/DISstressmarkRNG.h b/dis/Field/DISstressmarkRNG.h
new file mode 100644
index 0000000..4aa2620
--- /dev/null
+++ b/dis/Field/DISstressmarkRNG.h
@@ -0,0 +1,190 @@
+#include <math.h>
+
+#define IA 16807
+#define IM 2147483647
+#define AM (1.0/IM)
+#define IQ 127773
+#define IR 2836
+#define NTAB 32
+#define NDIV (1+(IM-1)/NTAB)
+#define EPS 1.2e-7
+#define RNMX (1.0-EPS)
+
+static long iy=0;
+static long iv[NTAB];
+static long iseed;
+
+int ABS(int x){
+  if (x>= 0) return x;
+  else
+    return (-x);
+}
+
+int sign(int x){
+  if (x >= 0) return 1;
+  else
+    return (-1);
+}
+
+int MAX(int x, int y){
+  if (x>= y) return x;
+  else
+    return y;
+}
+
+int MIN(int x, int y){
+  if (x<= y) return x;
+  else
+    return y;
+}
+
+void randInit(long idum)
+{
+  long j;
+  long k;
+
+  assert (idum <= 0);
+  assert (iy == 0);
+
+  iseed = idum;
+  if (-(iseed)<1){
+    iseed = 1;
+  }
+  else {
+    iseed = -(iseed);
+  }
+  for (j=NTAB+7; j>=0; j--){
+    k = (iseed)/IQ;
+    iseed = IA*(iseed-k*IQ)-IR*k;
+    if (iseed < 0){
+      iseed += IM;
+    }
+    if (j < NTAB){
+      iv[j] = iseed;
+    }
+  }
+  iy = iv[0];
+}
+
+float randNum()
+{
+  long j;
+  long k;
+  float temp;
+
+  assert (iy != 0);
+  
+  k = (iseed)/IQ;
+  iseed = IA*(iseed-k*IQ)-IR*k;
+
+  if (iseed < 0){
+    iseed += IM;
+  }
+  j = iy/NDIV; 
+  iy = iv[j];
+  iv[j] = iseed;
+
+  temp = AM * iy;
+
+  if (temp > RNMX){
+    return RNMX;
+  }
+  else {
+    return temp;
+  }
+}
+
+
+float randomFloat(float lowest_float, float highest_float)
+{
+ float value;
+ float range;
+
+assert (lowest_float < highest_float);
+
+range = highest_float - lowest_float;
+value = randNum()*(highest_float - lowest_float) + lowest_float;
+assert(value >= lowest_float);
+assert(value <= highest_float);
+
+return value;
+
+}
+
+float randomNonZeroFloat(float lowest_float, float highest_float, float epsilon)
+{
+
+  double range;
+  float value;
+
+
+  assert (lowest_float < 0);
+  assert (highest_float > 0);
+  assert (epsilon > 0);
+  assert ((epsilon < -lowest_float) && (epsilon < highest_float));
+  
+  range = highest_float - lowest_float;
+  value = (randNum() * range)+lowest_float;
+  
+  if (ABS(value) < epsilon) 
+    {
+      if (value > 0) value = value + epsilon;
+      else if (value < 0) value = value - epsilon;
+      
+    }
+
+  assert (value >= lowest_float);
+  assert (value <= highest_float);
+
+  return value;
+}
+
+unsigned int randomUInt(int lowest_uint, int highest_uint)
+{
+  float range;
+  unsigned int value;
+  float temp;
+
+  range =(float)(highest_uint - lowest_uint + 1);
+  temp = randNum();
+  value =(unsigned int)( floor(temp * range) + lowest_uint); 
+
+  assert (value >= lowest_uint);
+  assert (value <= highest_uint);
+
+  return value;
+}
+
+unsigned int randomNonZeroUInt(int lowest_uint, int highest_uint)
+{
+  float range;
+  unsigned int value;
+  float temp;
+
+  range =(float)(highest_uint - lowest_uint + 1);
+  value = 0;
+  while(value == 0){
+  temp = randNum();
+
+  value =(unsigned int)( floor(temp * range) + lowest_uint); 
+  }
+  
+  assert (value >= lowest_uint);
+  assert (value <= highest_uint);
+
+  return value;
+}
+
+int randInt(int lowest_uint,  int highest_uint)
+{
+  float range;
+  int value;
+  
+  range = highest_uint - lowest_uint + 1;
+  value = (int)(floor(randNum() * range) + lowest_uint);
+  
+  assert (value >= lowest_uint);
+  assert (value <= highest_uint);
+
+  return value;
+}                            
diff --git a/dis/Field/field.c b/dis/Field/field.c
new file mode 100644
index 0000000..8565e8c
--- /dev/null
+++ b/dis/Field/field.c
@@ -0,0 +1,153 @@
+/*
+ *  Sample code for the DIS Field Stressmark
+ *
+ * This source code is the completely correct source code based on 
+ * the example codes provided by Atlantic Aerospace Division, Titan 
+ * Systems Corporation, 2000.
+ * 
+ * If you just compile and generate the executables from this source 
+ * code, this code would be enough. However, if you wish to get a complete 
+ * understanding of this stressmark, it is strongly suggested that you
+ * read the Benchmark Analysis and Specifications Document Version 1.0
+ * before going on since the detailed comments are given in this documents.
+ * the comments are not repeated here.
+ */
+
+#include <stdio.h>
+#include <time.h>
+#include <assert.h>
+#include <stdlib.h>
+#include "DISstressmarkRNG.h"
+#include "extra.h"
+
+
+#define MIN_FIELD_SIZE 16
+#define MAX_FIELD_SIZE (16777216*4) // This has been quadrupled from original
+#define MIN_SEED -2147483647
+#define MAX_SEED -1
+#define MIN_MOD_OFFSET 0
+#define MAX_MOD_OFFSET 65535
+#define MIN_TOKENS 1
+#define MAX_TOKENS 256
+#define MIN_TOKEN_LENGTH 1
+#define MAX_TOKEN_LENGTH 8
+#define MIN_TOKEN_VALUE 0
+#define MAX_TOKEN_VALUE 255
+#define MAX_SUBFIELDS 256
+
+/*
+ * main()
+ */
+
+int main(int argc, char** argv){
+  SET_UP
+  unsigned char *field;
+  unsigned int f;
+  int seed;
+  int mod_offset;
+  unsigned int n;
+
+  time_t startTime;
+ 
+  struct tokenS{
+  unsigned char delimiter[MAX_TOKEN_LENGTH];
+  unsigned char length; 
+    struct statisticS{
+      unsigned int count;
+      unsigned char min;
+      unsigned char sum;
+    }stat[MAX_SUBFIELDS];
+    unsigned char subfields;
+  }token[MAX_TOKENS];
+
+  unsigned int l;
+
+  assert(fscanf(stdin, "%d %d %d %d", &f, &seed, &mod_offset, &n) == 4);
+
+  assert((f >= MIN_FIELD_SIZE) && (f <= MAX_FIELD_SIZE));
+  assert((seed >= MIN_SEED) && (seed <= MAX_SEED));
+  assert((mod_offset >= MIN_MOD_OFFSET) && (mod_offset <= MAX_MOD_OFFSET));
+  assert((n >= MIN_TOKENS) && (n <= MAX_TOKENS));
+    for (l=0; l<n; l++){
+      int x;
+      int index;
+      index = 0;
+      assert(fscanf(stdin,"%x", &x) == 1);
+      while(x!=0){
+        assert((x >= MIN_TOKEN_VALUE) && (x <= MAX_TOKEN_VALUE));
+          token[l].delimiter[index] = (unsigned char )x;
+        index++;
+        assert(fscanf(stdin,"%x", &x) == 1);
+      }
+      assert((index >= MIN_TOKEN_LENGTH) && (index <= MAX_TOKEN_LENGTH));
+        token[l].length = index;
+    }
+
+  if ((field = (unsigned char*)malloc(f*sizeof(unsigned char))) == NULL)
+    return (-1);
+
+  randInit(seed);
+  for (l =0; l<f; l++){
+    field[l] = randInt(MIN_TOKEN_VALUE, MAX_TOKEN_VALUE);
+  }
+  
+  startTime = time(NULL);
+  
+  for (l =0; l<n; l++){
+    START_LOOP
+    unsigned int index;
+    
+    token[l].subfields = 0;
+    token[l].stat[0].count = 0;
+    token[l].stat[0].sum = 0;
+    token[l].stat[0].min = MAX_TOKEN_VALUE;
+
+    index = 0;
+    while ((index < f) && (token[l].subfields < MAX_SUBFIELDS)){
+      unsigned char offset;
+      offset = 0;
+      while ((field[index+offset] == token[l].delimiter[offset]) &&
+             (offset < token[l].length)){
+        offset++;
+      }
+
+      if (offset == token[l].length){
+        for (offset=0; offset<token[l].length; offset++){
+          field[index+offset] = (field[index+offset] + 
+                                 field[(index+offset+mod_offset) % f])
+                                 %(MAX_TOKEN_VALUE+1);
+        }
+        index += token[l].length-1;
+        token[l].subfields++;
+        token[l].stat[token[l].subfields].count = 0;
+        token[l].stat[token[l].subfields].sum = 0;
+        token[l].stat[token[l].subfields].min = MAX_TOKEN_VALUE;
+      }                        
+                    
+        else{
+        token[l].stat[token[l].subfields].count++;
+        token[l].stat[token[l].subfields].sum += field[index];
+        if (token[l].stat[token[l].subfields].min > field[index])
+        token[l].stat[token[l].subfields].min = field[index];
+        }            
+        index++;
+    }
+    token[l].subfields++;
+    STOP_LOOP
+    }
+
+    startTime = time(NULL) - startTime;
+
+    for (l = 0; l< n; l++){
+      unsigned int ll;
+      fprintf(stdout, "%d subfields for token %d \n", token[l].subfields, l);
+      for ( ll =0; ll<token[l].subfields; ll++)
+        fprintf(stdout, "subfields %d:\tcount = %d\tmin= %x\tsum= %x\n",
+                ll, token[l].stat[ll].count, 
+                token[l].stat[ll].min, token[l].stat[ll].sum);
+    }
+    fprintf(stdout, "total time = %d seconds.\n", (int)startTime);
+    free(field);
+    WRITE_TO_FILE
+    return(0);
+  }