1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
function test_kernel_post()
more off;
rand('state', 123456);
n_tr = 10;
n_tt = 8;
m = 5;
K = 3;
X_tr = rand(n_tr, m);
Z_tr = [ones(n_tr, 1), X_tr];
set_matrix(Z_tr, "train->RAW", "train->r+1");
fprintf("train->Z = train->RAW;\n");
X_tt = rand(n_tt, m);
Z_tt = [ones(n_tt, 1), X_tt];
set_matrix(Z_tt, "test->RAW", "test->m+1");
fprintf("test->Z = test->RAW;\n");
gamma = 1.132;
K2 = zeros(n_tt, n_tr);
for ii=1:n_tt
for jj=1:n_tr
K2(ii, jj) = exp(-gamma * sum((X_tt(ii, :) - X_tr(jj, :)).^2));
end
end
M = X_tr;
Sigma = rand(m, 1);
set_matrix(Sigma, "train->Sigma", "1");
N = K2 * M * diag(Sigma)^(-2);
test_Z = [ones(size(N, 1), 1), N];
assert_matrix(test_Z, "test->Z", "test->r+1");
end
function set_matrix(A, name, cols)
for ii=1:size(A, 1)
for jj=1:size(A, 2)
fprintf("matrix_set(%s, %s, %i, %i, %.16f);\n", name, cols, ii-1, jj-1, A(ii, jj));
end
end
fprintf("\n");
end
function assert_matrix(A, name, cols)
for ii=1:size(A, 1)
for jj=1:size(A, 2)
fprintf(["mu_assert(fabs(matrix_get(%s, %s, %i, %i) -\n%.16f) <", ...
" eps,\n\"Incorrect %s at %i, %i\");\n"], name, cols, ...
ii-1, jj-1, A(ii, jj), name, ii-1, jj-1);
end
end
fprintf("\n");
end
|
