[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.