def run_once():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
filebasename = 'save/exact_' + partname + '_' + str(lambda1) + '_' + str(
lambda2) + '_' + lambda2_type
pt = AnalysisPhaseTransition(signal_size,
dict_size,
deltas,
rhos,
num_data,
numpy.inf, [solver],
oper_type=operator_type)
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6,
show=False,
basename=filebasename,
saveexts=['pdf', 'png'])
pt.compute_global_average_error(shape=[1],
thresh=1e-6,
textfilename=filebasename + '.txt')
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_type():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
types = ["value", "normalized_row"]
solvers = [
AnalysisBySynthesis(synth_solver, nullspace_multiplier=1, nullspace_multiplier_type=type) for type in types
]
pt = AnalysisPhaseTransition(
signal_size, dict_size, deltas, rhos, num_data, numpy.inf, solvers, oper_type=operator_type
)
filebasename = "save/exact_" + partname + "_lambdatype"
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_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 = [ApproximateMessagePassing(stoptol=1e-14, maxiter=10000, 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_uap():
signal_size, dict_size = 50, 70
deltas = numpy.arange(0.1, 1, 0.1)
rhos = numpy.arange(0.1, 1, 0.1)
print "Running analysis phase transition..."
pt = AnalysisPhaseTransition(signal_size, dict_size, deltas, rhos, 3,
[UnconstrainedAnalysisPursuit(1e-6, 1, 1),
GreedyAnalysisPursuit(1e-6)])
pt.run()
pt.plot()
print "Example finished."
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_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_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_once():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
filebasename = "save/exact _" + partname
solver = AnalysisBySynthesis(synth_solver, nullspace_multiplier_type="normalized_row")
pt = AnalysisPhaseTransition(
signal_size, dict_size, deltas, rhos, num_data, numpy.inf, [solver], oper_type=operator_type
)
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6, show=False, basename=filebasename, saveexts=["pdf", "png"])
# filebasename_synth = filebasename + '_synth'
# pt_synth = SynthesisPhaseTransition(signal_size, dict_size, deltas, rhos, num_data, numpy.inf, [synth_solver])
# pt_synth.run()
# pt_synth.savedata(filebasename_synth)
# pt_synth.plot(thresh=1e-6, show=False, basename=filebasename_synth, 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 "------"
def run_once():
time_start = datetime.datetime.now()
print "Start time: " + time_start.strftime("%Y-%m-%d --- %H:%M:%S:%f")
filebasename = 'save/exact_'+ partname + '_' + str(lambda1) + '_' + str(lambda2) + '_' + lambda2_type
pt = AnalysisPhaseTransition(signal_size, dict_size, deltas, rhos, num_data, numpy.inf, [solver], oper_type=operator_type)
pt.run()
pt.savedata(filebasename)
pt.plot(thresh=1e-6, show=False, basename=filebasename, saveexts=['pdf', 'png'])
pt.compute_global_average_error(shape=[1], thresh=1e-6, textfilename=filebasename+'.txt')
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():
signal_size, dict_size = 200, 240
deltas = numpy.arange(0.1, 1, 0.3)
rhos = numpy.arange(0.1, 1, 0.3)
print "Running analysis phase transition..."
pt = AnalysisPhaseTransition(
signal_size,
dict_size,
deltas,
rhos,
3,
[
GreedyAnalysisPursuit(1e-8),
# AnalysisL1Min(1e-8),
# AnalysisBySynthesis(L1Min(1e-6)),
# AnalysisBySynthesis(OrthogonalMatchingPursuit(1e-9, algorithm="sparsify_QR")),
# AnalysisBySynthesis(SmoothedL0(1e-8)),
# AnalysisBySynthesis(TwoStageThresholding(1e-8, maxiter=3000)),
# AnalysisBySynthesis(IterativeHardThresholding(1e-10, sparsity="real", maxiter=1000))
],
)
pt.run()
def run_lambda_vals():
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) 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'
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_uap():
signal_size, dict_size = 50, 70
deltas = numpy.arange(0.1, 1, 0.1)
rhos = numpy.arange(0.1, 1, 0.1)
print "Running analysis phase transition..."
pt = AnalysisPhaseTransition(signal_size, dict_size, deltas, rhos, 3, [
UnconstrainedAnalysisPursuit(1e-6, 1, 1),
GreedyAnalysisPursuit(1e-6)
])
pt.run()
pt.plot()
print "Example finished."
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)
