Ajout TP2

TPs/TP2/CODE/Makefile

+CC=gcc
+CFLAGS=-O3 -Wall
+LDFLAGS=-lm
+EXE=tp2.exe
+
+all : $(EXE)
+
+$(EXE) : tp2.o
+	$(CC) $(CFLAGS) -o $@ $< $(LDFLAGS) 
+
+%.o : %.c
+	$(CC) $(CFLAGS) -c -o $@ $<
+
+run : $(EXE)
+	./$< 10240
+
+proper : 
+	rm -f *.o
+
+clean :
+	rm -f *.o $(EXE)

TPs/TP2/CODE/tp2.c

+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define NSTREAMS 4
+
+void verif(float *out, int sz)
+{
+  float err = 0.;
+
+  for(int i = 0 ; i < sz ; i++)
+  {
+    err += abs(out[i] - exp( - abs(sin(i * 1.0)) ));
+  }
+
+  if (err/sz < 1.e-4)
+  {
+	  fprintf(stdout, "TEST PASSED (error %3.f < 1.e-4)\n", err/sz);
+  }
+  else
+  {
+	  fprintf(stderr, "TEST FAILED (error %3.f > 1.e-4)\n", err/sz);
+  }
+}
+
+void func(float *out, int size)
+{
+  for(int i = 0; i < size; ++i)
+  {
+    out[i] = exp( - abs(out[i]) );
+  }
+}
+
+int main(int argc, char** argv)
+{
+  int size = 1024;
+  if (argc == 2)
+  {
+	  size = atoi(argv[1]);
+  }
+
+  size *= NSTREAMS;
+
+  float *tab = NULL;
+  tab = (float*) malloc(sizeof(float) * size);
+
+  if(tab == NULL)
+  {
+    fprintf(stderr, "Bad allocation\n");
+    return -1;
+  }
+
+  for(int i = 0; i < size; ++i)
+  {
+    tab[i] = sin(i * 1.);
+  }
+
+  func(tab, size);
+
+  verif(tab, size);
+
+  free(tab);
+  return 0;
+}

TPs/TP2/CODE/tp2.cu

+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define NSTREAMS 4
+
+void verif(float *out, int sz)
+{
+  float err = 0.;
+
+  for(int i = 0 ; i < sz ; i++)
+  {
+    err += abs(out[i] - exp( - abs(sin(i * 1.0)) ));
+  }
+
+  if (err/sz < 1.e-4)
+  {
+	  fprintf(stdout, "TEST PASSED (error %3.f < 1.e-4)\n", err/sz);
+  }
+  else
+  {
+	  fprintf(stderr, "TEST FAILED (error %3.f > 1.e-4)\n", err/sz);
+  }
+}
+
+void func(float *out, int size)
+{
+  for(int i = 0; i < size; ++i)
+  {
+    out[i] = exp( - abs(out[i]) );
+  }
+}
+
+int main(int argc, char** argv)
+{
+  int size = 1024;
+  if (argc == 2)
+  {
+	  size = atoi(argv[1]);
+  }
+
+  size *= NSTREAMS;
+
+  float *tab = NULL;
+  tab = (float*) malloc(sizeof(float) * size);
+
+  float d_tab;
+  cudaMalloc(
+
+  if(tab == NULL)
+  {
+    fprintf(stderr, "Bad allocation\n");
+    return -1;
+  }
+
+  for(int i = 0; i < size; ++i)
+  {
+    tab[i] = sin(i * 1.);
+  }
+
+  func(tab, size);
+
+  verif(tab, size);
+
+  free(tab);
+  return 0;
+}