capidtmf.c

linux-2.6.15.6

/*
 *
  Copyright (c) Eicon Networks, 2002.
 *
  This source file is supplied for the use with
  Eicon Networks range of DIVA Server Adapters.
 *
  Eicon File Revision :    2.1
 *
  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, or (at your option)
  any later version.
 *
  This program is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY OF ANY KIND WHATSOEVER INCLUDING ANY
  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., 675 Mass Ave, Cambridge, MA 02139, USA.
 *
 */

#include "platform.h"


  
  





#include "capidtmf.h"

/* #define TRACE_ */

#define FILE_ "CAPIDTMF.C"

/*---------------------------------------------------------------------------*/


#define trace(a)



/*---------------------------------------------------------------------------*/

static short capidtmf_expand_table_alaw[0x0100] =
{
   -5504,   5504,   -344,    344, -22016,  22016,  -1376,   1376,
   -2752,   2752,    -88,     88, -11008,  11008,   -688,    688,
   -7552,   7552,   -472,    472, -30208,  30208,  -1888,   1888,
   -3776,   3776,   -216,    216, -15104,  15104,   -944,    944,
   -4480,   4480,   -280,    280, -17920,  17920,  -1120,   1120,
   -2240,   2240,    -24,     24,  -8960,   8960,   -560,    560,
   -6528,   6528,   -408,    408, -26112,  26112,  -1632,   1632,
   -3264,   3264,   -152,    152, -13056,  13056,   -816,    816,
   -6016,   6016,   -376,    376, -24064,  24064,  -1504,   1504,
   -3008,   3008,   -120,    120, -12032,  12032,   -752,    752,
   -8064,   8064,   -504,    504, -32256,  32256,  -2016,   2016,
   -4032,   4032,   -248,    248, -16128,  16128,  -1008,   1008,
   -4992,   4992,   -312,    312, -19968,  19968,  -1248,   1248,
   -2496,   2496,    -56,     56,  -9984,   9984,   -624,    624,
   -7040,   7040,   -440,    440, -28160,  28160,  -1760,   1760,
   -3520,   3520,   -184,    184, -14080,  14080,   -880,    880,
   -5248,   5248,   -328,    328, -20992,  20992,  -1312,   1312,
   -2624,   2624,    -72,     72, -10496,  10496,   -656,    656,
   -7296,   7296,   -456,    456, -29184,  29184,  -1824,   1824,
   -3648,   3648,   -200,    200, -14592,  14592,   -912,    912,
   -4224,   4224,   -264,    264, -16896,  16896,  -1056,   1056,
   -2112,   2112,     -8,      8,  -8448,   8448,   -528,    528,
   -6272,   6272,   -392,    392, -25088,  25088,  -1568,   1568,
   -3136,   3136,   -136,    136, -12544,  12544,   -784,    784,
   -5760,   5760,   -360,    360, -23040,  23040,  -1440,   1440,
   -2880,   2880,   -104,    104, -11520,  11520,   -720,    720,
   -7808,   7808,   -488,    488, -31232,  31232,  -1952,   1952,
   -3904,   3904,   -232,    232, -15616,  15616,   -976,    976,
   -4736,   4736,   -296,    296, -18944,  18944,  -1184,   1184,
   -2368,   2368,    -40,     40,  -9472,   9472,   -592,    592,
   -6784,   6784,   -424,    424, -27136,  27136,  -1696,   1696,
   -3392,   3392,   -168,    168, -13568,  13568,   -848,    848
};

