update problem sizes for examples/

Some of the problem sizes were too small (checked-in earlier inadvertently). Data set sizes are now reasonably large (don't fit in a typical L3 cache) Signed-off-by: Uday Bondhugula <udayreddy@gmail.com>
tajkhan · Jun 22, 2012 · f7beb37 · f7beb37
1 parent 6de9fc9
commit f7beb37
Show file tree

Hide file tree

Showing 16 changed files with 105 additions and 92 deletions.
diff --git a/examples/covcol/covcol.c b/examples/covcol/covcol.c
@@ -7,8 +7,8 @@
 
 #define my_sqrt_array(x,j) sqrt(x[j])
 
-#define N 1000
-#define M 1000
+#define N 1500
+#define M 1500
 
 float float_n = (float) N;
 float eps = 0.005;

diff --git a/examples/dsyr2k/Makefile b/examples/dsyr2k/Makefile
@@ -1,6 +1,6 @@
 SRC= dsyr2k
 
-PLCFLAGS = --context=100 #--unroll --ufactor=2
+PLCFLAGS = --unroll --ufactor=4
 TILEFLAGS = #--l2tile
 
 include ../common.mk

diff --git a/examples/dsyr2k/dsyr2k.c b/examples/dsyr2k/dsyr2k.c
@@ -3,7 +3,7 @@
 #include <math.h>
 #include <sys/time.h>
 
-#define NMAX 4096
+#define NMAX 2048
 
 #pragma declarations
 double a[NMAX][NMAX], b[NMAX][NMAX], c[NMAX][NMAX];

diff --git a/examples/dsyrk/dsyrk.c b/examples/dsyrk/dsyrk.c
@@ -5,7 +5,7 @@
 #include <unistd.h>
 #include <sys/time.h>
 
-#define NMAX 3000
+#define NMAX 2000
 #pragma declarations
 double a[NMAX][NMAX], c[NMAX][NMAX];
 #pragma enddeclarations

diff --git a/examples/fdtd-2d/Makefile b/examples/fdtd-2d/Makefile
@@ -1,6 +1,6 @@
 SRC=fdtd-2d
 
-PLCFLAGS += #--context=32  #--unroll
+PLCFLAGS += #--context=32  --unroll
 TILEFLAGS = #--l2tile
 
 include ../common.mk

diff --git a/examples/fdtd-2d/fdtd-2d.c b/examples/fdtd-2d/fdtd-2d.c
@@ -10,8 +10,8 @@
 #endif
 
 #define tmax 128
-#define nx 4000
-#define ny 4000
+#define nx 3000
+#define ny 3000
 
 #pragma declarations
 double ex[nx][ny+1];

diff --git a/examples/jacobi-1d-imper/Makefile b/examples/jacobi-1d-imper/Makefile
@@ -9,12 +9,12 @@ include ../common.mk
 $(SRC).sched.c: 
 	$(PLC) $(SRC).c --sched --tile --parallel $(TILEFLAGS) $(PLCFLAGS) -o $(SRC).sched.c
 
-sched: $(SRC).sched.c decls.h papi_defs.h util.h
+sched: $(SRC).sched.c papi_defs.h util.h
 	$(CC) $(OPT_FLAGS) -openmp -lm $(SRC).sched.c -o sched
 
-sched_test: $(SRC).sched.c decls.h papi_defs.h util.h
+sched_test: $(SRC).sched.c papi_defs.h util.h
 	$(CC) $(OPT_FLAGS) -openmp -lm $(SRC).sched.c -o sched_test -DTEST
 
-limlam: $(SRC).limlam.c decls.h papi_defs.h util.h
+limlam: $(SRC).limlam.c papi_defs.h util.h
 	$(CC) $(OPT_FLAGS) -openmp -lm $(SRC).limlam.c -o limlam
 
diff --git a/examples/jacobi-1d-imper/jacobi-1d-imper.c b/examples/jacobi-1d-imper/jacobi-1d-imper.c
@@ -4,8 +4,8 @@
 
 #include "util.h"
 
-#define N 1000
-#define T 50
+#define N 2000000
+#define T 1000
 
 #pragma declarations
 double a[N];

diff --git a/examples/jacobi-1d-imper/jacobi-1d-imper.limlam.c b/examples/jacobi-1d-imper/jacobi-1d-imper.limlam.c
@@ -3,7 +3,14 @@
 #include <math.h>
 
 #include <assert.h>
-#include "decls.h"
+
+#define N 2000000
+#define T 1000
+
+#pragma declarations
+double a[N];
+double b[N];
+#pragma enddeclarations
 
 #ifdef PERFCTR
 #include <papi.h>

diff --git a/examples/jacobi-1d-imper/jacobi-1d-imper.sched.c b/examples/jacobi-1d-imper/jacobi-1d-imper.sched.c
@@ -8,9 +8,37 @@
 #include <papi.h>
 #endif
 
-#include "decls.h"
+#define N 2000000
+#define T 1000
+
+#pragma declarations
+double a[N];
+double b[N];
+#pragma enddeclarations
+
+double t_start, t_end;
+
+void init_array()
+{
+    int j;
+
+    for (j=0; j<N; j++) {
+        a[j] = ((double)j)/N;
+    }
+}
+
+
+void print_array()
+{
+    int j;
+
+    for (j=0; j<N; j++) {
+        fprintf(stdout, "%lf ", a[j]);
+        if (j%80 == 20) fprintf(stdout, "\n");
+    }
+    fprintf(stdout, "\n");
+}
 
-#include "util.h"
 
 #define ceild(n,d)  ceil(((double)(n))/((double)(d)))
 #define floord(n,d) floor(((double)(n))/((double)(d)))

diff --git a/examples/jacobi-2d-imper/jacobi-2d-imper.c b/examples/jacobi-2d-imper/jacobi-2d-imper.c
@@ -5,7 +5,7 @@
 #include <assert.h>
 
 #define N 2000
-#define T 2000
+#define T 1000
 
 #pragma declarations
 double a[N][N];

diff --git a/examples/lu/lu.c b/examples/lu/lu.c
@@ -6,16 +6,62 @@
 
 #include <omp.h>
 
-#define N 1024
+#define N 2000
 
 #pragma declarations
-double a[N][N+13];
+double a[N][N];
 //double v_a[32][35];
 //double v_b[32][32];
 //double v_c[32][33];
 #pragma enddeclarations
 
-#include "util.h"
+#include <unistd.h>
+#include <sys/time.h>
+#include <math.h>
+
+#ifdef TIME
+#define IF_TIME(foo) foo;
+#else
+#define IF_TIME(foo)
+#endif
+
+void init_array()
+{
+    int i, j, k;
+
+    for (i=0; i<N; i++) {
+        for (j=0; j<N; j++) {
+            for (k=0; k<N; k++) {
+                a[i][j] += (i+k+1)*(k+j+1);//i==j?1:0;
+            }
+        }
+    }
+}
+
+
+void print_array()
+{
+    int i, j;
+
+    for (i=0; i<N; i++) {
+        for (j=0; j<N; j++) {
+            fprintf(stdout, "%lf ", round(a[i][j]));
+            if (j%80 == 79) fprintf(stdout, "\n");
+        }
+        fprintf(stdout, "\n");
+    }
+    fprintf(stdout, "\n");
+}
+
+double rtclock()
+{
+    struct timezone Tzp;
+    struct timeval Tp;
+    int stat;
+    stat = gettimeofday (&Tp, &Tzp);
+    if (stat != 0) printf("Error return from gettimeofday: %d",stat);
+    return(Tp.tv_sec + Tp.tv_usec*1.0e-6);
+}
 
 int main()
 {

diff --git a/examples/lu/util.h b/examples/lu/util.h
diff --git a/examples/seidel/seidel.c b/examples/seidel/seidel.c
@@ -4,8 +4,8 @@
 
 #include <assert.h>
 
-#define N 1000
-#define T 500
+#define N 2000
+#define T 200
 
 #pragma declarations
 double a[N][N];

diff --git a/examples/ssymm/ssymm.c b/examples/ssymm/ssymm.c
@@ -3,7 +3,7 @@
 #include <math.h>
 #include <sys/time.h>
 
-#define NMAX 500
+#define NMAX 2000
 
 #pragma declarations
 double a[NMAX][NMAX], b[NMAX][NMAX], c[NMAX][NMAX];

diff --git a/examples/strmm/strmm.c b/examples/strmm/strmm.c
@@ -3,7 +3,7 @@
 #include <math.h>
 #include <sys/time.h>
 
-#define NMAX 1000
+#define NMAX 2000
 
 #pragma declarations
 double a[NMAX][NMAX], b[NMAX][NMAX], c[NMAX][NMAX];