dedebenui/dispersionapp
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Better handling of no ZDW found

Browse Source
This commit is contained in:
Benoît Sierro
2023-03-21 10:52:21 +01:00
parent 9658839dc7
commit 6da66b5990
1 changed files with 16 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/dispersionapp/gui.py
View File

@@ -23,12 +23,12 @@ def app(config_file: os.PathLike | None = None):
with PlotApp( with PlotApp(
f"Dispersion design with {config.gas.title()}", f"Dispersion design with {config.gas.title()}",
core_diameter_um=np.linspace(50, 300), core_diameter_um=np.arange(50, 301, dtype=float),
pressure_mbar=np.geomspace(1, 2000), pressure_mbar=np.geomspace(1, 2000),
wall_thickness_um=np.geomspace(0.01, 10), wall_thickness_um=np.geomspace(0.01, 10),
n_tubes=np.arange(6, 16), n_tubes=np.arange(6, 16),
gap_um=np.linspace(1, 15), gap_um=np.arange(1, 15.5, 0.5),
t_fwhm_fs=np.linspace(10, 200, 96), t_fwhm_fs=np.arange(10, 201, dtype=float),
) as app: ) as app:
# initial setup # initial setup
app[0].horizontal_line("reference", 0, color="gray") app[0].horizontal_line("reference", 0, color="gray")
@@ -64,7 +64,7 @@ def app(config_file: os.PathLike | None = None):
disp = compute_capillary(core_diameter_um, pressure_mbar) disp = compute_capillary(core_diameter_um, pressure_mbar)
wl_zero_disp = sc.math.all_zeros(wl, disp) wl_zero_disp = sc.math.all_zeros(wl, disp)
if len(wl_zero_disp) == 0: if len(wl_zero_disp) == 0:
return raise ValueError("No zero dispersion found")
wl_zero_disp = wl_zero_disp[0] wl_zero_disp = wl_zero_disp[0]
with warnings.catch_warnings(): with warnings.catch_warnings():
@@ -138,15 +138,22 @@ def app(config_file: os.PathLike | None = None):
@app.update @app.update
def draw_energy_limit(core_diameter_um: float, pressure_mbar: float, t_fwhm_fs: float): def draw_energy_limit(core_diameter_um: float, pressure_mbar: float, t_fwhm_fs: float):
lim = compute_max_energy(core_diameter_um, pressure_mbar, t_fwhm_fs)
zdw = lim.wl_zero_disp * 1e9
app[0].set_line_data("zdw", [zdw, zdw], [-3, 3])
info_lines = [ info_lines = [
f"gas = {config.gas.title()}", f"gas = {config.gas.title()}",
f"rep rate = {config.rep_rate*1e-3}kHz", f"rep rate = {config.rep_rate*1e-3}kHz",
f"ZDW = {zdw:.0f}nm",
] ]
try:
lim = compute_max_energy(core_diameter_um, pressure_mbar, t_fwhm_fs)
except ValueError:
info_lines.append("No ZDW found")
app.info_label.setText("\n".join(info_lines))
return
zdw = lim.wl_zero_disp * 1e9
app[0].set_line_data("zdw", [zdw, zdw], [-3, 3])
info_lines.append(f"ZDW = {zdw:.0f}nm")
if lim.ion_lim > lim.sf_lim: if lim.ion_lim > lim.sf_lim:
power = lim.sf_lim * 1e3 * config.rep_rate power = lim.sf_lim * 1e3 * config.rep_rate
info_lines += [ info_lines += [