gauthier
c7518ecf42
git-svn-id: http://svn.eurecom.fr/openair4G/trunk@7885 818b1a75-f10b-46b9-bf7c-635c3b92a50f
649 lines
6.4 KiB
C
Executable File
649 lines
6.4 KiB
C
Executable File
/*******************************************************************************
|
|
OpenAirInterface
|
|
Copyright(c) 1999 - 2014 Eurecom
|
|
|
|
OpenAirInterface 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 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
|
|
OpenAirInterface 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 OpenAirInterface.The full GNU General Public License is
|
|
included in this distribution in the file called "COPYING". If not,
|
|
see <http://www.gnu.org/licenses/>.
|
|
|
|
Contact Information
|
|
OpenAirInterface Admin: openair_admin@eurecom.fr
|
|
OpenAirInterface Tech : openair_tech@eurecom.fr
|
|
OpenAirInterface Dev : openair4g-devel@lists.eurecom.fr
|
|
|
|
Address : Eurecom, Campus SophiaTech, 450 Route des Chappes, CS 50193 - 06904 Biot Sophia Antipolis cedex, FRANCE
|
|
|
|
*******************************************************************************/
|
|
#include "defs.h"
|
|
|
|
// Approximate 10*log10(x) in fixed point : x = 0...(2^32)-1
|
|
|
|
int8_t dB_table[256] = {
|
|
0,
|
|
3,
|
|
5,
|
|
6,
|
|
7,
|
|
8,
|
|
8,
|
|
9,
|
|
10,
|
|
10,
|
|
10,
|
|
11,
|
|
11,
|
|
11,
|
|
12,
|
|
12,
|
|
12,
|
|
13,
|
|
13,
|
|
13,
|
|
13,
|
|
13,
|
|
14,
|
|
14,
|
|
14,
|
|
14,
|
|
14,
|
|
14,
|
|
15,
|
|
15,
|
|
15,
|
|
15,
|
|
15,
|
|
15,
|
|
15,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
16,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
17,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
18,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
19,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
20,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
21,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
22,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
23,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24,
|
|
24
|
|
};
|
|
|
|
int16_t dB_table_times10[256] = {
|
|
0,
|
|
30,
|
|
47,
|
|
60,
|
|
69,
|
|
77,
|
|
84,
|
|
90,
|
|
95,
|
|
100,
|
|
104,
|
|
107,
|
|
111,
|
|
114,
|
|
117,
|
|
120,
|
|
123,
|
|
125,
|
|
127,
|
|
130,
|
|
132,
|
|
134,
|
|
136,
|
|
138,
|
|
139,
|
|
141,
|
|
143,
|
|
144,
|
|
146,
|
|
147,
|
|
149,
|
|
150,
|
|
151,
|
|
153,
|
|
154,
|
|
155,
|
|
156,
|
|
157,
|
|
159,
|
|
160,
|
|
161,
|
|
162,
|
|
163,
|
|
164,
|
|
165,
|
|
166,
|
|
167,
|
|
168,
|
|
169,
|
|
169,
|
|
170,
|
|
171,
|
|
172,
|
|
173,
|
|
174,
|
|
174,
|
|
175,
|
|
176,
|
|
177,
|
|
177,
|
|
178,
|
|
179,
|
|
179,
|
|
180,
|
|
181,
|
|
181,
|
|
182,
|
|
183,
|
|
183,
|
|
184,
|
|
185,
|
|
185,
|
|
186,
|
|
186,
|
|
187,
|
|
188,
|
|
188,
|
|
189,
|
|
189,
|
|
190,
|
|
190,
|
|
191,
|
|
191,
|
|
192,
|
|
192,
|
|
193,
|
|
193,
|
|
194,
|
|
194,
|
|
195,
|
|
195,
|
|
196,
|
|
196,
|
|
197,
|
|
197,
|
|
198,
|
|
198,
|
|
199,
|
|
199,
|
|
200,
|
|
200,
|
|
200,
|
|
201,
|
|
201,
|
|
202,
|
|
202,
|
|
202,
|
|
203,
|
|
203,
|
|
204,
|
|
204,
|
|
204,
|
|
205,
|
|
205,
|
|
206,
|
|
206,
|
|
206,
|
|
207,
|
|
207,
|
|
207,
|
|
208,
|
|
208,
|
|
208,
|
|
209,
|
|
209,
|
|
210,
|
|
210,
|
|
210,
|
|
211,
|
|
211,
|
|
211,
|
|
212,
|
|
212,
|
|
212,
|
|
213,
|
|
213,
|
|
213,
|
|
213,
|
|
214,
|
|
214,
|
|
214,
|
|
215,
|
|
215,
|
|
215,
|
|
216,
|
|
216,
|
|
216,
|
|
217,
|
|
217,
|
|
217,
|
|
217,
|
|
218,
|
|
218,
|
|
218,
|
|
219,
|
|
219,
|
|
219,
|
|
219,
|
|
220,
|
|
220,
|
|
220,
|
|
220,
|
|
221,
|
|
221,
|
|
221,
|
|
222,
|
|
222,
|
|
222,
|
|
222,
|
|
223,
|
|
223,
|
|
223,
|
|
223,
|
|
224,
|
|
224,
|
|
224,
|
|
224,
|
|
225,
|
|
225,
|
|
225,
|
|
225,
|
|
226,
|
|
226,
|
|
226,
|
|
226,
|
|
226,
|
|
227,
|
|
227,
|
|
227,
|
|
227,
|
|
228,
|
|
228,
|
|
228,
|
|
228,
|
|
229,
|
|
229,
|
|
229,
|
|
229,
|
|
229,
|
|
230,
|
|
230,
|
|
230,
|
|
230,
|
|
230,
|
|
231,
|
|
231,
|
|
231,
|
|
231,
|
|
232,
|
|
232,
|
|
232,
|
|
232,
|
|
232,
|
|
233,
|
|
233,
|
|
233,
|
|
233,
|
|
233,
|
|
234,
|
|
234,
|
|
234,
|
|
234,
|
|
234,
|
|
235,
|
|
235,
|
|
235,
|
|
235,
|
|
235,
|
|
235,
|
|
236,
|
|
236,
|
|
236,
|
|
236,
|
|
236,
|
|
237,
|
|
237,
|
|
237,
|
|
237,
|
|
237,
|
|
238,
|
|
238,
|
|
238,
|
|
238,
|
|
238,
|
|
238,
|
|
239,
|
|
239,
|
|
239,
|
|
239,
|
|
239,
|
|
239,
|
|
240,
|
|
240,
|
|
240,
|
|
240,
|
|
240
|
|
};
|
|
/*
|
|
int8_t dB_fixed(int x) {
|
|
|
|
int i=0,adj=0;
|
|
int8_t log10=0;
|
|
|
|
// find MSB
|
|
for (i=31;i>=0;i--) {
|
|
|
|
if ((x & (1<<i)) >= 1) {
|
|
log10 = 3*i;
|
|
i=0;
|
|
}
|
|
|
|
}
|
|
|
|
// look at next 2 MSBs and adjust between 0-2
|
|
if (i>1) {
|
|
adj = (x>>(i-2))&3;
|
|
if (adj == 1)
|
|
log10 += 1;
|
|
else if ((adj == 2) || (adj == 3))
|
|
log10 += 2;
|
|
}
|
|
else
|
|
log10 += (x&1 == 1) ? 1 : 0;
|
|
|
|
return(log10);
|
|
}
|
|
*/
|
|
|
|
int16_t dB_fixed_times10(uint32_t x)
|
|
{
|
|
int16_t dB_power=0;
|
|
|
|
|
|
if (x==0) {
|
|
dB_power = 0;
|
|
} else if ( (x&0xff000000) != 0 ) {
|
|
dB_power = dB_table_times10[((x>>24)&255)-1];
|
|
dB_power += 3*dB_table_times10[255];
|
|
} else if ( (x&0x00ff0000) != 0 ) {
|
|
dB_power = dB_table_times10[((x>>16)&255)-1];
|
|
dB_power += 2*dB_table_times10[255];
|
|
} else if ( (x&0x0000ff00) != 0 ) {
|
|
dB_power = dB_table_times10[((x>>8)&255)-1];
|
|
dB_power += dB_table_times10[255];
|
|
} else {
|
|
dB_power = dB_table_times10[(x&255)-1];
|
|
}
|
|
|
|
if (dB_power > 900)
|
|
return(900);
|
|
|
|
return dB_power;
|
|
}
|
|
|
|
int8_t dB_fixed(uint32_t x)
|
|
{
|
|
|
|
int8_t dB_power=0;
|
|
|
|
|
|
if (x==0) {
|
|
dB_power = 0;
|
|
} else if ( (x&0xff000000) != 0 ) {
|
|
dB_power = dB_table[((x>>24)&255)-1];
|
|
dB_power += 3*dB_table[255];
|
|
} else if ( (x&0x00ff0000) != 0 ) {
|
|
dB_power = dB_table[((x>>16)&255)-1];
|
|
dB_power += 2*dB_table[255];
|
|
} else if ( (x&0x0000ff00) != 0 ) {
|
|
dB_power = dB_table[((x>>8)&255)-1];
|
|
dB_power += dB_table[255];
|
|
} else {
|
|
dB_power = dB_table[(x&255)-1];
|
|
}
|
|
|
|
if (dB_power > 90)
|
|
return(90);
|
|
|
|
return dB_power;
|
|
}
|
|
|
|
int8_t dB_fixed2(uint32_t x, uint32_t y)
|
|
{
|
|
|
|
if ((x>0) && (y>0) )
|
|
if (x>y)
|
|
return(dB_fixed(x/y));
|
|
else
|
|
return(-dB_fixed(y/x));
|
|
else if (y==0)
|
|
return(127);
|
|
else if (x==0)
|
|
return(-128);
|
|
|
|
return(0);
|
|
}
|