legendfontsize = 24
cdfylabel = True #Label y axis as "CDF"
plot_dir = "plots"
plot_suffix = "SROM_pwlin_eol_CDF_m"
for m in sromsizes:
plot_suffix += "_" + str(m)
#Load / initialize target random variable from samples:
samples = np.genfromtxt(targetsamples)
target = SampleRandomVector(samples)
#Build up sromsize-to-SROM object map for plotting routine
sroms = OrderedDict()
for sromsize in sromsizes:
#Get EOL SROM Surrogate samples to make SampleRandomVector representation of CDF
eolsamplefile = "srom_eol_samples_m" + str(sromsize) + ".txt"
eolsamplefile = os.path.join(srom_dir, eolsamplefile)
eolsamples = np.genfromtxt(eolsamplefile)
sroms[sromsize] = SampleRandomVector(eolsamples)
Postprocessor.compare_srom_CDFs(sroms, target, plotdir="plots",
plotsuffix=plot_suffix, variablenames=varz, xlimits=xlimits, ylimits=ylimits, xticks=xticks,
cdfylabel=True, xaxispadding=xaxispadding,
axisfontsize=axisfontsize,
labelfontsize=labelfontsize,
legendfontsize=legendfontsize)
#Load / initialize target random variable from samples:
samples = np.genfromtxt(targetsamples)
target = SampleRandomVector(samples)
#Set x limits for each variable based on target:
xlimits = []
for i in range(target._dim):
lims = [np.min(samples[:, i]), np.max(samples[:, i])]
xlimits.append(lims)
#Build up sromsize-to-SROM object map for plotting routine
sroms = OrderedDict()
for sromsize in sromsizes:
#Generate SROM from file:
srom = SROM(sromsize, target._dim)
sromfile = "srom_m" + str(sromsize) + ".txt"
sromfile = os.path.join(srom_dir, sromfile)
srom.load_params(sromfile)
sroms[sromsize] = srom
Postprocessor.compare_srom_CDFs(sroms,
target,
plotdir="plots",
plotsuffix=plot_suffix,
variablenames=varz,
xlimits=xlimits,
xticks=xticks,
cdfylabel=cdfylabel)