Context Navigation

test_normpdf.m @ 86

Last change on this file since 86 was 37, checked in by (none), 15 years ago
Added original make3d
File size: 853 bytes

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.