/usr/share/octave/site/m/vlfeat/toolbox/demo/vl_demo_kmeans_init.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 | % VL_DEMO_KMEANS_INIT
if 0
I = vl_imtest('box') ;
[f, d] = vl_dsift(single(I) / 255, 'step', 1, 'size', 3) ;
X = double(d) ;
numCenters = 300 ;
else
X = [] ;
numCenters = 40 ;
for i=1:numCenters
X_ = randn(10,100) ;
X_(i) = X_(i) + 100 ;
X = [X X_] ;
end
numCenters = numCenters * 2 ;
end
numData = size(X, 2) ;
dimension = size(X, 1) ;
clear elaps energy ;
dimension = 128 ;
numData = 1000 ;
numCenters = 100 ;
numTrials = 5 ;
numIters = 100 ;
vl_twister(0) ;
for trial = 1:numTrials
tic ;
[C, idx, e] = vl_kmeans(X, numCenters, 'Verbose', ...
'Initialization', 'randsel', ...
'MaxNumIterations', numIters) ;
elaps.vlrs(trial) = toc ;
energy.vlrs(trial) = e ;
tic ;
[C_, idx_, e_] = vl_kmeans(X, numCenters, 'Verbose', ...
'Initialization', 'plusplus', ...
'MaxNumIterations', numIters) ;
elaps.vlpp(trial) = toc ;
energy.vlpp(trial) = e_ ;
end
figure(1) ; clf ;
subplot(1,2,1) ; title('Energy') ;
mu = [mean(elaps.vlrs) mean(elaps.vlpp)] ;
st = [std(elaps.vlrs) std(elaps.vlpp)] ;
bar(mu) ; hold on ;
errorbar(mu, st, 'linestyle', 'none', 'color', 'r', 'linewidth', 4) ;
title('Time') ;
subplot(1,2,2) ;
mu = [mean(energy.vlrs) mean(energy.vlpp)] ;
st = [std(energy.vlrs) std(energy.vlpp)] ;
bar(mu) ; hold on ;
errorbar(mu, st, 'linestyle', 'none', 'color', 'r', 'linewidth', 4) ;
title('Energy') ;
|