31 color cosi2 = color(cosi * cosi);
32 color one = color(1, 1, 1);
33 color tmp_f = eta * eta + k * k;
34 color tmp = tmp_f * cosi2;
35 color Rparl2 = (tmp - (2.0 * eta * cosi) + one) / (tmp + (2.0 * eta * cosi) + one);
36 color Rperp2 = (tmp_f - (2.0 * eta * cosi) + cosi2) / (tmp_f + (2.0 * eta * cosi) + cosi2);
37 return (Rparl2 + Rperp2) * 0.5;