3 files changed, 10 insertions, 12 deletions
diff --git a/bench/benchmark.cpp b/bench/benchmark.cpp
index 99fd097ba..39ed317aa 100644
--- a/bench/benchmark.cpp
+++ b/bench/benchmark.cpp
@@ -18,12 +18,11 @@ USING_PART_OF_NAMESPACE_EIGEN
 
 int main(int argc, char *argv[])
 {
-    Matrix<SCALAR,MATSIZE,MATSIZE> I;
+    Matrix<SCALAR,MATSIZE,MATSIZE> I = Matrix<SCALAR,MATSIZE,MATSIZE>::ones();
     Matrix<SCALAR,MATSIZE,MATSIZE> m;
     for(int i = 0; i < MATSIZE; i++)
         for(int j = 0; j < MATSIZE; j++)
         {
-            I(i,j) = (i==j);
             m(i,j) = (i+MATSIZE*j);
         }
     asm("#begin");
diff --git a/bench/benchmarkX.cpp b/bench/benchmarkX.cpp
index 590f2636b..a269f695a 100644
--- a/bench/benchmarkX.cpp
+++ b/bench/benchmarkX.cpp
@@ -1,33 +1,32 @@
-// g++ -O3 -DNDEBUG benchmarkX.cpp -o benchmarkX && time ./benchmarkX
-
+// g++ -fopenmp -I .. -O3 -DNDEBUG -finline-limit=1000 benchmarkX.cpp -o b && time ./b
 #include <Eigen/Core>
 
 using namespace std;
 USING_PART_OF_NAMESPACE_EIGEN
 
 #ifndef MATTYPE
-#define MATTYPE MatrixXd
+#define MATTYPE MatrixXLd
 #endif
 
 #ifndef MATSIZE
-#define MATSIZE 20
+#define MATSIZE 400
 #endif
 
 #ifndef REPEAT
-#define REPEAT 100000
+#define REPEAT 100
 #endif
 
 int main(int argc, char *argv[])
 {
-	MATTYPE I = MATTYPE::identity(MATSIZE,MATSIZE);
+	MATTYPE I = MATTYPE::ones(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+MATSIZE*j)/MATSIZE;
+		m(i,j) = (i+j+1)/(MATSIZE*MATSIZE);
 	}
 	for(int a = 0; a < REPEAT; a++)
 	{
-		m = I + 0.00005 * (m + m*m);
+		m = I + 0.0001 * (m + m*m);
 	}
 	cout << m(0,0) << endl;
 	return 0;
diff --git a/bench/benchmarkXL.cpp b/bench/benchmarkXcwise.cpp
index ff128da06..dd29743cd 100644
--- a/bench/benchmarkXL.cpp
+++ b/bench/benchmarkXcwise.cpp
@@ -19,11 +19,11 @@ USING_PART_OF_NAMESPACE_EIGEN
 
 int main(int argc, char *argv[])
 {
-	MATTYPE I = MATTYPE::identity(MATSIZE,MATSIZE);
+	MATTYPE I = MATTYPE::ones(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+MATSIZE*j)/MATSIZE;
+		m(i,j) = 0.1 * (i+j+1)/(MATSIZE*MATSIZE);
 	}
 	for(int a = 0; a < REPEAT; a++)
 	{