1 files changed, 68 insertions, 0 deletions
diff --git a/tests/src/test_gensvm_simplex.c b/tests/src/test_gensvm_simplex.c
new file mode 100644
index 0000000..c13c4ea
--- /dev/null
+++ b/tests/src/test_gensvm_simplex.c
@@ -0,0 +1,68 @@
+/**
+ * @file test_gensvm_simplex.c
+ * @author Gertjan van den Burg
+ * @date May, 2016
+ * @brief Unit tests for gensvm_simplex.c functions
+ */
+
+#include "minunit.h"
+#include "gensvm_simplex.h"
+
+char *test_simplex_1()
+{
+	double *U = Calloc(double, 2*1);
+
+	gensvm_simplex(2, U);
+
+	mu_assert(matrix_get(U, 1, 0, 0) == -0.5, "U(0, 0) incorrect.");
+	mu_assert(matrix_get(U, 1, 1, 0) == 0.5, "U(1, 0) incorrect.");
+
+	free(U);
+
+	return NULL;
+}
+
+char *test_simplex_2()
+{
+	double *U = Calloc(double, 4*3);
+
+	gensvm_simplex(4, U);
+
+	mu_assert(matrix_get(U, 3, 0, 0) == -0.5, "U(0, 0) incorrect.");
+	mu_assert(matrix_get(U, 3, 1, 0) == 0.5, "U(1, 0) incorrect.");
+	mu_assert(matrix_get(U, 3, 2, 0) == 0.0, "U(2, 0) incorrect.");
+	mu_assert(matrix_get(U, 3, 3, 0) == 0.0, "U(3, 0) incorrect.");
+
+	mu_assert(fabs(matrix_get(U, 3, 0, 1) - -0.5/sqrt(3)) < 1e-14,
+		  	"U(0, 1) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 1, 1) - -0.5/sqrt(3)) < 1e-14,
+		  	"U(1, 1) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 2, 1) - 1.0/sqrt(3)) < 1e-14,
+		  	"U(2, 1) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 3, 1) - 0.0) < 1e-14,
+		  	"U(3, 1) incorrect.");
+
+	mu_assert(fabs(matrix_get(U, 3, 0, 2) - -1.0/sqrt(24)) < 1e-14,
+		  	"U(0, 2) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 1, 2) - -1.0/sqrt(24)) < 1e-14,
+		  	"U(1, 2) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 2, 2) - -1.0/sqrt(24)) < 1e-14,
+		  	"U(2, 2) incorrect.");
+	mu_assert(fabs(matrix_get(U, 3, 3, 2) - 3.0/sqrt(24)) < 1e-14,
+		       	"U(3, 2) incorrect.");
+
+	free(U);
+
+	return NULL;
+}
+
+char *all_tests()
+{
+	mu_suite_start();
+	mu_run_test(test_simplex_1);
+	mu_run_test(test_simplex_2);
+
+	return NULL;
+}
+
+RUN_TESTS(all_tests);