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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
/*
* Copyright (C) 2012 Fons Adriaensen <fons@linuxaudio.org>
* Copyright (C) 2013-2019 Robin Gareus <robin@gareus.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, write to the Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*/
#include <math.h>
#include "ardour/iec1ppmdsp.h"
float Iec1ppmdsp::_w1;
float Iec1ppmdsp::_w2;
float Iec1ppmdsp::_w3;
float Iec1ppmdsp::_g;
Iec1ppmdsp::Iec1ppmdsp (void)
: _z1 (0)
, _z2 (0)
, _m (0)
, _res (true)
{}
Iec1ppmdsp::~Iec1ppmdsp (void) {}
void
Iec1ppmdsp::process (float const* p, int n)
{
float z1, z2, m, t;
z1 = _z1 > 20 ? 20 : (_z1 < 0 ? 0 : _z1);
z2 = _z2 > 20 ? 20 : (_z2 < 0 ? 0 : _z2);
m = _res ? 0: _m;
_res = false;
n /= 4;
while (n--) {
z1 *= _w3;
z2 *= _w3;
t = fabsf (*p++);
if (t > z1) z1 += _w1 * (t - z1);
if (t > z2) z2 += _w2 * (t - z2);
t = fabsf (*p++);
if (t > z1) z1 += _w1 * (t - z1);
if (t > z2) z2 += _w2 * (t - z2);
t = fabsf (*p++);
if (t > z1) z1 += _w1 * (t - z1);
if (t > z2) z2 += _w2 * (t - z2);
t = fabsf (*p++);
if (t > z1) z1 += _w1 * (t - z1);
if (t > z2) z2 += _w2 * (t - z2);
t = z1 + z2;
if (t > m) m = t;
}
_z1 = z1 + 1e-10f;
_z2 = z2 + 1e-10f;
_m = m;
}
float
Iec1ppmdsp::read (void)
{
_res = true;
return _g * _m;
}
void
Iec1ppmdsp::reset ()
{
_z1 = _z2 = _m = .0f;
_res = true;
}
void
Iec1ppmdsp::init (float fsamp)
{
_w1 = 450.0f / fsamp;
_w2 = 1300.0f / fsamp;
_w3 = 1.0f - 5.4f / fsamp;
_g = 0.5108f;
}
|
