Line | |
---|
1 | % timing test for normpdf
|
---|
2 |
|
---|
3 | niter = 1000;
|
---|
4 | d = 100;
|
---|
5 | A = rand(d);
|
---|
6 | A = A*A';
|
---|
7 | S = chol(A);
|
---|
8 | iA = inv(A);
|
---|
9 | iS = chol(iA);
|
---|
10 | x = rand(d,1);
|
---|
11 |
|
---|
12 | fprintf('time for standard normal:')
|
---|
13 | tic;for i = 1:niter normpdfln(x); end;toc
|
---|
14 | fprintf('time for single mean:')
|
---|
15 | tic;for i = 1:niter normpdfln(x,x(:,1)); end;toc
|
---|
16 | fprintf('time for zero mean, Cholesky argument:')
|
---|
17 | tic;for i = 1:niter normpdfln(x,[],S); end;toc
|
---|
18 | fprintf('time for Cholesky argument:')
|
---|
19 | tic;for i = 1:niter normpdfln(x,x,S); end;toc
|
---|
20 | fprintf('time for variance argument:')
|
---|
21 | tic;for i = 1:niter normpdfln(x,x,[],A); end;toc
|
---|
22 | fprintf('time for inverse Cholesky argument:')
|
---|
23 | tic;for i = 1:niter normpdfln(x,x,iS,'inv'); end;toc
|
---|
24 | fprintf('time for inverse variance argument:')
|
---|
25 | tic;for i = 1:niter normpdfln(x,x,'inv',iA); end;toc
|
---|
26 |
|
---|
27 | x=[0;1];
|
---|
28 | m=[0.1;1.1];
|
---|
29 | v=[0.1 0.1;0.1 0.5];
|
---|
30 | normpdf(x,m,[],v)
|
---|
Note: See
TracBrowser
for help on using the repository browser.