static short capidtmf_expand_table_ulaw[0x0100] =
{
  -32124,  32124,  -1884,   1884,  -7932,   7932,   -372,    372,
  -15996,  15996,   -876,    876,  -3900,   3900,   -120,    120,
  -23932,  23932,  -1372,   1372,  -5884,   5884,   -244,    244,
  -11900,  11900,   -620,    620,  -2876,   2876,    -56,     56,
  -28028,  28028,  -1628,   1628,  -6908,   6908,   -308,    308,
  -13948,  13948,   -748,    748,  -3388,   3388,    -88,     88,
  -19836,  19836,  -1116,   1116,  -4860,   4860,   -180,    180,
   -9852,   9852,   -492,    492,  -2364,   2364,    -24,     24,
  -30076,  30076,  -1756,   1756,  -7420,   7420,   -340,    340,
  -14972,  14972,   -812,    812,  -3644,   3644,   -104,    104,
  -21884,  21884,  -1244,   1244,  -5372,   5372,   -212,    212,
  -10876,  10876,   -556,    556,  -2620,   2620,    -40,     40,
  -25980,  25980,  -1500,   1500,  -6396,   6396,   -276,    276,
  -12924,  12924,   -684,    684,  -3132,   3132,    -72,     72,
  -17788,  17788,   -988,    988,  -4348,   4348,   -148,    148,
   -8828,   8828,   -428,    428,  -2108,   2108,     -8,      8,
  -31100,  31100,  -1820,   1820,  -7676,   7676,   -356,    356,
  -15484,  15484,   -844,    844,  -3772,   3772,   -112,    112,
  -22908,  22908,  -1308,   1308,  -5628,   5628,   -228,    228,
  -11388,  11388,   -588,    588,  -2748,   2748,    -48,     48,
  -27004,  27004,  -1564,   1564,  -6652,   6652,   -292,    292,
  -13436,  13436,   -716,    716,  -3260,   3260,    -80,     80,
  -18812,  18812,  -1052,   1052,  -4604,   4604,   -164,    164,
   -9340,   9340,   -460,    460,  -2236,   2236,    -16,     16,
  -29052,  29052,  -1692,   1692,  -7164,   7164,   -324,    324,
  -14460,  14460,   -780,    780,  -3516,   3516,    -96,     96,
  -20860,  20860,  -1180,   1180,  -5116,   5116,   -196,    196,
  -10364,  10364,   -524,    524,  -2492,   2492,    -32,     32,
  -24956,  24956,  -1436,   1436,  -6140,   6140,   -260,    260,
  -12412,  12412,   -652,    652,  -3004,   3004,    -64,     64,
  -16764,  16764,   -924,    924,  -4092,   4092,   -132,    132,
   -8316,   8316,   -396,    396,  -1980,   1980,      0,      0
};


/*---------------------------------------------------------------------------*/

static short capidtmf_recv_window_function[CAPIDTMF_RECV_ACCUMULATE_CYCLES] =
{
    -500L,   -999L,  -1499L,  -1998L,  -2496L,  -2994L,  -3491L,  -3988L,
   -4483L,  -4978L,  -5471L,  -5963L,  -6454L,  -6943L,  -7431L,  -7917L,
   -8401L,  -8883L,  -9363L,  -9840L, -10316L, -10789L, -11259L, -11727L,
  -12193L, -12655L, -13115L, -13571L, -14024L, -14474L, -14921L, -15364L,
  -15804L, -16240L, -16672L, -17100L, -17524L, -17944L, -18360L, -18772L,
  -19180L, -19583L, -19981L, -20375L, -20764L, -21148L, -21527L, -21901L,
  -22270L, -22634L, -22993L, -23346L, -23694L, -24037L, -24374L, -24705L,
  -25030L, -25350L, -25664L, -25971L, -26273L, -26568L, -26858L, -27141L,
  -27418L, -27688L, -27952L, -28210L, -28461L, -28705L, -28943L, -29174L,
  -29398L, -29615L, -29826L, -30029L, -30226L, -30415L, -30598L, -30773L,
  -30941L, -31102L, -31256L, -31402L, -31541L, -31673L, -31797L, -31914L,
  -32024L, -32126L, -32221L, -32308L, -32388L, -32460L, -32524L, -32581L,
  -32631L, -32673L, -32707L, -32734L, -32753L, -32764L, -32768L, -32764L,
  -32753L, -32734L, -32707L, -32673L, -32631L, -32581L, -32524L, -32460L,
  -32388L, -32308L, -32221L, -32126L, -32024L, -31914L, -31797L, -31673L,
  -31541L, -31402L, -31256L, -31102L, -30941L, -30773L, -30598L, -30415L,
  -30226L, -30029L, -29826L, -29615L, -29398L, -29174L, -28943L, -28705L,
  -28461L, -28210L, -27952L, -27688L, -27418L, -27141L, -26858L, -26568L,
  -26273L, -25971L, -25664L, -25350L, -25030L, -24705L, -24374L, -24037L,
  -23694L, -23346L, -22993L, -22634L, -22270L, -21901L, -21527L, -21148L,
  -20764L, -20375L, -19981L, -19583L, -19180L, -18772L, -18360L, -17944L,
  -17524L, -17100L, -16672L, -16240L, -15804L, -15364L, -14921L, -14474L,
  -14024L, -13571L, -13115L, -12655L, -12193L, -11727L, -11259L, -10789L,
  -10316L,  -9840L,  -9363L,  -8883L,  -8401L,  -7917L,  -7431L,  -6943L,
   -6454L,  -5963L,  -5471L,  -4978L,  -4483L,  -3988L,  -3491L,  -2994L,
   -2496L,  -1998L,  -1499L,   -999L,   -500L, 
};

