From 3307b4b70e0abd029e00f5bd1face00c2fe3c3d3 Mon Sep 17 00:00:00 2001 From: Marco Ratto Date: Fri, 29 Jun 2012 00:36:02 +0200 Subject: [PATCH] Use Sobol sequence with seed = 1 (to avoid [0 0 0 0 0 ... ] point ) --- matlab/gsa/stab_map_.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/matlab/gsa/stab_map_.m b/matlab/gsa/stab_map_.m index 049c659bb..18bb43755 100644 --- a/matlab/gsa/stab_map_.m +++ b/matlab/gsa/stab_map_.m @@ -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);