  We present the results of Monte Carlo simulation for a Kondo lattice model in
which itinerant electrons interact with Ising spins with spin-ice type
easy-axis anisotropy on a pyrochlore lattice. We demonstrate the efficiency of
the truncated polynomial expansion algorithm, which enables a large scale
simulation, in comparison with a conventional algorithm using the exact
diagonalization. Computing the sublattice magnetization, we show the
convergence of the data with increasing the number of polynomials and
truncation distance.
