def main():
loaddir = "../../../data/"
probdir = "prob_trans.inputs.sigma2_u1.100.0.sigma2_u2.1.0.sigma2_v.1.0.m.20.nofu1s.5.N.25000"
utils.setsizes()
params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt")
ymax = data[params==0.001]
plot(params, data/ymax, label="LMS")
#params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt")
#plot(params, data/ymax, label="K-Y")
params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt")
plot(params, data/ymax, label="IDBD")
params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt")
plot(params, data/ymax, label="Autostep")
params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt")
plot(params, data/ymax, label="RLS")
legend(loc='lower left')
xlim([10**-4, 1])
ylim([0, 1.2])
xscale("log")
xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
ylabel("$S$ (MSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
def main():
loaddir = "../../../data/"
probdir = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.1.0.m.20.o1s.5.N.50000"
utils.setsizes()
params, data = dataload(dir=loaddir+"lms/"+probdir, paramfile="stepsizes.txt")
ymax = 1
plot(params, data/ymax, label="LMS")
# params, data = dataload(dir=loaddir+"bmp/"+probdir, paramfile="metastepsizes.txt")
# plot(params, data/ymax, label="BMP")
params, data = dataload(dir=loaddir+"idbd/"+probdir, paramfile="metastepsizes.txt")
plot(params, data/ymax, label="IDBD")
params, data = dataload(dir=loaddir+"autostep/"+probdir, paramfile="metastepsizes.txt")
plot(params, data/ymax, label="Autostep")
params, data = dataload(dir=loaddir+"rls/"+probdir, paramfile="one_minus_gammas.txt")
plot(params, data/ymax, label="RLS")
legend(loc='lower left')
xlim([10**-11, 1])
ylim([0, 140])
xscale("log")
xlabel("tunable parameter", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
ylabel("RMSE", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())
def main():
utils.setsizes()
mpl.rcParams['lines.markeredgewidth'] = utils.geterrorbarlinewidth()
loaddir = "../../../data/"
probdir1 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.0.1.s2o1.0.0.m.20.o1s.5.N.50000"
probdir2 = "prob.s2u.1.0.a.0.9998.s2v.1.0.s2o1.100.0.s2o1.0.0.m.20.o1s.5.N.50000"
params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt")
ymax1 = data1[params==0.05]
ymax2 = data2[params==0.05]
#plot(params, (data1/ymax1 + data2/ymax2)/2.0)
params, data1, data2 = jointdataload(dir1=loaddir+"lms/"+probdir1, dir2=loaddir+"lms/"+probdir2, paramfile="stepsizes.txt")
plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="LMS")
params, data1, data2 = jointdataload(dir1=loaddir+"idbd/"+probdir1, dir2=loaddir+"idbd/"+probdir2, paramfile="metastepsizes.txt")
plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="IDBD")
params, data1, data2 = jointdataload(dir1=loaddir+"autostep/"+probdir1, dir2=loaddir+"autostep/"+probdir2, paramfile="metastepsizes.txt")
plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="Autostep")
#params, data1, data2 = jointdataload(dir1=loaddir+"bmp/"+probdir1, dir2=loaddir+"bmp/"+probdir2, paramfile="metastepsizes.txt")
#plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="K-Y")
params, data1, data2 = jointdataload(dir1=loaddir+"rls/"+probdir1, dir2=loaddir+"rls/"+probdir2, paramfile="one_minus_gammas.txt")
plot(params, (data1/ymax1 + data2/ymax2)/2.0, label="RLS")
legend(loc='lower left')
xlim([10**-9, 1])
ylim([0.6, 1.2])
xscale("log")
xlabel("tunable parameter", fontsize=utils.getxlabelsize())
ylabel("$S$ (RMSE relative to best LMS)", fontsize=utils.getxlabelsize(), labelpad=utils.getxlabelpadsize())