octave-vlfeat 0.9.17+dfsg0-6build1 / usr / share / octave / site / m / vlfeat / toolbox / xtest / vl_test_alldist2.m

This file is indexed.

/usr/share/octave/site/m/vlfeat/toolbox/xtest/vl_test_alldist2.m is in octave-vlfeat 0.9.17+dfsg0-6build1.

This file is owned by root:root, with mode 0o644.

The actual contents of the file can be viewed below. 

 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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
function results = vl_test_alldist2(varargin)
% VL_TEST_ALLDIST
vl_test_init ;

% TODO: test integer classes

function s = setup()
vl_twister('state', 0) ;
s.X = 3.1 * vl_twister(10,10) ;
s.Y = 4.7 * vl_twister(10,7) ;

function test_null_args(s)
vl_assert_equal(...
  vl_alldist2(zeros(15,12), zeros(15,0), 'kl2'), ...
  zeros(12,0)) ;

vl_assert_equal(...
  vl_alldist2(zeros(15,0), zeros(15,0), 'kl2'), ...
  zeros(0,0)) ;

vl_assert_equal(...
  vl_alldist2(zeros(15,0), zeros(15,12), 'kl2'), ...
  zeros(0,12)) ;

vl_assert_equal(...
  vl_alldist2(zeros(0,15), zeros(0,12), 'kl2'), ...
  zeros(15,12)) ;

function test_self(s)
vl_assert_almost_equal(...
  vl_alldist2(s.X, 'kl2'), ...
  makedist(@(x,y) x*y, s.X, s.X), ...
  1e-6) ;

function test_distances(s)
dists = {'chi2', 'l2', 'l1', 'hell', ...
         'kchi2', 'kl2', 'kl1', 'khell'} ;
distsEquiv = { ...
  @(x,y) (x-y)^2 / (x + y), ...
  @(x,y) (x-y)^2, ...
  @(x,y) abs(x-y), ...
  @(x,y) (sqrt(x) - sqrt(y))^2, ...
  @(x,y) 2 * (x*y) / (x + y), ...
  @(x,y) x*y, ...
  @(x,y) min(x,y), ...
  @(x,y) sqrt(x.*y)};
types = {'single', 'double', 'sparse'} ;

for simd = [0 1]
  for d = 1:length(dists)
    for t = 1:length(types)
      vl_simdctrl(simd) ;
      X = feval(str2func(types{t}), s.X) ;
      Y = feval(str2func(types{t}), s.Y) ;
      a = vl_alldist2(X,Y,dists{d}) ;
      b = makedist(distsEquiv{d},X,Y) ;
      vl_assert_almost_equal(a,b, ...
        1e-4, ...
        'alldist failed for dist=%s type=%s simd=%d', ...
        dists{d}, ...
        types{t}, ...
        simd) ;
    end
  end
end

function test_distance_kernel_pairs(s)
dists = {'chi2', 'l2', 'l1', 'hell'} ;
for d = 1:length(dists)
  dist = char(dists{d}) ;
  X = s.X ;
  Y = s.Y ;
  ker = ['k' dist] ;
  kxx = vl_alldist2(X,X,ker) ;
  kyy = vl_alldist2(Y,Y,ker) ;
  kxy = vl_alldist2(X,Y,ker) ;
  kxx = repmat(diag(kxx), 1, size(s.Y,2)) ;
  kyy = repmat(diag(kyy), 1, size(s.X,1))' ;
  d2  = vl_alldist2(X,Y,dist) ;
  vl_assert_almost_equal(d2, kxx + kyy - 2 * kxy, '1e-6') ;
end

function D = makedist(cmp,X,Y)
[d,m] = size(X) ;
[d,n] = size(Y) ;
D = zeros(m,n) ;
for i = 1:m
  for j = 1:n
    acc = 0 ;
    for k = 1:d
      acc = acc + cmp(X(k,i),Y(k,j)) ;
    end
    D(i,j) = acc ;
  end
end
conv = str2func(class(X)) ;
D = conv(D) ;

APT Browse - Built by Thomas Orozco.