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