FreeCalypso > hg > gsm-codec-lib
view libtwamr/log2.c @ 466:0c4e1bc06740
amrefr-encode-r: use the new DHF xform that matches T-Mobile
author | Mychaela Falconia <falcon@freecalypso.org> |
---|---|
date | Sun, 12 May 2024 21:44:59 +0000 |
parents | bfe74a9edd5a |
children |
line wrap: on
line source
/* ******************************************************************************** * * GSM AMR-NB speech codec R98 Version 7.6.0 December 12, 2001 * R99 Version 3.3.0 * REL-4 Version 4.1.0 * ******************************************************************************** * * File : log2.c * Purpose : Computes log2(L_x) * ******************************************************************************** */ /* ******************************************************************************** * MODULE INCLUDE FILE AND VERSION ID ******************************************************************************** */ #include "namespace.h" #include "log2.h" /* ******************************************************************************** * INCLUDE FILES ******************************************************************************** */ #include "typedef.h" #include "basic_op.h" #include "no_count.h" /* ******************************************************************************** * LOCAL VARIABLES AND TABLES ******************************************************************************** */ #include "log2.tab" /* Table for Log2() */ /* ******************************************************************************** * PUBLIC PROGRAM CODE ******************************************************************************** */ /************************************************************************* * * FUNCTION: Log2_norm() * * PURPOSE: Computes log2(L_x, exp), where L_x is positive and * normalized, and exp is the normalisation exponent * If L_x is negative or zero, the result is 0. * * DESCRIPTION: * The function Log2(L_x) is approximated by a table and linear * interpolation. The following steps are used to compute Log2(L_x) * * 1- exponent = 30-norm_exponent * 2- i = bit25-b31 of L_x; 32<=i<=63 (because of normalization). * 3- a = bit10-b24 * 4- i -=32 * 5- fraction = table[i]<<16 - (table[i] - table[i+1]) * a * 2 * *************************************************************************/ void Log2_norm ( Word32 L_x, /* (i) : input value (normalized) */ Word16 exp, /* (i) : norm_l (L_x) */ Word16 *exponent, /* (o) : Integer part of Log2. (range: 0<=val<=30) */ Word16 *fraction /* (o) : Fractional part of Log2. (range: 0<=val<1) */ ) { Word16 i, a, tmp; Word32 L_y; test (); if (L_x <= (Word32) 0) { *exponent = 0; move16 (); *fraction = 0; move16 (); return; } *exponent = sub (30, exp); move16 (); L_x = L_shr (L_x, 9); i = extract_h (L_x); /* Extract b25-b31 */ L_x = L_shr (L_x, 1); a = extract_l (L_x); /* Extract b10-b24 of fraction */ a = a & (Word16) 0x7fff; logic16 (); i = sub (i, 32); L_y = L_deposit_h (table[i]); /* table[i] << 16 */ tmp = sub (table[i], table[i + 1]); /* table[i] - table[i+1] */ L_y = L_msu (L_y, tmp, a); /* L_y -= tmp*a*2 */ *fraction = extract_h (L_y);move16 (); return; } /************************************************************************* * * FUNCTION: Log2() * * PURPOSE: Computes log2(L_x), where L_x is positive. * If L_x is negative or zero, the result is 0. * * DESCRIPTION: * normalizes L_x and then calls Log2_norm(). * *************************************************************************/ void Log2 ( Word32 L_x, /* (i) : input value */ Word16 *exponent, /* (o) : Integer part of Log2. (range: 0<=val<=30) */ Word16 *fraction /* (o) : Fractional part of Log2. (range: 0<=val<1) */ ) { Word16 exp; exp = norm_l (L_x); Log2_norm (L_shl (L_x, exp), exp, exponent, fraction); }