diff --git a/test/tests.cc b/test/tests.cc
new file mode 100644
index 0000000..ed070d3
--- /dev/null
+++ b/test/tests.cc
@@ -0,0 +1,8 @@
+
+#include <iostream>
+#include <cassert>
+using namespace std;
+
+int main()
+{
+}
diff --git a/util/Matrix.cc b/util/Matrix.cc
index 5d53153..a1e7d58 100644
--- a/util/Matrix.cc
+++ b/util/Matrix.cc
@@ -1,5 +1,8 @@
 
 #include "Matrix.h"
+#include <math.h>       /* fabs() */
+
+#define FP_EQ(x,y) (fabs((x)-(y)) < 0.00001)
 
 Matrix::Matrix()
 {
@@ -193,3 +196,7 @@ Vector operator*(const Matrix & m, const Vector & v)
     }
     return res;
 }
+
+bool operator==(const Matrix & m1, const Matrix & m2)
+{
+}
diff --git a/util/Matrix.h b/util/Matrix.h
index 7023a70..3c0bf80 100644
--- a/util/Matrix.h
+++ b/util/Matrix.h
@@ -24,6 +24,7 @@ class Matrix
 
 Matrix operator*(const Matrix & m1, const Matrix & m2);
 Vector operator*(const Matrix & m, const Vector & v);
+bool operator==(const Matrix & m1, const Matrix & m2);
 
 #endif