  An integral equation-based numerical method for scattering from
multi-dielectric cylinders is presented. Electromagnetic fields are represented
via layer potentials in terms of surface densities with physical
interpretations. The existence of null-field representations then adds superior
flexibility to the modeling. Local representations are used for fast field
evaluation at points away from their sources. Partially global representations,
constructed as to reduce the strength of kernel singularities, are used for
near-evaluations. A mix of local- and partially global representations is also
used to derive the system of integral equations from which the physical
densities are solved. Unique solvability is proven for the special case of
scattering from a homogeneous cylinder under rather general conditions. High
achievable accuracy is demonstrated for several examples found in the
literature.
