1
0
mirror of https://github.com/aircrack-ng/rtl8812au.git synced 2024-11-25 14:44:09 +00:00
rtl8812au/hal/phydm/phydm_math_lib.c
2022-12-19 22:53:39 +01:00

291 lines
6.1 KiB
C

/******************************************************************************
*
* Copyright(c) 2007 - 2017 Realtek Corporation.
*
* This program is free software; you can redistribute it and/or modify it
* under the terms of version 2 of the GNU General Public License as
* published by the Free Software Foundation.
*
* 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.
*
* The full GNU General Public License is included in this distribution in the
* file called LICENSE.
*
* Contact Information:
* wlanfae <wlanfae@realtek.com>
* Realtek Corporation, No. 2, Innovation Road II, Hsinchu Science Park,
* Hsinchu 300, Taiwan.
*
* Larry Finger <Larry.Finger@lwfinger.net>
*
*****************************************************************************/
/*@************************************************************
* include files
************************************************************/
#include "mp_precomp.h"
#include "phydm_precomp.h"
const u32 db_invert_table[12][8] = {
{10, 13, 16, 20, 25, 32, 40, 50}, /* @U(32,3) */
{64, 80, 101, 128, 160, 201, 256, 318}, /* @U(32,3) */
{401, 505, 635, 800, 1007, 1268, 1596, 2010}, /* @U(32,3) */
{316, 398, 501, 631, 794, 1000, 1259, 1585}, /* @U(32,0) */
{1995, 2512, 3162, 3981, 5012, 6310, 7943, 10000}, /* @U(32,0) */
{12589, 15849, 19953, 25119, 31623, 39811, 50119, 63098}, /* @U(32,0) */
{79433, 100000, 125893, 158489, 199526, 251189, 316228,
398107}, /* @U(32,0) */
{501187, 630957, 794328, 1000000, 1258925, 1584893, 1995262,
2511886}, /* @U(32,0) */
{3162278, 3981072, 5011872, 6309573, 7943282, 1000000, 12589254,
15848932}, /* @U(32,0) */
{19952623, 25118864, 31622777, 39810717, 50118723, 63095734,
79432823, 100000000}, /* @U(32,0) */
{125892541, 158489319, 199526232, 251188643, 316227766, 398107171,
501187234, 630957345}, /* @U(32,0) */
{794328235, 1000000000, 1258925412, 1584893192, 1995262315,
2511886432U, 3162277660U, 3981071706U} }; /* @U(32,0) */
/*Y = 10*log(X)*/
s32 odm_pwdb_conversion(s32 X, u32 total_bit, u32 decimal_bit)
{
s32 Y, integer = 0, decimal = 0;
u32 i;
if (X == 0)
X = 1; /* @log2(x), x can't be 0 */
for (i = (total_bit - 1); i > 0; i--) {
if (X & BIT(i)) {
integer = i;
if (i > 0) {
/*decimal is 0.5dB*3=1.5dB~=2dB */
decimal = (X & BIT(i - 1)) ? 2 : 0;
}
break;
}
}
Y = 3 * (integer - decimal_bit) + decimal; /* @10*log(x)=3*log2(x), */
return Y;
}
s32 odm_sign_conversion(s32 value, u32 total_bit)
{
if (value & BIT(total_bit - 1))
value -= BIT(total_bit);
return value;
}
/*threshold must form low to high*/
u16 phydm_find_intrvl(void *dm_void, u16 val, u16 *threshold, u16 th_len)
{
struct dm_struct *dm = (struct dm_struct *)dm_void;
u16 i = 0;
u16 ret_val = 0;
u16 max_th = threshold[th_len - 1];
for (i = 0; i < th_len; i++) {
if (val < threshold[i]) {
ret_val = i;
break;
} else if (val >= max_th) {
ret_val = th_len;
break;
}
}
return ret_val;
}
void phydm_seq_sorting(void *dm_void, u32 *value, u32 *rank_idx, u32 *idx_out,
u8 seq_length)
{
u8 i = 0, j = 0;
u32 tmp_a, tmp_b;
u32 tmp_idx_a, tmp_idx_b;
for (i = 0; i < seq_length; i++)
rank_idx[i] = i;
for (i = 0; i < (seq_length - 1); i++) {
for (j = 0; j < (seq_length - 1 - i); j++) {
tmp_a = value[j];
tmp_b = value[j + 1];
tmp_idx_a = rank_idx[j];
tmp_idx_b = rank_idx[j + 1];
if (tmp_a < tmp_b) {
value[j] = tmp_b;
value[j + 1] = tmp_a;
rank_idx[j] = tmp_idx_b;
rank_idx[j + 1] = tmp_idx_a;
}
}
}
for (i = 0; i < seq_length; i++)
idx_out[rank_idx[i]] = i + 1;
}
u32 odm_convert_to_db(u64 value)
{
u8 i;
u8 j;
u32 dB;
if (value >= db_invert_table[11][7])
return 96; /* @maximum 96 dB */
for (i = 0; i < 12; i++) {
if (i <= 2 && (value << FRAC_BITS) <= db_invert_table[i][7])
break;
else if (i > 2 && value <= db_invert_table[i][7])
break;
}
for (j = 0; j < 8; j++) {
if (i <= 2 && (value << FRAC_BITS) <= db_invert_table[i][j])
break;
else if (i > 2 && i < 12 && value <= db_invert_table[i][j])
break;
}
/*special cases*/
if (j == 0 && i == 0)
goto end;
if (i == 3 && j == 0) {
if (db_invert_table[3][0] - value >
value - (db_invert_table[2][7] >> FRAC_BITS)) {
i = 2;
j = 7;
}
goto end;
}
if (i < 3)
value = value << FRAC_BITS; /*@elements of row 0~2 shift left*/
/*compare difference to get precise dB*/
if (j == 0) {
if (db_invert_table[i][j] - value >
value - db_invert_table[i - 1][7]) {
i = i - 1;
j = 7;
}
} else {
if (db_invert_table[i][j] - value >
value - db_invert_table[i][j - 1]) {
j = j - 1;
}
}
end:
dB = (i << 3) + j + 1;
return dB;
}
u64 phydm_db_2_linear(u32 value)
{
u8 i = 0;
u8 j = 0;
u64 linear = 0;
value = value & 0xFF;
/* @1dB~96dB */
if (value > 96) {
value = 96;
} else if (value < 1) {
linear = 1;
return linear;
}
i = (u8)((value - 1) >> 3);
j = (u8)(value - 1) - (i << 3);
linear = db_invert_table[i][j];
if (i > 2)
linear = linear << FRAC_BITS;
return linear;
}
u16 phydm_show_fraction_num(u32 frac_val, u8 bit_num)
{
u8 i = 0;
u16 val = 0;
u16 base = 5000;
for (i = bit_num; i > 0; i--) {
if (frac_val & BIT(i - 1))
val += (base >> (bit_num - i));
}
return val;
}
u16 phydm_ones_num_in_bitmap(u64 val, u8 size)
{
u8 i = 0;
u8 ones_num = 0;
for (i = 0; i < size; i++) {
if (val & BIT(0))
ones_num++;
val = val >> 1;
}
return ones_num;
}
u64 phydm_gen_bitmask(u8 mask_num)
{
u8 i = 0;
u64 bitmask = 0;
if (mask_num > 64)
return 1;
for (i = 0; i < mask_num; i++)
bitmask = (bitmask << 1) | BIT(0);
return bitmask;
}
s32 phydm_cnvrt_2_sign(u32 val, u8 bit_num)
{
if (bit_num >= 32)
return (s32)val;
if (val & BIT(bit_num - 1)) /*Sign BIT*/
val -= (1 << bit_num); /*@2's*/
return val;
}
s64 phydm_cnvrt_2_sign_64(u64 val, u8 bit_num)
{
u64 one = 1;
s64 val_sign = (s64)val;
if (bit_num >= 64)
return (s64)val;
if (val & (one << (bit_num - 1))) /*Sign BIT*/
val_sign = val - (one << bit_num); /*@2's*/
return val_sign;
}