// Dca' = 2.Dca < Da ? // 2.Sca.Dca + Sca.(1 - Da) + Dca.(1 - Sa) // Sa.Da - 2.(Da - Dca).(Sa - Sca) + Sca.(1 - Da) + Dca.(1 - Sa) // Da' = Sa + Da - Sa.Da vec4 composite(vec4 src, vec4 dst) { vec4 result; result.rgb = mix(2 * src.rgb * dst.rgb + src.rgb * (1 - dst.a) + dst.rgb * (1 - src.a), src.a * dst.a - 2 * (dst.a - dst.rgb) * (src.a - src.rgb) + src.rgb * (1 - dst.a) + dst.rgb * (1 - src.a), step(dst.a, 2 * dst.rgb)); result.a = src.a + dst.a - src.a * dst.a; return result; }