6 files changed, 1063 insertions, 0 deletions

diff --git a/baseline/source/fmref/Changelog.txt b/baseline/source/fmref/Changelog.txt
new file mode 100644
index 0000000..e9e7055
--- /dev/null
+++ b/baseline/source/fmref/Changelog.txt
@@ -0,0 +1,5 @@
+2017-06-27
+- Remove unused variables
+- Introduce fmref_init and fmref_return functions.
+- Add prefix fmref_ to global variables.
+- Introduce dummy initialization in ieee754_rem_pio2f to please linter.
diff --git a/baseline/source/fmref/fmref.c b/baseline/source/fmref/fmref.c
new file mode 100644
index 0000000..6be1436
--- /dev/null
+++ b/baseline/source/fmref/fmref.c
@@ -0,0 +1,283 @@
+/*
+* fmref.c: C reference implementation of FM Radio
+* David Maze <dmaze@cag.lcs.mit.edu>
+* $Id: fmref.c,v 1.2 2010-10-04 21:21:26 garus Exp $
+*/
+
+
+#include "../extra.h"
+#include "wcclibm.h"
+#ifndef M_PI
+#define M_PI        3.1415926535897932384626433832795
+#endif
+
+// Defines
+#define SAMPLING_RATE 250000000
+#define CUTOFF_FREQUENCY 108000000
+#define NUM_TAPS 64
+#define MAX_AMPLITUDE 27000.0
+#define BANDWIDTH 10000
+#define DECIMATION 4 
+/* Must be at least NUM_TAPS+1: */
+#define IN_BUFFER_LEN 200
+#define EQUALIZER_BANDS 10
+
+
+// Type declarations
+typedef struct FloatBuffer
+{
+  float buff[IN_BUFFER_LEN];
+  int rpos, rlen;
+}
+FloatBuffer;
+/* Low pass filter: */
+typedef struct LPFData
+{
+  float coeff[NUM_TAPS];
+  float freq;
+  int taps, decimation;
+}
+LPFData;
+typedef struct EqualizerData
+{
+  LPFData lpf[EQUALIZER_BANDS + 1];
+  FloatBuffer fb[EQUALIZER_BANDS + 1];
+  float gain[EQUALIZER_BANDS];
+}
+EqualizerData;
+
+// Global vars
+float fmref_lpf_coeff[NUM_TAPS];
+float fmref_eq_cutoffs[EQUALIZER_BANDS + 1] =
+  { 55.000004f, 77.78174f, 110.00001f, 155.56354f, 220.00002f, 311.12695f,
+    440.00003f, 622.25415f, 880.00006f, 1244.5078f, 1760.0001f };
+static int fmref_numiters = 2;
+
+// Forward declarations
+void fmref_fb_compact(FloatBuffer *fb);
+int fmref_fb_ensure_writable(FloatBuffer *fb, int amount);
+void fmref_get_floats(FloatBuffer *fb);
+void fmref_init_lpf_data(LPFData *data, float freq, int taps, int decimation);
+void fmref_run_lpf(FloatBuffer *fbin, FloatBuffer *fbout, LPFData *data);
+void fmref_run_demod(FloatBuffer *fbin, FloatBuffer *fbout);
+void fmref_init_equalizer(EqualizerData *data);
+void fmref_run_equalizer(FloatBuffer *fbin, FloatBuffer *fbout, EqualizerData *data);
+void fmref_main(void);
+
+void fmref_init(void)
+{
+  // dummy init function
+}
+
+int fmref_return(void)
+{
+  // dummy return value
+  return 0;
+}
+
+int main(int argc, char **argv){
+
+        SET_UP
+        for(jobsComplete=-1; jobsComplete<maxJobs; jobsComplete++){
+                START_LOOP
+                fmref_init();
+                fmref_main();
+                STOP_LOOP
+        }
+        WRITE_TO_FILE
+  return fmref_return();
+}
+
+FloatBuffer fmref_fb1, fmref_fb2, fmref_fb3, fmref_fb4;
+LPFData fmref_lpf_data;
+
+void fmref_main(void)
+{
+  int i;
+  EqualizerData eq_data;
+
+  fmref_fb1.rpos = fmref_fb1.rlen = 0;
+  fmref_fb2.rpos = fmref_fb2.rlen = 0;
+  fmref_fb3.rpos = fmref_fb3.rlen = 0;
+  fmref_fb4.rpos = fmref_fb4.rlen = 0;
+
+  fmref_init_lpf_data(&fmref_lpf_data, CUTOFF_FREQUENCY, NUM_TAPS, DECIMATION);
+  fmref_init_equalizer(&eq_data);
+
+  /* Startup: */
+  fmref_get_floats(&fmref_fb1);
+  /* LPF needs at least NUM_TAPS+1 inputs; fmref_get_floats is fine. */
+  fmref_run_lpf(&fmref_fb1, &fmref_fb2, &fmref_lpf_data);
+  /* run_demod needs 1 input, OK here. */
+  /* run_equalizer needs 51 inputs (same reason as for LPF).  This means
+   * running the pipeline up to demod 50 times in advance: */
+  _Pragma( "loopbound min 64 max 64" )
+  for (i = 0; i < 64; i++) {
+    if (fmref_fb1.rlen - fmref_fb1.rpos < NUM_TAPS + 1)
+      fmref_get_floats(&fmref_fb1);
+    fmref_run_lpf(&fmref_fb1, &fmref_fb2, &fmref_lpf_data);
+    fmref_run_demod(&fmref_fb2, &fmref_fb3);
+  }
+
+  /* Main loop: */
+  _Pragma( "loopbound min 2 max 2" )
+  while (fmref_numiters-- > 0) {
+    /* The low-pass filter will need NUM_TAPS+1 items; read them if we
+     * need to. */
+    if (fmref_fb1.rlen - fmref_fb1.rpos < NUM_TAPS + 1)
+      fmref_get_floats(&fmref_fb1);
+    fmref_run_lpf(&fmref_fb1, &fmref_fb2, &fmref_lpf_data);
+    fmref_run_demod(&fmref_fb2, &fmref_fb3);
+    fmref_run_equalizer(&fmref_fb3, &fmref_fb4, &eq_data);
+    
+  }
+}
+
+void fmref_fb_compact(FloatBuffer *fb)
+{
+  
+  int i;
+  char *source;
+  char *target;
+  target = (char*)(fb->buff);
+  source = (char*)(fb->buff + fb->rpos);
+  _Pragma( "loopbound min 0 max 60" )
+  for (i = 0; i < fb->rlen - fb->rpos; i++) {
+    target[i] = source[i];
+  }
+  fb->rlen -= fb->rpos;
+  fb->rpos = 0;
+}
+
+int fmref_fb_ensure_writable(FloatBuffer *fb, int amount)
+{
+  int available = IN_BUFFER_LEN - fb->rlen;
+  if (available >= amount)
+    return 1;
+
+  /* Nope, not enough room, move current contents back to the beginning. */
+  fmref_fb_compact(fb);
+
+  available = IN_BUFFER_LEN - fb->rlen;
+  if (available >= amount)
+    return 1;
+
+  return 0;
+}
+
+void fmref_get_floats(FloatBuffer *fb)
+{
+  static int x = 0;
+  fmref_fb_compact(fb);
+
+  /* Fill the remaining space in fb with 1.0. */
+  _Pragma( "loopbound min 200 max 200" )
+  while (fb->rlen < IN_BUFFER_LEN) {
+    fb->buff[fb->rlen++] = (float)x;
+    x++;
+  }
+}
+
+void fmref_init_lpf_data(LPFData *data, float freq, int taps, int decimation)
+{
+  /* Assume that CUTOFF_FREQUENCY is non-zero.  See comments in
+   * StreamIt LowPassFilter.java for origin. */
+  float w = 2 * M_PI * freq / SAMPLING_RATE;
+  int i;
+  float m = taps - 1.0f;
+
+  data->freq = freq;
+  data->taps = taps;
+  data->decimation = decimation;
+
+  _Pragma( "loopbound min 64 ma