def run_once():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
solver = AnalysisBySynthesis(synth_solver)
pt = AnalysisPhaseTransition(signal_size,
dict_size,
deltas,
rhos,
num_data,
numpy.inf, [solver],
oper_type=operator_type)
filebasename = 'save/exact _' + partname
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6,
show=False,
basename=filebasename,
saveexts=['pdf', 'png'])
time_end = datetime.datetime.now()
print "End time: " + time_end.strftime("%Y-%m-%d --- %H:%M:%S:%f")
print "Elapsed: " + str((time_end - time_start).seconds) + " seconds"
print "------"
def run_analysis():
signal_size, dict_size = 50, 60
deltas = numpy.arange(0.1, 0.99, 0.1)
rhos = numpy.arange(0.1, 0.9, 0.1)
print "Running analysis phase transition..."
pt = AnalysisPhaseTransition(
signal_size, dict_size, deltas, rhos, 3, numpy.inf, [
AnalysisL1Min(1e-6),
AnalysisBySynthesis(L1Min(1e-6)),
AnalysisBySynthesis(
OrthogonalMatchingPursuit(1e-10, algorithm="sparsify_QR")),
AnalysisBySynthesis(SmoothedL0(1e-8)),
GreedyAnalysisPursuit(1e-6)
])
pt.run(processes=1)
pt.plot(thresh=1e-6)
def run_debias_type():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
debias_types = [False, True, 1e-6, 20, "real", "all"]
synth_solvers = [
IterativeHardThresholding(1e-20,
sparsity="real",
maxiter=100000,
debias=debias) for debias in debias_types
]
solvers = [
AnalysisBySynthesis(synth_solver,
nullspace_multiplier=1,
nullspace_multiplier_type="normalized_row")
for synth_solver in synth_solvers
]
pt = AnalysisPhaseTransition(signal_size,
dict_size,
deltas,
rhos,
num_data,
numpy.inf,
solvers,
oper_type=operator_type)
filebasename = 'save/exact_' + partname + '_debiastype'
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6,
subplot=True,
show=False,
basename=filebasename,
saveexts=['pdf', 'png'])
pt.plot_global_error(shape=(1, len(solvers)),
thresh=1e-6,
show=False,
basename=filebasename + '_globalerr',
saveexts=['pdf', 'png'])
time_end = datetime.datetime.now()
print "End time: " + time_end.strftime("%Y-%m-%d --- %H:%M:%S:%f")
print "Elapsed: " + str((time_end - time_start).seconds) + " seconds"
print "------"
def run_lambda_vals_normrow():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
lambdas = numpy.logspace(-4, 4, 9)
solvers = [
AnalysisBySynthesis(synth_solver,
nullspace_multiplier=lmbd,
nullspace_multiplier_type="normalized_row")
for lmbd in lambdas
]
pt = AnalysisPhaseTransition(signal_size,
dict_size,
deltas,
rhos,
num_data,
numpy.inf,
solvers,
oper_type=operator_type)
filebasename = 'save/exact_' + partname + '_lambdavals_normrow'
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6,
subplot=True,
show=False,
basename=filebasename,
saveexts=['pdf', 'png'])
pt.plot_global_error(shape=(1, len(solvers)),
thresh=1e-6,
show=False,
basename=filebasename + '_globalerr',
saveexts=['pdf', 'png'])
time_end = datetime.datetime.now()
print "End time: " + time_end.strftime("%Y-%m-%d --- %H:%M:%S:%f")
print "Elapsed: " + str((time_end - time_start).seconds) + " seconds"
print "------"