filter.h
上传用户:wstnjxml
上传日期:2014-04-03
资源大小:7248k
文件大小:4k
源码类别:

Windows CE

开发平台:

C/C++

  1. /*
  2.  * filter.h
  3.  *
  4.  * Description:  TTAv1 filter functions
  5.  * Developed by: Alexander Djourik <ald@true-audio.com>
  6.  *               Pavel Zhilin <pzh@true-audio.com>
  7.  *
  8.  * Copyright (c) 2004 True Audio Software. All rights reserved.
  9.  *
  10.  */
  11. /*
  12.  * Redistribution and use in source and binary forms, with or without
  13.  * modification, are permitted provided that the following conditions
  14.  * are met:
  15.  *
  16.  * 1. Redistributions of source code must retain the above copyright
  17.  *    notice, this list of conditions and the following disclaimer.
  18.  * 2. Redistributions in binary form must reproduce the above copyright
  19.  *    notice, this list of conditions and the following disclaimer in the
  20.  *    documentation and/or other materials provided with the distribution.
  21.  * 3. Neither the name of the True Audio Software nor the names of its
  22.  *    contributors may be used to endorse or promote products derived
  23.  *    from this software without specific prior written permission.
  24.  *
  25.  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
  26.  * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
  27.  * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
  28.  * A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
  29.  * OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
  30.  * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
  31.  * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
  32.  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
  33.  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  34.  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
  35.  * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  36.  */
  37. #ifndef FILTER_H
  38. #define FILTER_H
  39. ///////// Filter Settings //////////
  40. static int flt_set[3] = {10, 9, 10};
  41. __inline void
  42. memshl (register int *pA, register int *pB) {
  43. *pA++ = *pB++;
  44. *pA++ = *pB++;
  45. *pA++ = *pB++;
  46. *pA++ = *pB++;
  47. *pA++ = *pB++;
  48. *pA++ = *pB++;
  49. *pA++ = *pB++;
  50. *pA   = *pB;
  51. }
  52. __inline void
  53. hybrid_filter (fltst *fs, int *in) {
  54. register int *pA = fs->dl;
  55. register int *pB = fs->qm;
  56. register int *pM = fs->dx;
  57. register int sum = fs->round;
  58. if (!fs->error) {
  59. sum += *pA++ * *pB, pB++;
  60. sum += *pA++ * *pB, pB++;
  61. sum += *pA++ * *pB, pB++;
  62. sum += *pA++ * *pB, pB++;
  63. sum += *pA++ * *pB, pB++;
  64. sum += *pA++ * *pB, pB++;
  65. sum += *pA++ * *pB, pB++;
  66. sum += *pA++ * *pB, pB++; pM += 8;
  67. } else if (fs->error < 0) {
  68. sum += *pA++ * (*pB -= *pM++), pB++;
  69. sum += *pA++ * (*pB -= *pM++), pB++;
  70. sum += *pA++ * (*pB -= *pM++), pB++;
  71. sum += *pA++ * (*pB -= *pM++), pB++;
  72. sum += *pA++ * (*pB -= *pM++), pB++;
  73. sum += *pA++ * (*pB -= *pM++), pB++;
  74. sum += *pA++ * (*pB -= *pM++), pB++;
  75. sum += *pA++ * (*pB -= *pM++), pB++;
  76. } else {
  77. sum += *pA++ * (*pB += *pM++), pB++;
  78. sum += *pA++ * (*pB += *pM++), pB++;
  79. sum += *pA++ * (*pB += *pM++), pB++;
  80. sum += *pA++ * (*pB += *pM++), pB++;
  81. sum += *pA++ * (*pB += *pM++), pB++;
  82. sum += *pA++ * (*pB += *pM++), pB++;
  83. sum += *pA++ * (*pB += *pM++), pB++;
  84. sum += *pA++ * (*pB += *pM++), pB++;
  85. }
  86. *(pM-0) = ((*(pA-1) >> 30) | 1) << 2;
  87. *(pM-1) = ((*(pA-2) >> 30) | 1) << 1;
  88. *(pM-2) = ((*(pA-3) >> 30) | 1) << 1;
  89. *(pM-3) = ((*(pA-4) >> 30) | 1);
  90. fs->error = *in;
  91. *in += (sum >> fs->shift);
  92. *pA = *in;
  93. *(pA-1) = *(pA-0) - *(pA-1);
  94. *(pA-2) = *(pA-1) - *(pA-2);
  95. *(pA-3) = *(pA-2) - *(pA-3);
  96. memshl (fs->dl, fs->dl + 1);
  97. memshl (fs->dx, fs->dx + 1);
  98. }
  99. void
  100. filter_init (fltst *fs, int shift) {
  101. memset (fs, 0, sizeof(fltst));
  102. fs->shift = shift;
  103. fs->round = 1 << (shift - 1);
  104. }
  105. #endif /* FILTER_H */