  We perform a comparison between two radiative transfer algorithms commonly
employed in hydrodynamical calculations of star formation: grey flux limited
diffusion and the hybrid scheme, in addition we compare these algorithms to
results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like
density structures the hybrid scheme performs significantly better than the FLD
method in the optically thin regions, with comparable results in optically
thick regions. In the case of a forming high mass star we find the FLD method
significantly underestimates the radiation pressure by a factor of ~100.
