1 files changed, 59 insertions, 0 deletions
diff --git a/tests/aux/test_kernel_compute.m b/tests/aux/test_kernel_compute.m
new file mode 100644
index 0000000..7cc9388
--- /dev/null
+++ b/tests/aux/test_kernel_compute.m
@@ -0,0 +1,59 @@
+
+function kernel_tests(kerneltype)
+
+rand('state', 123456);
+more off;
+
+n = 10;
+m = 3;
+
+X = rand(n, m);
+Z = [ones(n, 1), X];
+
+for ii=1:n
+  for jj=1:m+1
+    fprintf("matrix_set(data->RAW, data->m+1, %i, %i, %.16f);\n", ii-1, jj-1, Z(ii, jj));
+  end
+end
+
+
+K = zeros(n, n);
+if strcmp(kerneltype, 'poly')
+  # Polynomial kernel
+  # (gamma * <x_1, x_2> + c)^d
+  gamma = 1.5;
+  c = 3.0;
+  d = 1.78;
+  
+  for ii=1:n
+    for jj=1:n
+      K(ii, jj) = (gamma * (X(ii, :) * X(jj, :)') + c)^d;
+    end
+  end
+elseif strcmp(kerneltype, 'rbf')
+  # RBF kernel
+  # exp(-gamma * norm(x1 - x2)^2)
+  gamma = 0.348
+  for ii=1:n
+    for jj=1:n
+      K(ii, jj) = exp(-gamma * sum((X(ii, :) - X(jj, :)).^2));
+    end
+  end
+elseif strcmp(kerneltype, 'sigmoid')
+  # Sigmoid kernel
+  # tanh(gamma * <x_1, x_2> + c)
+  gamma = 1.23;
+  c = 1.6;
+  for ii=1:n
+    for jj=1:n
+      K(ii, jj) = tanh(gamma * (X(ii, :) * X(jj, :)') + c);
+    end
+  end
+end
+
+fprintf('\n');
+for ii=1:n
+  for jj=1:n
+    fprintf("mu_assert(fabs(matrix_get(K, data->n, %i, %i) -\n %.16f) < eps,\n\"Incorrect K at %i, %i\");\n", ii-1, jj-1, K(ii, jj), ii-1, jj-1);
+  end
+end
+\ No newline at end of file