> The beauty here is that eml(x,y) is a pure, continuous analytical function with no hidden branching whatsoever.
They use the complex version of logarithm, that has a lot of branching problems.
They use the complex version of logarithm, that has a lot of branching problems.
So it isn't exploiting the branching for computation.
For example, IIRC ln( -inf.0 + y * i ) = ´+inf.0 + pi * sign(y)