Radiation from stars and black holes strongly affects their surroundings and plays a crucial role in topics such as stellar atmospheres, the interstellar medium, star formation, galaxy formation, supernovae, and cosmology.  It is a well-studied problem, but its treatment in multidimensional calculations is difficult because of the dependence on seven variables — three spatial, two angular, frequency, and time.  We have successfully implemented and extensively utilized a method of solving the radiative transfer equation inside the adaptive mesh refinement code, Enzo.  We are still actively developing better algorithms to solve this equation more accurately and efficiently in numerical simulations.