This work is concerned with numerical methods for matrix eigenvalue problems that are nonlinear in the eigenvalue parameter. In particular, we focus on eigenvalue problems for which the evaluation of the matrix valued function is computationally expensive. Such problems arise, e.g., from boundary integral formulations of elliptic PDE-eigenvalue problems and typically exclude the use of establis...