This is a tested implementation of the Hadamard transform . It should be quite fast since it uses the fast Hadamard transform algorithm and is implemented using python C api. You can apply the Hadamard transform to ndarray of dimensions 1, 2 or 3. It requires that input shapes are power of two since otherwise the Hadamard transform is not defined.
You need numpy for fht to run.