static byte capidtmf_leading_zeroes_table[0x100] =
{
  8, 7, 6, 6, 5, 5, 5, 5, 4, 4, 4, 4, 4, 4, 4, 4,
  3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
};

#define capidtmf_byte_leading_zeroes(b)  (capidtmf_leading_zeroes_table[(BYTE)(b)])
#define capidtmf_word_leading_zeroes(w)  (((w) & 0xff00) ? capidtmf_leading_zeroes_table[(w) >> 8] : 8 + capidtmf_leading_zeroes_table[(w)])
#define capidtmf_dword_leading_zeroes(d)  (((d) & 0xffff0000L) ?    (((d) & 0xff000000L) ? capidtmf_leading_zeroes_table[(d) >> 24] : 8 + capidtmf_leading_zeroes_table[(d) >> 16]) :    (((d) & 0xff00) ? 16 + capidtmf_leading_zeroes_table[(d) >> 8] : 24 + capidtmf_leading_zeroes_table[(d)]))


/*---------------------------------------------------------------------------*/


static void capidtmf_goertzel_loop (long *buffer, long *coeffs, short *sample, long count)
{
  int i, j;
  long c, d, q0, q1, q2;

  for (i = 0; i < CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT - 1; i++)
  {
    q1 = buffer[i];
    q2 = buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT];
    d = coeffs[i] >> 1;
    c = d << 1;
    if (c >= 0)
    {
      for (j = 0; j < count; j++)
      {
        q0 = sample[j] - q2 + (c * (q1 >> 16)) + (((dword)(((dword) d) * ((dword)(q1 & 0xffff)))) >> 15);
        q2 = q1;
        q1 = q0;
      }
    }
    else
    {
      c = -c;
      d = -d;
      for (j = 0; j < count; j++)
      {
        q0 = sample[j] - q2 - ((c * (q1 >> 16)) + (((dword)(((dword) d) * ((dword)(q1 & 0xffff)))) >> 15));
        q2 = q1;
        q1 = q0;
      }
    }
    buffer[i] = q1;
    buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] = q2;
  }
  q1 = buffer[i];
  q2 = buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT];
  c = (coeffs[i] >> 1) << 1;
  if (c >= 0)
  {
    for (j = 0; j < count; j++)
    {
      q0 = sample[j] - q2 + (c * (q1 >> 16)) + (((dword)(((dword)(c >> 1)) * ((dword)(q1 & 0xffff)))) >> 15);
      q2 = q1;
      q1 = q0;
      c -= CAPIDTMF_RECV_FUNDAMENTAL_DECREMENT;
    }
  }
  else
  {
    c = -c;
    for (j = 0; j < count; j++)
    {
      q0 = sample[j] - q2 - ((c * (q1 >> 16)) + (((dword)(((dword)(c >> 1)) * ((dword)(q1 & 0xffff)))) >> 15));
      q2 = q1;
      q1 = q0;
      c += CAPIDTMF_RECV_FUNDAMENTAL_DECREMENT;
    }
  }
  coeffs[i] = c;
  buffer[i] = q1;
  buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] = q2;
}


