Several different numerical methods for calculating diffuse radiation viewfactors are described. Each is applied to a range of surface configurations, from almost completely unobstructed to a dense set of intersecting surfaces. The speed, accuracy, and unique characteristics are discussed in order to define optimal methods for different surface geometries.