Exact computation of the color function for triangular element interfaces