static void capidtmf_goertzel_result (long *buffer, long *coeffs)
{
  int i;
  long d, e, q1, q2, lo, mid, hi;
  dword k;

  for (i = 0; i < CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT; i++)
  {
    q1 = buffer[i];
    q2 = buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT];
    d = coeffs[i] >> 1;
    if (d >= 0)
      d = ((d << 1) * (-q1 >> 16)) + (((dword)(((dword) d) * ((dword)(-q1 & 0xffff)))) >> 15);
    else
      d = ((-d << 1) * (-q1 >> 16)) + (((dword)(((dword) -d) * ((dword)(-q1 & 0xffff)))) >> 15);
    e = (q2 >= 0) ? q2 : -q2;
    if (d >= 0)
    {
      k = ((dword)(d & 0xffff)) * ((dword)(e & 0xffff));
      lo = k & 0xffff;
      mid = k >> 16;
      k = ((dword)(d >> 16)) * ((dword)(e & 0xffff));
      mid += k & 0xffff;
      hi = k >> 16;
      k = ((dword)(d & 0xffff)) * ((dword)(e >> 16));
      mid += k & 0xffff;
      hi += k >> 16;
      hi += ((dword)(d >> 16)) * ((dword)(e >> 16));
    }
    else
    {
      d = -d;
      k = ((dword)(d & 0xffff)) * ((dword)(e & 0xffff));
      lo = -((long)(k & 0xffff));
      mid = -((long)(k >> 16));
      k = ((dword)(d >> 16)) * ((dword)(e & 0xffff));
      mid -= k & 0xffff;
      hi = -((long)(k >> 16));
      k = ((dword)(d & 0xffff)) * ((dword)(e >> 16));
      mid -= k & 0xffff;
      hi -= k >> 16;
      hi -= ((dword)(d >> 16)) * ((dword)(e >> 16));
    }
    if (q2 < 0)
    {
      lo = -lo;
      mid = -mid;
      hi = -hi;
    }
    d = (q1 >= 0) ? q1 : -q1;
    k = ((dword)(d & 0xffff)) * ((dword)(d & 0xffff));
    lo += k & 0xffff;
    mid += k >> 16;
    k = ((dword)(d >> 16)) * ((dword)(d & 0xffff));
    mid += (k & 0xffff) << 1;
    hi += (k >> 16) << 1;
    hi += ((dword)(d >> 16)) * ((dword)(d >> 16));
    d = (q2 >= 0) ? q2 : -q2;
    k = ((dword)(d & 0xffff)) * ((dword)(d & 0xffff));
    lo += k & 0xffff;
    mid += k >> 16;
    k = ((dword)(d >> 16)) * ((dword)(d & 0xffff));
    mid += (k & 0xffff) << 1;
    hi += (k >> 16) << 1;
    hi += ((dword)(d >> 16)) * ((dword)(d >> 16));
    mid += lo >> 16;
    hi += mid >> 16;
    buffer[i] = (lo & 0xffff) | (mid << 16);
    buffer[i + CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] = hi;
  }
}


/*---------------------------------------------------------------------------*/

#define CAPIDTMF_RECV_GUARD_SNR_INDEX_697     0
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_770     1
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_852     2
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_941     3
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1209    4
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1336    5
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1477    6
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1633    7
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_635     8
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1010    9
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1140    10
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1272    11
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1405    12
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1555    13
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1715    14
#define CAPIDTMF_RECV_GUARD_SNR_INDEX_1875    15

#define CAPIDTMF_RECV_GUARD_SNR_DONTCARE      0xc000
#define CAPIDTMF_RECV_NO_DIGIT                0xff
#define CAPIDTMF_RECV_TIME_GRANULARITY        (CAPIDTMF_RECV_ACCUMULATE_CYCLES + 1)