This is a sndkit DSP algorithm. A more up-to-date version can be found on the sndkit wiki here.
DSP
sndkit