PySDM_examples.Bartman_2020_MasterThesis.fig_4_adaptive_sdm

View Source

 1import os
 2
 3from matplotlib import pyplot as plt
 4from PySDM_examples.Shima_et_al_2009.example import run
 5from PySDM_examples.Shima_et_al_2009.settings import Settings
 6from PySDM_examples.Shima_et_al_2009.spectrum_plotter import SpectrumPlotter
 7
 8
 9def main(plot: bool = True, save: str = None):
10    n_sds = [13, 15, 17]
11    dts = [10, 5, 1, "adaptive"]
12    iters = 10
13    base_time = None
14
15    plt.ioff()
16    fig, axs = plt.subplots(
17        len(dts), len(n_sds), sharex=True, sharey=True, figsize=(10, 10)
18    )
19
20    for i, dt in enumerate(dts):
21        for j, n_sd in enumerate(n_sds):
22            outputs = []
23            exec_time = 0
24            for _ in range(iters):
25                settings = Settings()
26
27                settings.n_sd = 2**n_sd
28                settings.dt = dt if dt != "adaptive" else 10
29                settings.adaptive = dt == "adaptive"
30
31                states, exec_time = run(settings)
32                outputs.append(states)
33            mean_time = exec_time / iters
34            if base_time is None:
35                base_time = mean_time
36            norm_time = mean_time / base_time
37            mean_output = {}
38            for key in outputs[0].keys():
39                mean_output[key] = sum((output[key] for output in outputs)) / len(
40                    outputs
41                )
42
43            plotter = SpectrumPlotter(settings, legend=False)
44            plotter.fig = fig
45            plotter.ax = axs[i, j]
46            plotter.smooth = True
47            for step, vals in mean_output.items():
48                plotter.plot(vals, step * settings.dt)
49
50            plotter.ylabel = (
51                r"$\bf{dt: " + str(dt) + "}$\ndm/dlnr [g/m^3/(unit dr/r)]"
52                if j == 0
53                else None
54            )
55            plotter.xlabel = (
56                "particle radius [µm]\n" + r"$\bf{n_{sd}: 2^{" + str(n_sd) + "}}$"
57                if i == len(dts) - 1
58                else None
59            )
60            plotter.title = f"norm. time: {norm_time:.2f}; " + plotter.title
61            plotter.finished = False
62            plotter.finish()
63    if save is not None:
64        n_sd = settings.n_sd
65        plotter.save(save + "/" + f"{n_sd}_shima_fig_2" + "." + plotter.format)
66    if plot:
67        plotter.show()
68
69
70if __name__ == "__main__":
71    main(plot="CI" not in os.environ, save=".")

def main(plot: bool = True, save: str = None): View Source

10def main(plot: bool = True, save: str = None):
11    n_sds = [13, 15, 17]
12    dts = [10, 5, 1, "adaptive"]
13    iters = 10
14    base_time = None
15
16    plt.ioff()
17    fig, axs = plt.subplots(
18        len(dts), len(n_sds), sharex=True, sharey=True, figsize=(10, 10)
19    )
20
21    for i, dt in enumerate(dts):
22        for j, n_sd in enumerate(n_sds):
23            outputs = []
24            exec_time = 0
25            for _ in range(iters):
26                settings = Settings()
27
28                settings.n_sd = 2**n_sd
29                settings.dt = dt if dt != "adaptive" else 10
30                settings.adaptive = dt == "adaptive"
31
32                states, exec_time = run(settings)
33                outputs.append(states)
34            mean_time = exec_time / iters
35            if base_time is None:
36                base_time = mean_time
37            norm_time = mean_time / base_time
38            mean_output = {}
39            for key in outputs[0].keys():
40                mean_output[key] = sum((output[key] for output in outputs)) / len(
41                    outputs
42                )
43
44            plotter = SpectrumPlotter(settings, legend=False)
45            plotter.fig = fig
46            plotter.ax = axs[i, j]
47            plotter.smooth = True
48            for step, vals in mean_output.items():
49                plotter.plot(vals, step * settings.dt)
50
51            plotter.ylabel = (
52                r"$\bf{dt: " + str(dt) + "}$\ndm/dlnr [g/m^3/(unit dr/r)]"
53                if j == 0
54                else None
55            )
56            plotter.xlabel = (
57                "particle radius [µm]\n" + r"$\bf{n_{sd}: 2^{" + str(n_sd) + "}}$"
58                if i == len(dts) - 1
59                else None
60            )
61            plotter.title = f"norm. time: {norm_time:.2f}; " + plotter.title
62            plotter.finished = False
63            plotter.finish()
64    if save is not None:
65        n_sd = settings.n_sd
66        plotter.save(save + "/" + f"{n_sd}_shima_fig_2" + "." + plotter.format)
67    if plot:
68        plotter.show()