Rev | Line | |
---|
[37] | 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.