+#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 isPiOMPforCritical_Ok(int n);\r
+\r
+/*--------------------------------------*\\r
+ |* Private *|\r
+ \*-------------------------------------*/\r
+\r
+static double piOMPforCritique(int n);\r
+\r
+/*----------------------------------------------------------------------*\\r
+ |* Implementation *|\r
+ \*---------------------------------------------------------------------*/\r
+\r
+/*--------------------------------------*\\r
+ |* Public *|\r
+ \*-------------------------------------*/\r
+\r
+bool isPiOMPforCritical_Ok(int n)\r
+ {\r
+ return isAlgoPI_OK(piOMPforCritique, n, "Pi OMP for critique");\r
+ }\r
+\r
+/*--------------------------------------*\\r
+ |* Private *|\r
+ \*-------------------------------------*/\r
+\r
+/**\r
+ * synchronisation couteuse!\r
+ */\r
+double piOMPforCritique(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 critical(sum_fpi)\r
+ sum += fpi(xi);\r
+ }\r
+\r
+ return sum * DX;\r
+ }\r
+\r
+/*----------------------------------------------------------------------*\\r
+ |* End *|\r
+ \*---------------------------------------------------------------------*/\r
+\r