def calc_luminance(cos_a, cos_b, cos_phi, cos_theta, sin_a, sin_b, sin_phi, sin_theta):
    """
    Luminance should be equal to:
    cosphi*costheta*sinB - cosA*costheta*sinphi -
    sinA*sintheta + cosB*(cosA*sintheta - costheta*sinA*sinphi);
    :param cos_a:
    :param cos_b:
    :param cos_phi:
    :param cos_theta:
    :param sin_a:
    :param sin_b:
    :param sin_phi:
    :param sin_theta:
    :return:
    """
    return cos_phi * cos_theta * sin_b - \
           cos_a * cos_theta * sin_phi - \
           sin_a * sin_theta + cos_b * (
                   cos_a * sin_theta - cos_theta * sin_a * sin_phi
           )