0001 function [output] = generateRandomSample(model, n)
0002
0003 if (nargin < 1)
0004 disp 'function [output] = generateRandomSample(model, n)';
0005 return;
0006 end
0007 if (nargin < 2)
0008 n = 5000;
0009 end
0010
0011 m.A = model.S;
0012 m.lb = model.lb;
0013 m.ub = model.ub;
0014
0015 m = gpSampler(m,10,[],0,0);
0016 m.warmupPts = goodInitialPoint(model, n);
0017 mf = 1;
0018 while (mf > .52)
0019 [m,mf] = gpSampler(m,[],[],200,300);
0020
0021 end
0022 m.warmupPts = m.points;
0023 m = rmfield(m, 'points');
0024 mf = 1;
0025 while (mf > .52)
0026 [m,mf] = gpSampler(m,[],[],200,300);
0027
0028 end
0029 m.warmupPts = m.points;
0030 m = rmfield(m, 'points');
0031 mf = 1;
0032 while (mf > .52)
0033 [m,mf] = gpSampler(m,[],[],200,300);
0034
0035 end
0036
0037 output.point = m.points;
0038 output.mf = mf;
0039