Commit d671d684 authored by alexander.dornheim's avatar alexander.dornheim
Browse files

Minor changes to Siemens files

parent f27bc5a8
......@@ -9,4 +9,4 @@ Siemens = Phase(Siemens_nonneg_in.new_config)
Siemens.solve()
Siemens.compare_to_matlab()
#Siemens.compare_data_to_matlab()
#JWST.show()
Siemens.show()
......@@ -7,6 +7,6 @@ from phase import Phase
Siemens = Phase(Siemens_real_in.new_config)
Siemens.solve()
Siemens.compare_to_matlab()
#Siemens.compare_to_matlab()
#Siemens.compare_data_to_matlab()
Siemens.show()
......@@ -66,6 +66,7 @@ def Siemens_processor(config):
config['amplitude'] = config['amplitude']/norm(config['amplitude'],'fro')*config['norm_rt_data']
config['supp_phase'] = np.nonzero(np.ones((config['Ny'],config['Nx'])))
config['illumination_phase'] = np.nonzero(np.ones((config['Ny'],config['Nx'])))
elif config['object'] == 'complex':
# put some phase across S
points = config['Nx']
......@@ -86,6 +87,7 @@ def Siemens_processor(config):
config['data']=M**M
config['data_zeros'] = np.where(M==0)
config['support_idx'] = np.nonzero(S)
elif config['object'] == 'phase':
# put some phase across S
points = config['Nx']
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment