From 73e3fa18b11c962dffcbed80a03d522bd38e3fb2 Mon Sep 17 00:00:00 2001 From: Jan Petykiewicz Date: Tue, 17 Oct 2017 13:00:46 -0700 Subject: [PATCH] fix cylindrical operator --- fdfd_tools/waveguide.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/fdfd_tools/waveguide.py b/fdfd_tools/waveguide.py index 85e96bf..1566a74 100644 --- a/fdfd_tools/waveguide.py +++ b/fdfd_tools/waveguide.py @@ -340,8 +340,8 @@ def cylindrical_operator(omega: complex, rx = r0 + numpy.cumsum(dxes[0][0]) ry = r0 + dxes[0][0]/2.0 + numpy.cumsum(dxes[1][0]) - tx = 1 + rx/r0 - ty = 1 + ry/r0 + tx = rx/r0 + ty = ry/r0 Tx = sparse.diags(vec(tx[:, None].repeat(dxes[0][1].size, axis=1))) Ty = sparse.diags(vec(ty[:, None].repeat(dxes[1][1].size, axis=1)))