--- /dev/null
+#include <omp.h>\r
+#include "00_pi_tools.h"\r
+#include "OmpTools.h"\r
+\r
+/*----------------------------------------------------------------------*\\r
+ |* Declaration *|\r
+ \*---------------------------------------------------------------------*/\r
+\r
+/*--------------------------------------*\\r
+ |* Imported *|\r
+ \*-------------------------------------*/\r
+\r
+\r
+\r
+/*--------------------------------------*\\r
+ |* Public *|\r
+ \*-------------------------------------*/\r
+\r
+bool isPiOMPforAtomic_Ok(int n);\r
+\r
+/*--------------------------------------*\\r
+ |* Private *|\r
+ \*-------------------------------------*/\r
+\r
+static double piOMPforAtomic(int n);\r
+\r
+/*----------------------------------------------------------------------*\\r
+ |* Implementation *|\r
+ \*---------------------------------------------------------------------*/\r
+\r
+/*--------------------------------------*\\r
+ |* Public *|\r
+ \*-------------------------------------*/\r
+\r
+bool isPiOMPforAtomic_Ok(int n)\r
+ {\r
+ return isAlgoPI_OK(piOMPforAtomic, n, "Pi OMP for atomic");\r
+ }\r
+\r
+/*--------------------------------------*\\r
+ |* Private *|\r
+ \*-------------------------------------*/\r
+\r
+/**\r
+ * synchronisation couteuse!\r
+ */\r
+double piOMPforAtomic(int n)\r
+ {\r
+ double sum = 0;\r
+ double xi;\r
+ const double DX = 1.0/(double)n;\r
+\r
+ #pragma omp parralel for private(xi)\r
+ for (int i = 0; i < n; i++)\r
+ {\r
+ xi = i*DX;\r
+ #pragma omp atomic\r
+ sum += fpi(xi); // Uniquement l'operateur += est atomi, 'fpi' est exécuté en parralèle.\r
+ }\r
+\r
+ return sum * DX;\r
+ }\r
+\r
+/*----------------------------------------------------------------------*\\r
+ |* End *|\r
+ \*---------------------------------------------------------------------*/\r
+\r