Exemplo n.º 1
0
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")
Exemplo n.º 5
0
#%% 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 == "":
Exemplo n.º 6
0
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]