swh:1:snp:49cd9498d6cccc5e78252c27dcb645bcf7bf0c91
Raw File
Tip revision: 19f949f52599ba7c3f67a5897ac6be14bfcb1200 authored by Linus Torvalds on 18 February 2013, 23:58:34 UTC
Linux 3.8
Tip revision: 19f949f
capidtmf.c

/*
 *
 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);
}


/*---------------------------------------------------------------------------*/
back to top