Modeling solar heat gains of complex fenestration systems using bidirectional scattering distribution functions