stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added unitary tests for qmc_sequence mex file.

time-shift
Stéphane Adjemian (Charybdis) 2011-12-22 18:16:24 +01:00
parent 199eb98aeb
commit 627b24ddc6
1 changed files with 50 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

51
matlab/qmc_sequence.m
View File

@ -81,4 +81,53 @@
%$ end
%$
%$ T = all(t);
%@eof:1
%@eof:1
%@test:2
%$ t = ones(3,1);
%$
%$ d = 2;
%$ n = 100;
%$ s = int64(0);
%$
%$ [draws1, S] = qmc_sequence(d,s,0,n);
%$ [draws2, Q] = qmc_sequence(d,S,0,n);
%$ [draws3, P] = qmc_sequence(d,s,0,2*n);
%$
%$ t(1) = dyn_assert(s,int64(0));
%$ t(2) = dyn_assert(P,Q);
%$ t(3) = dyn_assert([draws1,draws2],draws3);
%$ T = all(t);
%@eof:2
%@test:3
%$ t = ones(3,1);
%$
%$ d = 2;
%$ n = 100;
%$ s = int64(0);
%$
%$ [draws1, S] = qmc_sequence(d,s,0,n,[0 , 2; -1, 2]);
%$ [draws2, Q] = qmc_sequence(d,s,0,n);
%$
%$ draws3 = draws2;
%$ draws3(1,:) = 2*draws2(1,:);
%$ draws3(2,:) = 3*draws2(2,:)-1;
%$ t(1) = dyn_assert(S,Q);
%$ t(2) = dyn_assert(draws1,draws3);
%$ T = all(t);
%@eof:3
%@test:4
%$ t = ones(3,1);
%$
%$ d = 2;
%$ n = 100;
%$ s = int64(0);
%$
%$ [draws, S] = qmc_sequence(d,s,2,n,.1);
%$ draws
%$ t(1) = 1;
%$ t(2) = 1;
%$ T = all(t);
%@eof:4