if server: path_to_exec_file, num_mpi_threads, path_to_maxent = parameters.get_server_run( ) else: path_to_exec_file, num_mpi_threads, path_to_maxent = parameters.get_local_run( ) ############################################# # # # SPECIFICATION OF A MODEL # # # ############################################# parameters.set_model_parameters() lattice_type, beta, U, hartree_shift, Nk, num_of_neighbours, t, Coulomb, mu, particle_hole_symm, \ sweeps, time_limit, delta_mix, lambda_mix, max_it_num, start_it = parameters.get_model_parameters() # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # C T - H Y B S E G M E N T C A L C U L A T I O N # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # run CT-HYB SEGMENT solver number_of_fermionic_freqs = 1024 number_of_fermionic_freqs_for_fourier = 1000 # Because of noise we cut the tail of Delta (fermionic hybr. function) # off and make a Fouriet transform into the \tau - space by the first frequencies with smooth data. number_of_bosonic_frequencies = 1024 number_of_discrete_tau_points = 4096 # Friedrich - 4096 start_time = time.time()
KG.write(str(np.round(k_x, 6))) KG.write('\t') KG.write(str(np.round(k_y,6))) KG.write('\t') KG.write(str(np.round(interaction[i],6))) KG.write('\n') i += 1 GM.close() MK.close() KG.close() return interaction, kpoints_coordinates """ main part """ filename = 'Xw.dat' freq, X_w = iteration_cycle.read_real_function(filename) number_of_real_points = len(X_w) parameters.set_model_parameters() lattice_type, beta, U, hartree_shift, Nk, num_of_neighbours, t, V, mu, particle_hole_symm, sweeps, time_limit, delta_mix, lambda_mix, number_of_iterations, start_from_it = \ parameters.get_model_parameters() # 1. Selection G - M - K - G test, kpoints_coordinates = interaction_GMKG('test', Nk, lattice_type, V)
if server: path_to_exec_file, num_mpi_threads, path_to_maxent = parameters.get_server_run( ) else: path_to_exec_file, num_mpi_threads, path_to_maxent = parameters.get_local_run( ) ################################################# # SPECIFICATION OF THE MODEL # # (1) Copy the file set_parameters_change.py # # (2) Change the name as parameters.py # # (3) Set parameters # ################################################# parameters.set_model_parameters() lattice_type, beta, U, hartree_shift, Nk, num_of_neighbours, t, Coulomb, mu, particle_hole_symm, \ sweeps, time_limit, delta_mix, lambda_mix, number_of_max_it, start_it = parameters.get_model_parameters() # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # C T - H Y B S E G M E N T C A L C U L A T I O N # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # run CT-HYB SEGMENT solver number_of_fermionic_freqs = 1024 number_of_fermionic_freqs_for_fourier = 1024 # Because of noise we cut the tail of Delta (fermionic hybr. function) # off and make a Fouriet transform into the \tau - space by the first frequencies with smooth data. number_of_bosonic_frequencies = 1024 number_of_discrete_tau_points = 4096 # Friedrich - 4096 start_time = time.time()