mathdp31.c
上传用户:cxx_68
上传日期:2021-02-21
资源大小:161k
文件大小:7k
源码类别:

语音压缩

开发平台:

Visual C++

  1. /*
  2. 2.4 kbps MELP Proposed Federal Standard speech coder
  3. Fixed-point C code, version 1.0
  4. Copyright (c) 1998, Texas Instruments, Inc.  
  5. Texas Instruments has intellectual property rights on the MELP
  6. algorithm.  The Texas Instruments contact for licensing issues for
  7. commercial and non-government use is William Gordon, Director,
  8. Government Contracts, Texas Instruments Incorporated, Semiconductor
  9. Group (phone 972 480 7442).
  10. The fixed-point version of the voice codec Mixed Excitation Linear
  11. Prediction (MELP) is based on specifications on the C-language software
  12. simulation contained in GSM 06.06 which is protected by copyright and
  13. is the property of the European Telecommunications Standards Institute
  14. (ETSI). This standard is available from the ETSI publication office
  15. tel. +33 (0)4 92 94 42 58. ETSI has granted a license to United States
  16. Department of Defense to use the C-language software simulation contained
  17. in GSM 06.06 for the purposes of the development of a fixed-point
  18. version of the voice codec Mixed Excitation Linear Prediction (MELP).
  19. Requests for authorization to make other use of the GSM 06.06 or
  20. otherwise distribute or modify them need to be addressed to the ETSI
  21. Secretariat fax: +33 493 65 47 16.
  22. */
  23. /***************************************************************************
  24.  *
  25.  *   File Name:  mathdp31.c
  26.  *
  27.  *   Purpose:  Contains functions increased-precision arithmetic operations.
  28.  *
  29.  *      Below is a listing of all the functions in this file.  There
  30.  *      is no interdependence among the functions.
  31.  *
  32.  *      L_mpy_ls()
  33.  *      L_mpy_ll()
  34.  *      isLwLimit()
  35.  *      isSwLimit()
  36.  *
  37.  ***************************************************************************/
  38. /*_________________________________________________________________________
  39.  |                                                                         |
  40.  |                              Include Files                              |
  41.  |_________________________________________________________________________|
  42. */
  43. #include "mathhalf.h"
  44. #include "typedefs.h"
  45. /****************************************************************************
  46.  *
  47.  *     FUNCTION NAME: isLwLimit
  48.  *
  49.  *     PURPOSE:
  50.  *
  51.  *        Check to see if the input Longword is at the
  52.  *        upper or lower limit of its range.  i.e.
  53.  *        0x7fff ffff or -0x8000 0000
  54.  *
  55.  *        Ostensibly this is a check for an overflow.
  56.  *        This does not truly mean an overflow occurred,
  57.  *        it means the value reached is the
  58.  *        maximum/minimum value representable.  It may
  59.  *        have come about due to an overflow.
  60.  *
  61.  *     INPUTS:
  62.  *
  63.  *       L_In               A Longword input variable
  64.  *
  65.  *
  66.  *     OUTPUTS:             none
  67.  *
  68.  *     RETURN VALUE:        1 if input == 0x7fff ffff or -0x8000 0000
  69.  *                          0 otherwise
  70.  *
  71.  *     KEYWORDS: saturation, limit
  72.  *
  73.  ***************************************************************************/
  74. short  isLwLimit(Longword L_In)
  75. {
  76.   Longword L_ls;
  77.   short  siOut;
  78.   if (L_In != 0)
  79.   {
  80.     L_ls = L_shl(L_In, 1);
  81.     if (L_sub(L_In, L_ls) == 0)
  82.       siOut = 1;
  83.     else
  84.       siOut = 0;
  85.   }
  86.   else
  87.   {
  88.     siOut = 0;
  89.   }
  90.   return (siOut);
  91. }
  92. /****************************************************************************
  93.  *
  94.  *     FUNCTION NAME: isSwLimit
  95.  *
  96.  *     PURPOSE:
  97.  *
  98.  *        Check to see if the input Shortword is at the
  99.  *        upper or lower limit of its range.  i.e.
  100.  *        0x7fff or -0x8000
  101.  *
  102.  *        Ostensibly this is a check for an overflow.
  103.  *        This does not truly mean an overflow occurred,
  104.  *        it means the value reached is the
  105.  *        maximum/minimum value representable.  It may
  106.  *        have come about due to an overflow.
  107.  *
  108.  *     INPUTS:
  109.  *
  110.  *       swIn               A Shortword input variable
  111.  *
  112.  *
  113.  *     OUTPUTS:             none
  114.  *
  115.  *     RETURN VALUE:        1 if input == 0x7fff or -0x8000
  116.  *                          0 otherwise
  117.  *
  118.  *     KEYWORDS: saturation, limit
  119.  *
  120.  ***************************************************************************/
  121. short  isSwLimit(Shortword swIn)
  122. {
  123.   Shortword swls;
  124.   short  siOut;
  125.   if (swIn != 0)
  126.   {
  127.     swls = shl(swIn, 1);
  128.     if (sub(swIn, swls) == 0)          /* logical compare outputs 1/0 */
  129.       siOut = 1;
  130.     else
  131.       siOut = 0;
  132.   }
  133.   else
  134.   {
  135.     siOut = 0;
  136.   }
  137.   return (siOut);
  138. }
  139. /****************************************************************************
  140.  *
  141.  *     FUNCTION NAME: L_mpy_ll
  142.  *
  143.  *     PURPOSE:    Multiply a 32 bit number (L_var1) and a 32 bit number
  144.  *                 (L_var2), and return a 32 bit result.
  145.  *
  146.  *     INPUTS:
  147.  *
  148.  *       L_var1             A Longword input variable
  149.  *
  150.  *       L_var2             A Longword input variable
  151.  *
  152.  *     OUTPUTS:             none
  153.  *
  154.  *     IMPLEMENTATION:
  155.  *
  156.  *        Performs a 31x31 bit multiply, Complexity=24 Ops.
  157.  *
  158.  *        Let x1x0, or y1y0, be the two constituent halves
  159.  *        of a 32 bit number.  This function performs the
  160.  *        following:
  161.  *
  162.  *        low = ((x0 >> 1)*(y0 >> 1)) >> 16     (low * low)
  163.  *        mid1 = [(x1 * (y0 >> 1)) >> 1 ]       (high * low)
  164.  *        mid2 = [(y1 * (x0 >> 1)) >> 1]        (high * low)
  165.  *        mid =  (mid1 + low + mid2) >> 14      (sum so far)
  166.  *        output = (y1*x1) + mid                (high * high)
  167.  *
  168.  *
  169.  *     RETURN VALUE:        A Longword value
  170.  *
  171.  *     KEYWORDS: mult,mpy,multiplication
  172.  *
  173.  ***************************************************************************/
  174. Longword L_mpy_ll(Longword L_var1, Longword L_var2)
  175. {
  176.   Shortword swLow1,
  177.          swLow2,
  178.          swHigh1,
  179.          swHigh2;
  180.   Longword L_varOut,
  181.          L_low,
  182.          L_mid1,
  183.          L_mid2,
  184.          L_mid;
  185.   swLow1 = shr(extract_l(L_var1), 1);
  186.   swLow1 = SW_MAX & swLow1;
  187.   swLow2 = shr(extract_l(L_var2), 1);
  188.   swLow2 = SW_MAX & swLow2;
  189.   swHigh1 = extract_h(L_var1);
  190.   swHigh2 = extract_h(L_var2);
  191.   L_low = L_mult(swLow1, swLow2);
  192.   L_low = L_shr(L_low, 16);
  193.   L_mid1 = L_mult(swHigh1, swLow2);
  194.   L_mid1 = L_shr(L_mid1, 1);
  195.   L_mid = L_add(L_mid1, L_low);
  196.   L_mid2 = L_mult(swHigh2, swLow1);
  197.   L_mid2 = L_shr(L_mid2, 1);
  198.   L_mid = L_add(L_mid, L_mid2);
  199.   L_mid = L_shr(L_mid, 14);
  200.   L_varOut = L_mac(L_mid, swHigh1, swHigh2);
  201.   return (L_varOut);
  202. }
  203. /****************************************************************************
  204.  *
  205.  *     FUNCTION NAME: L_mpy_ls
  206.  *
  207.  *     PURPOSE:    Multiply a 32 bit number (L_var2) and a 16 bit
  208.  *                 number (var1) returning a 32 bit result. L_var2
  209.  *                 is truncated to 31 bits prior to executing the
  210.  *                 multiply.
  211.  *
  212.  *     INPUTS:
  213.  *
  214.  *       L_var2             A Longword input variable
  215.  *
  216.  *       var1               A Shortword input variable
  217.  *
  218.  *     OUTPUTS:             none
  219.  *
  220.  *     RETURN VALUE:        A Longword value
  221.  *
  222.  *     KEYWORDS: mult,mpy,multiplication
  223.  *
  224.  ***************************************************************************/
  225. Longword L_mpy_ls(Longword L_var2, Shortword var1)
  226. {
  227.   Longword L_varOut;
  228.   Shortword swtemp;
  229.   swtemp = shr(extract_l(L_var2), 1);
  230.   swtemp = (short) 32767 & (short) swtemp;
  231.   L_varOut = L_mult(var1, swtemp);
  232.   L_varOut = L_shr(L_varOut, 15);
  233.   L_varOut = L_mac(L_varOut, var1, extract_h(L_var2));
  234.   return (L_varOut);
  235. }