initial guess a bit higher than highest epsilon, fix comparison bench…

…mark
HelgeGehring · Jun 17, 2023 · c5fbaf2 · c5fbaf2
1 parent ce867b7
commit c5fbaf2
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/docs/benchmarks/mode_solver.py b/docs/benchmarks/mode_solver.py
@@ -87,7 +87,7 @@
     for subdomain, n in {"core": 3.44, "box": 3.40, "clad": 1}.items():
         epsilon[basis0.get_dofs(elements=subdomain)] = n**2
 
-    modes = compute_modes(basis0, epsilon, wavelength=1.15, num_modes=1, order=2)
+    modes = compute_modes(basis0, epsilon, wavelength=1.15, num_modes=1, order=2, solver="slepc")
     neff_values_femwell_slepc.append(np.real(modes[0].n_eff))
 
     modes = compute_modes(basis0, epsilon, wavelength=1.15, num_modes=1, order=2, solver="scipy")

diff --git a/femwell/maxwell/waveguide.py b/femwell/maxwell/waveguide.py
@@ -241,7 +241,7 @@ def bform(e_t, e_z, v_t, v_z, w):
     if n_guess:
         sigma = sigma = k0**2 * n_guess**2
     else:
-        sigma = sigma = k0**2 * np.max(epsilon_r)
+        sigma = sigma = k0**2 * np.max(epsilon_r) * 1.1
 
     if metallic_boundaries:
         lams, xs = solve(