Suppose that $L$ is a parabolic differential operator, is there a stable numerical scheme for finding eigen-functions of $L$? Here the domain of interest is a closed, bounded, convex subset of $\mathbb{R}^d$.
In fact, I'm looking for the eigen-function corresponding to the minimal eigenvalue, with Dirichlet boundary condition on the boundary of $D$.