Skip to content
Snippets Groups Projects
Commit f14282db authored by Marco Ratto's avatar Marco Ratto
Browse files

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

(cherry picked from commit 3307b4b7)
parent 9e4670f1
Branches
Tags
No related merge requests found
......@@ -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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment