File: itu468detect.cc

package info (click to toggle)
jnoisemeter 0.2.2-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 204 kB
  • sloc: cpp: 1,132; makefile: 34
file content (62 lines) | stat: -rw-r--r-- 1,473 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
// ----------------------------------------------------------------------------
//
//  Copyright (C) 2010-2018 Fons Adriaensen <fons@linuxaudio.org>
//    
//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2 of the License, or
//  (at your option) any later version.
//
//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with this program.  If not, see <http://www.gnu.org/licenses/>.
//
// ----------------------------------------------------------------------------


#include <math.h>
#include "itu468detect.h"


int Itu468detect::init (int fsamp)
{
    reset ();

    _a1 = 670.0f / fsamp;
    _b1 = 3.50f / fsamp;
    _a2 = 6.60f / fsamp;
    _b2 = 0.65f / fsamp;

    return 0;
}


void Itu468detect::reset (void)
{
    _z1 = _z2 = 0;
}


void Itu468detect::process (int n, const float *inp)
{
    float x, z1, z2;

    z1 = _z1;
    z2 = _z2;

    while (n--)
    {
	x = fabsf (*inp++) + 1e-30f;
	z1 -= z1 * _b1;
	if (x > z1) z1 += _a1 * (x - z1);
	z2 -= z2 * _b2;
	if (z1 > z2) z2 += _a2 * (z1 - z2);
    }

    _z1 = z1;
    _z2 = z2;
}