Use Sobol sequence with seed = 1 (to avoid [0 0 0 0 0 ... ] point )

time-shift
Marco Ratto 2012-06-29 00:36:02 +02:00
parent aa62568a44
commit 3307b4b70e
1 changed files with 1 additions and 1 deletions

View File

@ -127,7 +127,7 @@ if fload==0,
% Nsam=size(lpmat,1);
else
if np<52 && ilptau>0,
[lpmat] = qmc_sequence(np, int64(0), 0, Nsam)';
[lpmat] = qmc_sequence(np, int64(1), 0, Nsam)';
if np>30 || ilptau==2, % scrambled lptau
for j=1:np,
lpmat(:,j)=lpmat(randperm(Nsam),j);