blatt 2 daxpy

blatt02/blatt02.md

+# Blatt 2
+## Aufgabe 1
+``` build/ARM/gem5.opt --debug-flags=O3PipeView --debug-start=1 --debug-file=trace.out configs/example/se.py --cpu-type=detailed --caches -c ~/advanced_computer_architecture/exercises/blatt02/exec/daxpy ```

blatt02/exec/Makefile

+
+all:
+	arm-linux-gnueabihf-gcc -static daxpy.c -o daxpy
+
+clean:
+	rm daxpy

blatt02/exec/Makefile~

+FILE1 = basicmath_small.c   rad2deg.c  cubic.c   isqrt.c
+FILE2 = basicmath_large.c   rad2deg.c  cubic.c   isqrt.c
+
+all: basicmath_small basicmath_large
+
+basicmath_small: ${FILE1} Makefile
+	arm-linux-gnueabihf-gcc -static -O3 ${FILE1} -o basicmath_small -lm
+basicmath_large: ${FILE2} Makefile
+	arm-linux-gnueabihf-gcc -static -O3 ${FILE2} -o basicmath_large -lm
+
+clean:
+	rm -rf basicmath_small basicmath_large output*

blatt02/exec/daxpy.c

+#include<stdio.h>
+#define N 8  // size of vector
+#define R 10 // number of repetitions
+
+
+// daxpy: Y = aX + Y (n size vector, a constant, X,Y vectors)
+void daxpy(double* x, double* y, double a, int n)
+{
+  int i;
+  for (i = 0; i < n; i++) {
+    y[i] = a * x[i] + y[i];
+  }
+}
+
+// pretty print vector
+void print_vector(double* vec, int n)
+{
+  int i;
+  for (i=0; i<n-1; i++) {
+    printf("%f,\t",vec[i]);
+  }
+  printf("%f",vec[n-1]);
+  printf("\n");
+}
+
+void main()
+{
+  static double x[N], y[N];
+  double a = 5.5;
+  int i;
+  
+  for (i = 0; i<N; i++) {
+    x[i] = 7.833667;
+    y[i] = 0;
+  }
+
+  print_vector(y, N);
+
+  for (i = 0; i<R; i++) {
+    daxpy(x,y,a,N);
+  }
+
+  print_vector(y,N);
+}