Chapel

CHAPEL/Stencil/Makefile

@@ -0,0 +1,25 @@
+# Chapel Makefile
+
+# TODO Find a way to make this fit into the existing Makefile system
+
+OPTFLAGS = --fast -sassertNoSlicing
+MODULEFLAGS = --module-dir ../../include/
+
+# Debugging
+DEBUGFLAGS =
+ifdef DEBUG
+DEBUGFLAGS = --savec c
+endif
+
+CHPL = chpl
+
+SOURCES = stencil-serial.chpl stencil-shared.chpl
+EXECUTABLE = $(SOURCES:.chpl=)
+
+all: $(EXECUTABLE)
+
+%: %.chpl
+	$(CHPL) $(OPTFLAGS) $(MODULEFLAGS) $(DEBUGFLAGS) $*.chpl -o $*
+
+clean:
+	rm -f $(EXECUTABLE)

CHAPEL/Stencil/stencil-serial.chpl

@@ -0,0 +1,102 @@
+// Chapel's serial stencil implementation
+use PRK;
+
+// Note: Defaulting to STAR stencil (defines weight)
+// Configurable runtime constants
+config const n: int = 100;
+config const iterations: int = 100;
+config const debug: bool = false;
+
+// Compile constants
+param R = 2;
+param coefx = 1.0;
+param coefy = 1.0;
+param epsilon = 1.e-8;
+
+// Runtime constants
+const activePoints = (n-2*R)*(n-2*R);
+const stencilSize = 4*R + 1;
+
+// Timer
+var timer: Timer;
+
+// Domains
+const Dom = {0.. # n, 0.. # n},
+      InnerDom = Dom.expand(-R),
+       W = {-R..R, -R..R};
+
+// Arrays
+var input, output: [Dom] real = 0.0;
+var weight: [W] real = 0.0;
+
+for i in 1..R do {
+  const element = 1.0 / (2.0*i*R);
+  weight[0, i]  =  element;
+  weight[i, 0]  =  element;
+  weight[-i, 0] = -element;
+  weight[0, -i] = -element;
+}
+
+// Initialize the input and output arrays
+serial do [(i, j) in Dom] input[i,j] = coefx*i+coefy*j;
+
+// Print information before main loop
+writeln("Parallel Research Kernels Version ", PRKVERSION);
+writeln("Serial stencil execution on 2D grid");
+writeln("Grid size            = ", n);
+writeln("Radius of stencil    = ", R);
+writeln("Type of stencil      = star"); // Temporarily hard-coded
+writeln("Data type            = double precision");
+writeln("Untiled");                     // Temporarily hard-coded
+writeln("Number of iterations = ", iterations);
+
+for iteration in 0..iterations do {
+  // Start timer after warmup iteration
+  if (iteration == 1) {
+    timer.start();
+  }
+
+  for (i,j) in InnerDom {
+    for jj in -R..R  do output[i, j] += weight[0, jj] * input[i, j+jj];
+    for ii in -R..-1 do output[i, j] += weight[ii, 0] * input[i+ii, j];
+    for ii in 1..R   do output[i, j] += weight[ii, 0] * input[i+ii, j];
+  }
+
+
+  // Add constant to solution to force refresh of neighbor data, if any
+  for (i,j) in Dom do input[i,j] += 1.0;
+
+} // end of iterations
+
+timer.stop();
+
+// Timings
+var stencilTime = timer.elapsed();
+writeln("stencil_time: ", stencilTime);
+
+// Compute L1 norm in parallel
+var norm = + reduce abs(output);
+
+norm /= activePoints;
+
+/*******************************************************************************
+** Analyze and output results.
+********************************************************************************/
+
+// Verify correctness
+var referenceNorm = (iterations + 1) * (coefx + coefy);
+
+if abs(norm-referenceNorm) > epsilon then {
+  writeln("ERROR: L1 norm = ", norm, ", Reference L1 norm = ", referenceNorm);
+  exit(1);
+} else {
+  writeln("Solution validates");
+  if debug then {
+    writeln("L1 norm = ", norm, ", Reference L1 norm = ", referenceNorm);
+  }
+}
+
+var flops = (2*stencilSize + 1) * activePoints;
+var avgTime = stencilTime / iterations;
+writeln("Rate (MFlops/s): ", 1.0E-06 * flops/avgTime,
+        "  Avg time (s): ", avgTime);

CHAPEL/Stencil/stencil-shared.chpl

@@ -0,0 +1,105 @@
+// Chapel's shared-memory parallel stencil
+use PRK;
+
+// Note: Defaulting to STAR stencil (defines weight)
+// Configurable runtime constants
+config const n: int = 100;
+config const iterations: int = 100;
+config const debug: bool = false;
+
+// Compile constants
+param R = 2;
+param coefx = 1.0;
+param coefy = 1.0;
+param epsilon = 1.e-8;
+
+// Runtime constants
+const activePoints = (n-2*R)*(n-2*R);
+const stencilSize = 4*R + 1;
+
+// Timer
+var timer: Timer;
+
+// Domains
+const    Dom = {0.. # n, 0.. # n},
+ProblemSpace = Dom,
+    InnerDom = ProblemSpace.expand(-R),
+           W = {-R..R, -R..R};
+
+// Arrays (initialized to zeros)
+var input, output:  [ProblemSpace] real = 0.0;
+var weight: [W] real = 0.0;
+
+forall i in 1..R do {
+  const element = 1.0 / (2.0*i*R);
+  weight[0, i]  =  element;
+  weight[i, 0]  =  element;
+  weight[-i, 0] = -element;
+  weight[0, -i] = -element;
+}
+
+// Initialize the input and output arrays
+[(i, j) in ProblemSpace] input[i,j] = coefx*i+coefy*j;
+
+// Print information before main loop
+writeln("Parallel Research Kernels Version ", PRKVERSION);
+writeln("Serial stencil execution on 2D grid");
+writeln("Grid size            = ", n);
+writeln("Radius of stencil    = ", R);
+writeln("Type of stencil      = star"); // Temporarily hard-coded
+writeln("Data type            = double precision");
+writeln("Untiled");                     // Temporarily hard-coded
+writeln("Number of iterations = ", iterations);
+
+for iteration in 0..iterations {
+
+  // Start timer after warmup iteration
+  if (iteration == 1) {
+    timer.start();
+  }
+
+  forall (i,j) in ProblemSpace {
+    for jj in -R..R  do output[i, j] += weight[0, jj] * input[i, j+jj];
+    for ii in -R..-1 do output[i, j] += weight[ii, 0] * input[i+ii, j];
+    for ii in 1..R   do output[i, j] += weight[ii, 0] * input[i+ii, j];
+  }
+
+  // Add constant to solution to force refresh of neighbor data, if any
+  forall (i,j) in ProblemSpace {
+    input[i, j] += 1.0;
+  }
+
+} // end of iterations
+
+timer.stop();
+
+// Timings
+var stencilTime = timer.elapsed();
+writeln("stencil_time: ", stencilTime);
+
+// Compute L1 norm in parallel
+var norm = + reduce abs(output);
+
+norm /= activePoints;
+
+/*******************************************************************************
+** Analyze and output results.
+********************************************************************************/
+
+// Verify correctness
+var referenceNorm = (iterations + 1) * (coefx + coefy);
+
+if (abs(norm-referenceNorm) > epsilon) {
+  writeln("ERROR: L1 norm = ", norm, ", Reference L1 norm = ", referenceNorm);
+  exit(1);
+} else {
+  writeln("Solution validates");
+  if debug {
+    writeln("L1 norm = ", norm, ", Reference L1 norm = ", referenceNorm);
+  }
+}
+
+var flops = (2*stencilSize + 1) * activePoints;
+var avgTime = stencilTime / iterations;
+writeln("Rate (MFlops/s): ", 1.0E-06 * flops/avgTime,
+        "  Avg time (s): ", avgTime);

CHAPEL/Synch_p2p/Makefile

@@ -0,0 +1,25 @@
+# Chapel Makefile
+
+# TODO Find a way to make this fit into the existing Makefile system
+
+OPTFLAGS = --fast -sassertNoSlicing
+MODULEFLAGS = --module-dir ../../include/
+
+# Debugging
+DEBUGFLAGS =
+ifdef DEBUG
+DEBUGFLAGS = --savec c
+endif
+
+CHPL = chpl
+
+SOURCES = p2p-serial.chpl
+EXECUTABLE = $(SOURCES:.chpl=)
+
+all: $(EXECUTABLE)
+
+%: %.chpl
+	$(CHPL) $(OPTFLAGS) $(MODULEFLAGS) $(DEBUGFLAGS) $*.chpl -o $*
+
+clean:
+	rm -f $(EXECUTABLE)

CHAPEL/Synch_p2p/p2p-serial.chpl

@@ -0,0 +1,78 @@
+// Chapel's serial implementation of synch_p2p
+
+use PRK;
+
+config const iterations = 1000;
+config const m = 1000;
+config const n = 1000;
+config const debug: bool = false;
+
+param timer: Timer;
+
+if (iterations < 1) {
+  writeln("ERROR: iterations must be >= 1 : ", iterations);
+  exit(1);
+}
+
+if (m < 1 || n < 1) {
+  writeln("ERROR: grid dimensions must be positive:", m, ", ", n);
+  exit(1);
+}
+
+// Initialize and zero out vector
+const mrange = 0 .. # m,
+      nrange = 0 .. # n,
+      Dom = {mrange, nrange},
+      Dom1 = {1 ..(m-1), 1 .. (n-1)};
+
+var vector : [Dom] real = 0.0;
+
+writeln("Parallel Research Kernels version ", PRKVERSION);
+writeln("Serial pipeline execution on 2D grid");
+writeln("Grid sizes                = ", m, ", ", n);
+writeln("Number of iterations      = ", iterations);
+
+// Set boundary values (top and left side of grid)
+serial do [j in nrange] vector[0,j] = j;
+serial do [i in mrange] vector[i,0] = i;
+
+for iteration in 0 .. iterations {
+
+  // Start timer after warmup iteration
+  if (iteration == 1) then timer.start();
+
+  for (i,j) in Dom1 {
+    vector[i, j] = vector[i-1, j] + vector[i, j-1] - vector[i-1, j-1];
+  }
+
+  // Copy bottom right corner value to top left, creating dependency
+  vector[0, 0] = -vector[m-1, n-1];
+}
+
+timer.stop();
+
+var pipelineTime = timer.elapsed();
+
+// Analyze and output results
+
+// Error threshold
+const epsilon = 1.e-8;
+
+// Verify correctness, using bottom right value of array
+var cornerValue : real = (iterations + 1)*(n + m - 2);
+if (abs(vector[m-1, n-1] - cornerValue) / cornerValue > epsilon) {
+  writeln("ERROR: checksum ", vector[m-1, n-1], " does not match verification \
+      value", cornerValue);
+  exit(1);
+}
+
+if (debug) {
+  writeln("Solution validates; verification value = ", cornerValue);
+} else {
+  writeln("Solution validates");
+}
+
+var avgTime = pipelineTime / iterations;
+
+writeln("Rate (MFlops/s): ", 1.0e-6 * 2 * ((m-1)*(n-1)) / avgTime,
+        " Avg time (s): ", avgTime);

CHAPEL/Transpose/Makefile

@@ -0,0 +1,25 @@
+# Chapel Makefile
+
+# TODO Find a way to make this fit into the existing Makefile system
+
+OPTFLAGS = --fast -sassertNoSlicing
+MODULEFLAGS = --module-dir ../../include/
+
+# Debugging
+DEBUGFLAGS =
+ifdef DEBUG
+DEBUGFLAGS = --savec c
+endif
+
+CHPL = chpl
+
+SOURCES = transpose-serial.chpl transpose-shared.chpl
+EXECUTABLE = $(SOURCES:.chpl=)
+
+all: $(EXECUTABLE)
+
+%: %.chpl
+	$(CHPL) $(OPTFLAGS) $(MODULEFLAGS) $(DEBUGFLAGS) $*.chpl -o $*
+
+clean:
+	rm -f $(EXECUTABLE)