From 11e74795cb045a5364f13bc5e5447fd699c77ea3 Mon Sep 17 00:00:00 2001 From: sj90101 Date: Fri, 13 Jan 2017 10:44:09 +0800 Subject: [PATCH 1/2] [waveguide] Fix normal<->tangent confusion --- pytential/symbolic/pde/maxwell/waveguide.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/pytential/symbolic/pde/maxwell/waveguide.py b/pytential/symbolic/pde/maxwell/waveguide.py index e417584e..aa86013d 100644 --- a/pytential/symbolic/pde/maxwell/waveguide.py +++ b/pytential/symbolic/pde/maxwell/waveguide.py @@ -227,9 +227,9 @@ class SecondKindInfZMuellerOperator(L2WeightedPDEOperator): 1j * (Tt(0, phi) - Tt(1, phi)) - 1j * ( n0**2 * tangent.scalar_product( - S(0, normal * phi)) + S(0, tangent * phi)) - n1**2 * tangent.scalar_product( - S(1, normal * phi)))) + S(1, tangent * phi)))) a23 = +1*sym.cse(a23_expr(phi3), "a23") a41 = -1*sym.cse(a23_expr(phi1), "a41") -- GitLab From c1bed950b9a91c339612b67b7e5269459771bc05 Mon Sep 17 00:00:00 2001 From: sj90101 Date: Fri, 13 Jan 2017 16:02:03 +0800 Subject: [PATCH 2/2] Adjust mode find radius --- examples/find-photonic-mode-sk.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/find-photonic-mode-sk.py b/examples/find-photonic-mode-sk.py index 0e65f390..ba52dffb 100644 --- a/examples/find-photonic-mode-sk.py +++ b/examples/find-photonic-mode-sk.py @@ -82,7 +82,7 @@ def find_mode(): mesh = make_curve_mesh(curve_f, np.linspace(0, 1, nelements+1), target_order) - circle_radius = 50*k_vacuum + circle_radius = 4*k_vacuum mesh = affine_map(mesh, A=circle_radius*np.eye(2)) from meshmode.discretization import Discretization -- GitLab