Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.