diff options
Diffstat (limited to 'examples/testmnist.m')
-rw-r--r-- | examples/testmnist.m | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/examples/testmnist.m b/examples/testmnist.m new file mode 100644 index 0000000..38625a7 --- /dev/null +++ b/examples/testmnist.m | |||
@@ -0,0 +1,29 @@ | |||
1 | #! /usr/bin/octave -qf | ||
2 | % measuring Octave computing times | ||
3 | |||
4 | t0=time(); | ||
5 | load mnist.txt | ||
6 | disp("load"); | ||
7 | disp(time()-t0) | ||
8 | |||
9 | |||
10 | x = mnist(:,1:784); | ||
11 | d = mnist(:,785); | ||
12 | |||
13 | |||
14 | t0=time(); | ||
15 | xc = x - repmat(mean(x),rows(x),1); | ||
16 | disp("x - repmat(mean(x),rows(x),1)"); | ||
17 | disp(time()-t0) | ||
18 | |||
19 | t0=time(); | ||
20 | mc = (xc'*xc)/rows(x); | ||
21 | disp("(xc'*xc)/rows(x)"); | ||
22 | disp(time()-t0) | ||
23 | |||
24 | t0=time(); | ||
25 | [v,l]=eig(mc); | ||
26 | disp("eig"); | ||
27 | disp(time()-t0) | ||
28 | |||
29 | disp(flipud(diag(l))(1:10)); \ No newline at end of file | ||