GELU really is like magic:
UNARY(GELU, b / 2 * (1 + tanh(.7978845 * (b + .044715 * b * b * b))))
UNARY(GELU, b / 2 * (1 + tanh(.7978845 * (b + .044715 * b * b * b))))
In contrast, the fast inverse square root really exploits the bit representation of a floating point input to cheaply compute an initial guess.