[Add] Heaviside Operator

Context

This will be useful for SNNs, as the relationship between spiking S and potential U can be described as :

S[t] = \Theta(U[t] - U_{\rm thr}) where \Theta(\cdot) is the Heaviside step function.

See also pytorch operator

Fixes #173.