Add mtdev-matching test program

Measures matching speed and provides some historically difficult
test cases.

Signed-off-by: Henrik Rydberg <rydberg@euromail.se>

3 files changed, 146 insertions, 1 deletions

diff --git a/.gitignore b/.gitignore
index 132e096..fbe141a 100644
--- a/.gitignore
+++ b/.gitignore
@@ -77,4 +77,5 @@ core
 #
 test/mtdev-test
 test/mtdev-mapgen
+test/mtdev-matching
 patches
diff --git a/test/Makefile.am b/test/Makefile.am
index b1be019..12df79a 100644
--- a/test/Makefile.am
+++ b/test/Makefile.am
@@ -1,9 +1,12 @@
-noinst_PROGRAMS = mtdev-mapgen
+noinst_PROGRAMS = mtdev-mapgen mtdev-matching
 bin_PROGRAMS = mtdev-test
 
 INCLUDES=-I$(top_srcdir)/include/
 
 mtdev_mapgen_SOURCES = mtdev-mapgen.c
 
+mtdev_matching_SOURCES = mtdev-matching.c
+mtdev_matching_LDFLAGS = -L$(top_builddir)/src/.libs/ -lmtdev
+
 mtdev_test_SOURCES = mtdev-test.c
 mtdev_test_LDFLAGS = -lmtdev -L$(top_builddir)/src/.libs/
diff --git a/test/mtdev-matching.c b/test/mtdev-matching.c
new file mode 100644
index 0000000..dc68994
--- /dev/null
+++ b/test/mtdev-matching.c
@@ -0,0 +1,141 @@
+/*****************************************************************************
+ *
+ * mtdev - Multitouch Protocol Translation Library (MIT license)
+ *
+ * Copyright (C) 2010 Henrik Rydberg <rydberg@euromail.se>
+ * Copyright (C) 2010 Canonical Ltd.
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice (including the next
+ * paragraph) shall be included in all copies or substantial portions of the
+ * Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
+ * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ *
+ ****************************************************************************/
+
+#include <src/match.h>
+#include <stdio.h>
+#include <time.h>
+
+#define ITS 1000000
+
+static void test1()
+{
+	int A[] = {
+		1013,
+		3030660,
+		3559354,
+		12505925,
+		19008450,
+		6946421,
+		6118613,
+		698020,
+		3021800,
+		1017,
+		37573,
+		3242018,
+		8152794,
+		1266053,
+		942941,
+		462820,
+	};
+	int index[DIM_FINGER], i;
+	mtdev_match(index, A, 4, 4);
+	for (i = 0; i < 4; i++)
+		printf("match[%d] = %d\n", i, index[i]);
+}
+
+static void test2()
+{
+	int A[] = {
+		0,
+		4534330,
+		22653552,
+		12252500,
+		685352,
+		4534330,
+		0,
+		9619317,
+		28409530,
+		6710170,
+		22653552,
+		9619317,
+		0,
+		47015292,
+		29788572,
+		2809040,
+		10428866,
+		38615920,
+		17732500,
+		719528,
+		12113945,
+		28196220,
+		46778656,
+		405,
+		14175493,
+	};
+	int index[DIM_FINGER], i;
+	mtdev_match(index, A, 5, 5);
+	for (i = 0; i < 5; i++)
+		printf("match[%d] = %d\n", i, index[i]);
+}
+
+static void speed1()
+{
+	/* column-by-column matrix */
+	int A[DIM2_FINGER];
+	int x1[DIM_FINGER] = { 1, 5, 2, 3, 4, 5, 6, 7, 8 };
+	int y1[DIM_FINGER] = { 1, 5, 2, 3, 4, 6, 6, 7, 8 };
+	int x2[DIM_FINGER] = { 1.1, 3, 2, 4, 5, 6, 7, 8 };
+	int y2[DIM_FINGER] = { 1, 3, 2, 4, 5, 6, 7, 8 };
+	int index[DIM_FINGER];
+	int n1 = 4;
+	int n2 = 7;
+
+	int i, j;
+
+	for (i = 0; i < n1; i++) {
+		for (j = 0; j < n2; j++) {
+			A[i + n1 * j] =
+				(x1[i] - x2[j]) * (x1[i] - x2[j]) +
+				(y1[i] - y2[j]) * (y1[i] - y2[j]);
+		}
+	}
+
+	clock_t t1 = clock();
+	for (i = 0; i < ITS; i++)
+		mtdev_match(index, A, n1, n2);
+	clock_t t2 = clock();
+
+	printf("%lf matches per second\n",
+	       ITS * ((float)CLOCKS_PER_SEC / (t2 - t1)));
+
+	for (i = 0; i < n1; i++)
+		printf("match[%d] = %d\n", i, index[i]);
+
+}
+
+int main(int argc, char *argv[])
+{
+	printf("test1\n");
+	test1();
+	printf("test2\n");
+	test2();
+	printf("speed1\n");
+	speed1();
+	printf("done\n");
+	return 0;
+}