“…In its original formulation [5], the DSA method was only applicable to problems involving infinite circular and rectangular cylinders, without magnetic contrast. Since then, multiple extensions have been conceived: to determine the per-unit-of-length parameters of transmission lines including semiconducting materials [6], to triangular cross-sections [7], to threedimensional (3-D) cylinders and cuboids [8,9,10], and, very recently, to situations including combined dielectric and magnetic contrast [11]. The traditional approach to construct the DSA operator requires the Dirichlet eigenfunctions of the modeled structure, enforcing a restriction to canonical shapes, while the extension presented in [7] involves specific modifications to avoid a significant Gibbs effect, and requires a composition of multiple triangular parts to model arbitrary polygonal shapes.…”