diff --git a/eo/test/CMakeLists.txt b/eo/test/CMakeLists.txt
index fe28b21f5..8d7956bab 100644
--- a/eo/test/CMakeLists.txt
+++ b/eo/test/CMakeLists.txt
@@ -93,6 +93,18 @@ ELSEIF(ENABLE_CMAKE_TESTING)
     INSTALL(TARGETS ${test} RUNTIME DESTINATION share/eo/test COMPONENT test)
   ENDFOREACH (test)
 
+  SET(RESOURCES
+    boxplot.py
+    )
+
+  FOREACH(file ${RESOURCES})
+    EXECUTE_PROCESS(
+      COMMAND ${CMAKE_COMMAND} -E copy_if_different
+      ${CMAKE_CURRENT_SOURCE_DIR}/${file}
+      ${CMAKE_CURRENT_BINARY_DIR}/${file}
+      )
+  ENDFOREACH(file)
+
 ENDIF(ENABLE_MINIMAL_CMAKE_TESTING)
 
 ######################################################################################
diff --git a/eo/test/boxplot.py b/eo/test/boxplot.py
new file mode 100755
index 000000000..0590e4e97
--- /dev/null
+++ b/eo/test/boxplot.py
@@ -0,0 +1,10 @@
+#!/usr/bin/env python
+
+from pylab import *
+import sys
+
+DEFAULT_RESULTS_NAME = 'results.txt'
+
+if __name__ == '__main__':
+    boxplot( [ [ float(value) for value in line.split() ] for line in open( DEFAULT_RESULTS_NAME if len(sys.argv) < 2 else sys.argv[1] ).readlines() ] )
+    show()
diff --git a/eo/test/t-openmp.cpp b/eo/test/t-openmp.cpp
index 0e5f72388..4ccc54ad6 100644
--- a/eo/test/t-openmp.cpp
+++ b/eo/test/t-openmp.cpp
@@ -2,6 +2,10 @@
 // t-openmp.cpp
 //-----------------------------------------------------------------------------
 
+#include <fstream>
+#include <sstream>
+#include <climits>
+
 #include <utils/eoLogger.h>
 #include <utils/eoParserLogger.h>
 
@@ -24,13 +28,21 @@ typedef eoReal< eoMinimizingFitness > EOT;
 int main(int ac, char** av)
 {
     eoParserLogger parser(ac, av);
-    eoState state;
 
-    eoRealInitBounded<EOT>& init = make_genotype( parser, state, EOT() );
-    eoPop< EOT >& pop = make_pop( parser, state, init );
-    eoPop< EOT >& pop2 = make_pop( parser, state, init );
-    eoEvalFuncPtr< EOT, double, const std::vector< double >& > mainEval( real_value );
-    eoEvalFuncCounter<EOT> eval( mainEval );
+    unsigned int popMin = parser.getORcreateParam((unsigned int)1, "popMin", "Population Min", 'p', "Evolution Engine").value();
+    unsigned int popMax = parser.getORcreateParam((unsigned int)100, "popMax", "Population Max", 'P', "Evolution Engine").value();
+
+    unsigned int dimMin = parser.getORcreateParam((unsigned int)1, "dimMin", "Dimension Min", 'd', "Evolution Engine").value();
+    unsigned int dimMax = parser.getORcreateParam((unsigned int)100, "dimMax", "Dimension Max", 'D', "Evolution Engine").value();
+
+    unsigned int nRun = parser.getORcreateParam((unsigned int)100, "nRun", "Number of runs", 'r', "Evolution Engine").value();
+
+    std::string resultsFileName = parser.getORcreateParam(std::string("results.txt"), "resultsFileName", "Results file name", 0, "Evolution Engine").value();
+
+    double threshold = parser.getORcreateParam((double)3.0, "threshold", "Threshold of max speedup", 0, "Evolution Engine").value();
+
+    uint32_t seedParam = parser.getORcreateParam((uint32_t)0, "seed", "Random number seed", 0).value();
+    if (seedParam == 0) { seedParam = time(0); }
 
     if (parser.userNeedsHelp())
 	{
@@ -41,18 +53,63 @@ int main(int ac, char** av)
     make_help(parser);
     make_verbose(parser);
 
-    double ts1 = omp_get_wtime();
-    apply< EOT >( eval, pop );
-    //sleep(1);
-    double ts2 = omp_get_wtime();
+    rng.reseed( seedParam );
 
-    double tp1 = omp_get_wtime();
-    omp_apply< EOT >( eval, pop2 );
-    //sleep(1);
-    double tp2 = omp_get_wtime();
+    eoEvalFuncPtr< EOT, double, const std::vector< double >& > mainEval( real_value );
+    eoEvalFuncCounter<EOT> eval( mainEval );
 
-    eo::log << "Ts = " << ts2 - ts1 << std::endl;
-    eo::log << "Tp = " << tp2 - tp1 << std::endl;
+    eoUniformGenerator< double > gen(-5, 5);
+
+    std::ostringstream ss;
+    ss << resultsFileName << "-p" << popMin << "-P" << popMax << "-d" << dimMin << "-D" << dimMax << "-r" << nRun << "-t" << threshold << "-s" << seedParam;
+    std::ofstream resultsFile( ss.str().c_str() );
+
+    for ( size_t p = popMin; p <= popMax; ++p )
+	{
+	    for ( size_t d = dimMin; d <= dimMax; ++d )
+	    	{
+		    eo::log << eo::logging << p << 'x' << d << std::endl;
+
+		    for ( size_t r = 0; r < nRun; ++r )
+			{
+			    eoInitFixedLength< EOT > init( d, gen );
+
+			    double Ts;
+			    double Tp;
+
+			    // sequential scope
+			    {
+				eoPop< EOT > pop( p, init );
+				double t1 = omp_get_wtime();
+				apply< EOT >(eval, pop);
+				double t2 = omp_get_wtime();
+				Ts = t2 - t1;
+			    }
+
+			    // parallel scope
+			    {
+				eoPop< EOT > pop( p, init );
+				double t1 = omp_get_wtime();
+				omp_apply< EOT >(eval, pop);
+				double t2 = omp_get_wtime();
+				Tp = t2 - t1;
+			    }
+
+			    if ( ( Ts / Tp ) > threshold ) { continue; }
+
+			    resultsFile << Ts / Tp << ' ';
+
+			    eo::log << eo::debug;
+			    eo::log << "Ts = " << Ts << std::endl;
+			    eo::log << "Tp = " << Tp << std::endl;
+			    eo::log << "S_p = " << Ts / Tp << std::endl;
+			} // end of runs
+
+		    resultsFile << std::endl;
+
+		} // end of dimension
+
+	} // end of population
 
     return 0;
 }