A new capability that enables annual simulation of optically complex fenestration systems has been added to Radiance. The method relies on bidirectional scattering distribution function (BSDF) input data, which are used in an efficient matrix calculation to compute time-step performance given TMY data. The objective of this study was to explain the value of this capability to designers and developers of innovative daylighting systems and to demonstrate its speed and accuracy via comparisons of simulated to measured illuminance data for a daylight-redirecting optical louver system. The method was shown to provide valid results that accurately replicate real-world conditions with an absolute mean bias error below 13% and a root mean square error below 23%. Routine application of this new capability will not be hindered by slow computational speed for illuminance calculations. Instead, the capability will be dependent on the availability of BSDF data for daylighting, shading and fenestration systems.