A Multigrid Homotopy Method For Computing Eigenfunctions of Differential Operators With Two-Dimensional Heterogeneity