Page 1 of 1
What DSP code is this ASM rms calc?
Posted: Wed Aug 15, 2018 5:29 am
by nix
Hiya 8D
I have been intensively using Bootsy's fast RMS asm code.
I was wondering if anyone would mind writing it in DSP for me?
Here's the code->
Cheers, I looked at the Wiki,
but the DSP implementation might help me understand.
It's a wonderful algo-
ta Bootsy!
Re: What DSP code is this ASM rms calc?
Posted: Wed Aug 15, 2018 1:45 pm
by martinvicanek
Code: Select all
streamin in;
streamin c;
streamout out;
float rms;
rms = rms + c*(in*in - rms);
out = sqrt(rms);
You may get denormals (in both the ASM and DSP) if you feed it with synthetisised audio. Easy fix:
Code: Select all
streamin in;
streamin c;
streamout out;
float rms;
float denorm=1e-11;
rms = rms + c*(in*in - rms) + denorm - denorm;
out = sqrt(rms);
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 1:24 am
by tulamide
This is interesting for another reason.
Martin, as a guru of DSP, what do you think how long will it take until this code outputs meaningful values? I compare it to the origin of rms, which is taking the root mean square of the whole audio. On old consoles and early computer measurings of realtime audio, they measured in windows of certain sizes. This code seems to not have any window size, so I wonder about it - thus the question.
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 3:53 am
by nix
Thanks so much Martin!
will test it.
'c' is the window Tula, you need the original samples conversion in the first upload too
2205 samples @ 44.1k is accurate enough. longer is better precision, and shorter is not delayed as much
Cheers folk!
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 5:28 am
by martinvicanek
tulamide wrote:how long will it take until this code outputs meaningful values? [...] This code seems to not have any window size
The code contains a so-called leaky integrator for the squared signal. You can think of it as a sliding average with an asymmetric (exponential) window. As Nix points out, the c parameter (or rather its inverse) denotes the effective window size.
To your question: Yes, there will be an exponentially decaying transient at the beginning. After a few window lengths the transient is gone and you will get "meaningful" readings.
A rectangular sliding window would require only one window length to reach steady state. A sliding window consumes a bit more CPU than a leaky integrator.
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 10:07 am
by nix
So the square root makes it more accurate maybe?
-also, could we rectify instead of- in*in
then we wouldn't need the square root?
ta Martin for helping with the integrator concept
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 11:34 am
by martinvicanek
nix wrote:could we rectify instead of- in*in
then we wouldn't need the square root?
Sure, but then it wouldnt be exactly RMS. In many cases it does not make a difference, though, and you can get away with a rectifier.
Re: What DSP code is this ASM rms calc?
Posted: Thu Aug 16, 2018 1:55 pm
by tulamide
Ah, thank you both! I understand the code better now, esp. the difference between rectangular window and leaky integrator. I come to the conclusion that I often use techniques in Ruby (my peak meter works with rectangular window), that I don't know the names of (or even think I came up with it first
).