float cosSqOmega = cosOmega * cosOmega;
[branch]
if (cosSqOmega > sinSqSigma)
if (cosSqOmega >= sinSqSigma)
{
irradiance = sinSqSigma * saturate(cosOmega);
}