Computing leaky modes of optical fibers using a FEAST algorithm for polynomial eigenproblems

An efficient contour integral technique to approximate a cluster of nonlinear eigenvalues of a polynomial eigenproblem, circumventing certain large