Merge Gael's experimental OpenMP parallelization support into Assign.h.

1 files changed, 18 insertions, 6 deletions

diff --git a/bench/benchmarkX.cpp b/bench/benchmarkX.cpp
index 09173e1ed..590f2636b 100644
--- a/bench/benchmarkX.cpp
+++ b/bench/benchmarkX.cpp
@@ -5,18 +5,30 @@
 using namespace std;
 USING_PART_OF_NAMESPACE_EIGEN
 
+#ifndef MATTYPE
+#define MATTYPE MatrixXd
+#endif
+
+#ifndef MATSIZE
+#define MATSIZE 20
+#endif
+
+#ifndef REPEAT
+#define REPEAT 100000
+#endif
+
 int main(int argc, char *argv[])
 {
-	MatrixXd I = MatrixXd::identity(20,20);
-	MatrixXd m(20,20);
-	for(int i = 0; i < 20; i++) for(int j = 0; j < 20; j++)
+	MATTYPE I = MATTYPE::identity(MATSIZE,MATSIZE);
+	MATTYPE m(MATSIZE,MATSIZE);
+	for(int i = 0; i < MATSIZE; i++) for(int j = 0; j < MATSIZE; j++)
 	{
-		m(i,j) = 0.1 * (i+20*j);
+		m(i,j) = 0.1 * (i+MATSIZE*j)/MATSIZE;
 	}
-	for(int a = 0; a < 100000; a++)
+	for(int a = 0; a < REPEAT; a++)
 	{
 		m = I + 0.00005 * (m + m*m);
 	}
-	cout << m << endl;
+	cout << m(0,0) << endl;
 	return 0;
 }