home = vs.get_home() #%% LOAD VACUUM DATA vacuum_path = os.path.join(home, vacuum_folder) vacuum_file = lambda f, s: os.path.join(vacuum_path, f, s) vacuum_series = os.listdir(vacuum_path) vacuum_series = vu.filter_by_string_must(vacuum_series, "SC") vacuum_series = vu.sort_by_number(vacuum_series, 0) vacuum_data = [] vacuum_params = [] for s in vacuum_series: vacuum_data.append(np.loadtxt(vacuum_file(s, "Results.txt"))) vacuum_params.append(vs.retrieve_footer(vacuum_file(s, "Results.txt"))) vacuum_header = vs.retrieve_header(vacuum_file(s, "Results.txt")) vacuum_params = [vu.fix_params_dict(p) for p in vacuum_params] #%% LOAD WATER DATA water_path = os.path.join(home, water_folder) water_file = lambda f, s: os.path.join(water_path, f, s) water_series = os.listdir(water_path) water_series = vu.filter_by_string_must(water_series, "AllWater") water_series = vu.sort_by_number(water_series, 0) water_data = [] water_params = []
series_mustnt): path.append(os.path.join(home, f)) file.append(lambda f, s: os.path.join(path[-1], f, s)) series.append(os.listdir(path[-1])) series[-1] = vu.filter_to_only_directories(series[-1]) series[-1] = vu.filter_by_string_must(series[-1], sm) if smn != "": series[-1] = vu.filter_by_string_must(series[-1], smn, False) series[-1] = sf(series[-1]) data.append([]) params.append([]) for s in series[-1]: data[-1].append(np.loadtxt(file[-1](s, "Results.txt"))) params[-1].append(vs.retrieve_footer(file[-1](s, "Results.txt"))) header.append(vs.retrieve_header(file[-1](s, "Results.txt"))) for i in range(len(params[-1])): if not isinstance(params[-1][i], dict): params[-1][i] = vu.fix_params_dict(params[-1][i]) r = [] from_um_factor = [] resolution = [] paper = [] index = [] for p in params: r.append([pi["r"] for pi in p]) from_um_factor.append([pi["from_um_factor"] for pi in p]) resolution.append([pi["resolution"] for pi in p])
folder = "AuMieSphere/AuMie/8)TestParNP/2nd" home = vs.get_home() nprocesses = np.arange(8) + 1 #%% LOAD DATA path = os.path.join(home, folder) file = lambda f, s: os.path.join(path, f, s) series = os.listdir(path) series.sort() params = [] for s in series: params.append(vs.retrieve_footer(file(s, "Results.txt"))) params = [vu.fix_params_dict(p) for p in params] enlapsed = [p["enlapsed"] for p in params] total = np.array([sum(e) for e in enlapsed]) enlapsed = np.array(enlapsed) #%% PLOT plt.figure() # plt.plot(nprocesses, total_0, 'ro-', label="Serie 1") plt.plot(nprocesses, total, 'bo-', label="Serie 2") plt.xlabel("Número de subprocesos") plt.ylabel("Tiempo total (s)") plt.legend(["Serie 1", "Serie 2"])
""" import v_meep as vm import v_save as vs import numpy as np import os sysname = vs.get_sys_name() syshome = vs.get_sys_home() home = vs.get_home() #%% flux_file = os.path.join(home, "FluxData/FluxDataDirectory.txt") flux_dir = vs.retrieve_footer(flux_file) #%% data_files = flux_dir["path"] found_files = [] found_params = [] for df in data_files: try: found_params.append(vs.retrieve_footer(os.path.join(df, "Results.txt"))) found_files.append(df) except: print("No file found at ", df) print(f"Only found {100*len(found_files)/len(data_files):.2f}% of the files")
#%% PARAMETERS # Saving directories series = ["2020111101", "2020111103"] folder = "AuMieResults" home = "/home/vall/Documents/Thesis/ThesisPython/" #%% LOAD DATA file = lambda f, s: os.path.join(home, folder, "{}Results".format(s), f) data = [] params = [] for s in series: data.append(np.loadtxt(file("MidFlux.txt", s))) params.append(vs.retrieve_footer(file("MidFlux.txt", s))) header = vs.retrieve_header(file("MidFlux.txt", s)) r = [p["r"] for p in params] resolution = [p["resolution"] for p in params] from_um_factor = [p["from_um_factor"] for p in params] until_after_sources = [p["until_after_sources"] for p in params] #%% PLOT new_series = ["F8"] #["", "F1", "F2"] def proposed_factor(series, r, resolution, from_um_factor, until_after_sources): if series == "":
series = [ "2020102601", "2020102708", "2020102709", *["202010271" + str(i) for i in range(3)] ] folder = "MieResults" home = "/home/vall/Documents/Thesis/ThesisPython/MeepTutorial/" #%% params = [] for s in series: path = os.path.join(home, folder, "{}Results".format(s)) file = lambda f: os.path.join(path, f) p = vs.retrieve_footer(file("MidFlux.txt")) f = h5.File(file("MidField.h5"), "r+") for a in p.keys(): f["Ez"].attrs[a] = p[a] f.close() params.append(p) #%% resolution = [p["resolution"] for p in params] pml_width = [p["pml_width"] for p in params] air_width = [p["air_width"] for p in params] cutoff = [p["cutoff"] for p in params] until_after_sources = [p["until_after_sources"] for p in params]