text: 
  /*
   *
   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)
  
  #define CAPIDTMF_RECV_INDICATION_DIGIT        0x0001
  
  static long capidtmf_recv_goertzel_coef_table[CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] =
  {
  	0xda97L * 2,  /* 697 Hz (Low group 697 Hz) */
  	0xd299L * 2,  /* 770 Hz (Low group 770 Hz) */
  	0xc8cbL * 2,  /* 852 Hz (Low group 852 Hz) */
  	0xbd36L * 2,  /* 941 Hz (Low group 941 Hz) */
  	0x9501L * 2,  /* 1209 Hz (High group 1209 Hz) */
  	0x7f89L * 2,  /* 1336 Hz (High group 1336 Hz) */
  	0x6639L * 2,  /* 1477 Hz (High group 1477 Hz) */
  	0x48c6L * 2,  /* 1633 Hz (High group 1633 Hz) */
  	0xe14cL * 2,  /* 630 Hz (Lower guard of low group 631 Hz) */
  	0xb2e0L * 2,  /* 1015 Hz (Upper guard of low group 1039 Hz) */
  	0xa1a0L * 2,  /* 1130 Hz (Lower guard of high group 1140 Hz) */
  	0x8a87L * 2,  /* 1272 Hz (Guard between 1209 Hz and 1336 Hz: 1271 Hz) */
  	0x7353L * 2,  /* 1405 Hz (2nd harmonics of 697 Hz and guard between 1336 Hz and 1477 Hz: 1405 Hz) */
  	0x583bL * 2,  /* 1552 Hz (2nd harmonics of 770 Hz and guard between 1477 Hz and 1633 Hz: 1553 Hz) */
  	0x37d8L * 2,  /* 1720 Hz (2nd harmonics of 852 Hz and upper guard of high group: 1715 Hz) */
  	0x0000L * 2   /* 100-630 Hz (fundamentals) */
  };
  
  
  static word capidtmf_recv_guard_snr_low_table[CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] =
  {
  	14,                                    /* Low group peak versus 697 Hz */
  	14,                                    /* Low group peak versus 770 Hz */
  	16,                                    /* Low group peak versus 852 Hz */
  	16,                                    /* Low group peak versus 941 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1209 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1336 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1477 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1633 Hz */
  	14,                                    /* Low group peak versus 635 Hz */
  	16,                                    /* Low group peak versus 1010 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1140 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* Low group peak versus 1272 Hz */
  	DSPDTMF_RX_HARMONICS_SEL_DEFAULT - 8,  /* Low group peak versus 1405 Hz */
  	DSPDTMF_RX_HARMONICS_SEL_DEFAULT - 4,  /* Low group peak versus 1555 Hz */
  	DSPDTMF_RX_HARMONICS_SEL_DEFAULT - 4,  /* Low group peak versus 1715 Hz */
  	12                                     /* Low group peak versus 100-630 Hz */
  };
  
  
  static word capidtmf_recv_guard_snr_high_table[CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT] =
  {
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 697 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 770 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 852 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 941 Hz */
  	20,                                    /* High group peak versus 1209 Hz */
  	20,                                    /* High group peak versus 1336 Hz */
  	20,                                    /* High group peak versus 1477 Hz */
  	20,                                    /* High group peak versus 1633 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 635 Hz */
  	CAPIDTMF_RECV_GUARD_SNR_DONTCARE,      /* High group peak versus 1010 Hz */
  	16,                                    /* High group peak versus 1140 Hz */
  	4,                                     /* High group peak versus 1272 Hz */
  	6,                                     /* High group peak versus 1405 Hz */
  	8,                                     /* High group peak versus 1555 Hz */
  	16,                                    /* High group peak versus 1715 Hz */
  	12                                     /* High group peak versus 100-630 Hz */
  };
  
  
  /*---------------------------------------------------------------------------*/
  
  static void capidtmf_recv_init(t_capidtmf_state *p_state)
  {
  	p_state->recv.min_gap_duration = 1;
  	p_state->recv.min_digit_duration = 1;
  
  	p_state->recv.cycle_counter = 0;
  	p_state->recv.current_digit_on_time = 0;
  	p_state->recv.current_digit_off_time = 0;
  	p_state->recv.current_digit_value = CAPIDTMF_RECV_NO_DIGIT;
  
  	p_state->recv.digit_write_pos = 0;
  	p_state->recv.digit_read_pos = 0;
  	p_state->recv.indication_state = 0;
  	p_state->recv.indication_state_ack = 0;
  	p_state->recv.state = CAPIDTMF_RECV_STATE_IDLE;
  }
  
  
  void capidtmf_recv_enable(t_capidtmf_state *p_state, word min_digit_duration, word min_gap_duration)
  {
  	p_state->recv.indication_state_ack &= CAPIDTMF_RECV_INDICATION_DIGIT;
  	p_state->recv.min_digit_duration = (word)(((((dword) min_digit_duration) * 8) +
  						   ((dword)(CAPIDTMF_RECV_TIME_GRANULARITY / 2))) / ((dword) CAPIDTMF_RECV_TIME_GRANULARITY));
  	if (p_state->recv.min_digit_duration <= 1)
  		p_state->recv.min_digit_duration = 1;
  	else
  		(p_state->recv.min_digit_duration)--;
  	p_state->recv.min_gap_duration =
  		(word)((((dword) min_gap_duration) * 8) / ((dword) CAPIDTMF_RECV_TIME_GRANULARITY));
  	if (p_state->recv.min_gap_duration <= 1)
  		p_state->recv.min_gap_duration = 1;
  	else
  		(p_state->recv.min_gap_duration)--;
  	p_state->recv.state |= CAPIDTMF_RECV_STATE_DTMF_ACTIVE;
  }
  
  
  void capidtmf_recv_disable(t_capidtmf_state *p_state)
  {
  	p_state->recv.state &= ~CAPIDTMF_RECV_STATE_DTMF_ACTIVE;
  	if (p_state->recv.state == CAPIDTMF_RECV_STATE_IDLE)
  		capidtmf_recv_init(p_state);
  	else
  	{
  		p_state->recv.cycle_counter = 0;
  		p_state->recv.current_digit_on_time = 0;
  		p_state->recv.current_digit_off_time = 0;
  		p_state->recv.current_digit_value = CAPIDTMF_RECV_NO_DIGIT;
  	}
  }
  
  
  word capidtmf_recv_indication(t_capidtmf_state *p_state, byte *buffer)
  {
  	word i, j, k, flags;
  
  	flags = p_state->recv.indication_state ^ p_state->recv.indication_state_ack;
  	p_state->recv.indication_state_ack ^= flags & CAPIDTMF_RECV_INDICATION_DIGIT;
  	if (p_state->recv.digit_write_pos != p_state->recv.digit_read_pos)
  	{
  		i = 0;
  		k = p_state->recv.digit_write_pos;
  		j = p_state->recv.digit_read_pos;
  		do
  		{
  			buffer[i++] = p_state->recv.digit_buffer[j];
  			j = (j == CAPIDTMF_RECV_DIGIT_BUFFER_SIZE - 1) ? 0 : j + 1;
  		} while (j != k);
  		p_state->recv.digit_read_pos = k;
  		return (i);
  	}
  	p_state->recv.indication_state_ack ^= flags;
  	return (0);
  }
  
  
  #define CAPIDTMF_RECV_WINDOWED_SAMPLES  32
  
  void capidtmf_recv_block(t_capidtmf_state *p_state, byte *buffer, word length)
  {
  	byte result_digit;
  	word sample_number, cycle_counter, n, i;
  	word low_peak, high_peak;
  	dword lo, hi;
  	byte   *p;
  	short *q;
  	byte goertzel_result_buffer[CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT];
  	short windowed_sample_buffer[CAPIDTMF_RECV_WINDOWED_SAMPLES];
  
  
  	if (p_state->recv.state & CAPIDTMF_RECV_STATE_DTMF_ACTIVE)
  	{
  		cycle_counter = p_state->recv.cycle_counter;
  		sample_number = 0;
  		while (sample_number < length)
  		{
  			if (cycle_counter < CAPIDTMF_RECV_ACCUMULATE_CYCLES)
  			{
  				if (cycle_counter == 0)
  				{
  					for (i = 0; i < CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT; i++)
  					{
  						p_state->recv.goertzel_buffer[0][i] = 0;
  						p_state->recv.goertzel_buffer[1][i] = 0;
  					}
  				}
  				n = CAPIDTMF_RECV_ACCUMULATE_CYCLES - cycle_counter;
  				if (n > length - sample_number)
  					n = length - sample_number;
  				if (n > CAPIDTMF_RECV_WINDOWED_SAMPLES)
  					n = CAPIDTMF_RECV_WINDOWED_SAMPLES;
  				p = buffer + sample_number;
  				q = capidtmf_recv_window_function + cycle_counter;
  				if (p_state->ulaw)
  				{
  					for (i = 0; i < n; i++)
  					{
  						windowed_sample_buffer[i] =
  							(short)((capidtmf_expand_table_ulaw[p[i]] * ((long)(q[i]))) >> 15);
  					}
  				}
  				else
  				{
  					for (i = 0; i < n; i++)
  					{
  						windowed_sample_buffer[i] =
  							(short)((capidtmf_expand_table_alaw[p[i]] * ((long)(q[i]))) >> 15);
  					}
  				}
  				capidtmf_recv_goertzel_coef_table[CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT - 1] = CAPIDTMF_RECV_FUNDAMENTAL_OFFSET;
  				capidtmf_goertzel_loop(p_state->recv.goertzel_buffer[0],
  						       capidtmf_recv_goertzel_coef_table, windowed_sample_buffer, n);
  				cycle_counter += n;
  				sample_number += n;
  			}
  			else
  			{
  				capidtmf_goertzel_result(p_state->recv.goertzel_buffer[0],
  							 capidtmf_recv_goertzel_coef_table);
  				for (i = 0; i < CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT; i++)
  				{
  					lo = (dword)(p_state->recv.goertzel_buffer[0][i]);
  					hi = (dword)(p_state->recv.goertzel_buffer[1][i]);
  					if (hi != 0)
  					{
  						n = capidtmf_dword_leading_zeroes(hi);
  						hi = (hi << n) | (lo >> (32 - n));
  					}
  					else
  					{
  						n = capidtmf_dword_leading_zeroes(lo);
  						hi = lo << n;
  						n += 32;
  					}
  					n = 195 - 3 * n;
  					if (hi >= 0xcb300000L)
  						n += 2;
  					else if (hi >= 0xa1450000L)
  						n++;
  					goertzel_result_buffer[i] = (byte) n;
  				}
  				low_peak = DSPDTMF_RX_SENSITIVITY_LOW_DEFAULT;
  				result_digit = CAPIDTMF_RECV_NO_DIGIT;
  				for (i = 0; i < CAPIDTMF_LOW_GROUP_FREQUENCIES; i++)
  				{
  					if (goertzel_result_buffer[i] > low_peak)
  					{
  						low_peak = goertzel_result_buffer[i];
  						result_digit = (byte) i;
  					}
  				}
  				high_peak = DSPDTMF_RX_SENSITIVITY_HIGH_DEFAULT;
  				n = CAPIDTMF_RECV_NO_DIGIT;
  				for (i = CAPIDTMF_LOW_GROUP_FREQUENCIES; i < CAPIDTMF_RECV_BASE_FREQUENCY_COUNT; i++)
  				{
  					if (goertzel_result_buffer[i] > high_peak)
  					{
  						high_peak = goertzel_result_buffer[i];
  						n = (i - CAPIDTMF_LOW_GROUP_FREQUENCIES) << 2;
  					}
  				}
  				result_digit |= (byte) n;
  				if (low_peak + DSPDTMF_RX_HIGH_EXCEEDING_LOW_DEFAULT < high_peak)
  					result_digit = CAPIDTMF_RECV_NO_DIGIT;
  				if (high_peak + DSPDTMF_RX_LOW_EXCEEDING_HIGH_DEFAULT < low_peak)
  					result_digit = CAPIDTMF_RECV_NO_DIGIT;
  				n = 0;
  				for (i = 0; i < CAPIDTMF_RECV_TOTAL_FREQUENCY_COUNT; i++)
  				{
  					if ((((short)(low_peak - goertzel_result_buffer[i] - capidtmf_recv_guard_snr_low_table[i])) < 0)
  					    || (((short)(high_peak - goertzel_result_buffer[i] - capidtmf_recv_guard_snr_high_table[i])) < 0))
  					{
  						n++;
  					}
  				}
  				if (n != 2)
  					result_digit = CAPIDTMF_RECV_NO_DIGIT;
  
  				if (result_digit == CAPIDTMF_RECV_NO_DIGIT)
  				{
  					if (p_state->recv.current_digit_on_time != 0)
  					{
  						if (++(p_state->recv.current_digit_off_time) >= p_state->recv.min_gap_duration)
  						{
  							p_state->recv.current_digit_on_time = 0;
  							p_state->recv.current_digit_off_time = 0;
  						}
  					}
  					else
  					{
  						if (p_state->recv.current_digit_off_time != 0)
  							(p_state->recv.current_digit_off_time)--;
  					}
  				}
  				else
  				{
  					if ((p_state->recv.current_digit_on_time == 0)
  					    && (p_state->recv.current_digit_off_time != 0))
  					{
  						(p_state->recv.current_digit_off_time)--;
  					}
  					else
  					{
  						n = p_state->recv.current_digit_off_time;
  						if ((p_state->recv.current_digit_on_time != 0)
  						    && (result_digit != p_state->recv.current_digit_value))
  						{
  							p_state->recv.current_digit_on_time = 0;
  							n = 0;
  						}
  						p_state->recv.current_digit_value = result_digit;
  						p_state->recv.current_digit_off_time = 0;
  						if (p_state->recv.current_digit_on_time != 0xffff)
  						{
  							p_state->recv.current_digit_on_time += n + 1;
  							if (p_state->recv.current_digit_on_time >= p_state->recv.min_digit_duration)
  							{
  								p_state->recv.current_digit_on_time = 0xffff;
  								i = (p_state->recv.digit_write_pos == CAPIDTMF_RECV_DIGIT_BUFFER_SIZE - 1) ?
  									0 : p_state->recv.digit_write_pos + 1;
  								if (i == p_state->recv.digit_read_pos)
  								{
  									trace(dprintf("%s,%d: Receive digit overrun",
  										      (char *)(FILE_), __LINE__));
  								}
  								else
  								{
  									p_state->recv.digit_buffer[p_state->recv.digit_write_pos] = result_digit;
  									p_state->recv.digit_write_pos = i;
  									p_state->recv.indication_state =
  										(p_state->recv.indication_state & ~CAPIDTMF_RECV_INDICATION_DIGIT) |
  										(~p_state->recv.indication_state_ack & CAPIDTMF_RECV_INDICATION_DIGIT);
  								}
  							}
  						}
  					}
  				}
  				cycle_counter = 0;
  				sample_number++;
  			}
  		}
  		p_state->recv.cycle_counter = cycle_counter;
  	}
  }
  
  
  void capidtmf_init(t_capidtmf_state *p_state, byte ulaw)
  {
  	p_state->ulaw = ulaw;
  	capidtmf_recv_init(p_state);
  }
  
  
  /*---------------------------------------------------------------------------*/
  