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
|
function [V] = test_optimize()
clear;
more off;
rand('state', 902183);
n = 8;
m = 3;
K = 4;
X = rand(n, m);
set_matrix(X);
y = [2 1 3 2 3 2 4 1];
set_vector(y);
p = 1.2143;
kappa = 0.90298;
lambda = 0.00219038;
epsilon = 1e-15;
rho = rand(n, 1);
set_vector(rho);
[W, t] = msvmmaj(X, y, rho, p, kappa, lambda, epsilon);
V = [t'; W];
end
function set_matrix(A)
for i=1:size(A, 1)
for j=1:size(A, 2)
fprintf('matrix_set(A, %i, %i, %i, %.16f);\n', size(A, 2), i-1, j-1, A(i, j));
end
end
end
function set_vector(a)
for i=1:numel(a)
fprintf('A[%i] = %.16f;\n', i-1, a(i));
end
end
|
