1e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard/* 2e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** Copyright 2003-2010, VisualOn, Inc. 3e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** 4e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** Licensed under the Apache License, Version 2.0 (the "License"); 5e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** you may not use this file except in compliance with the License. 6e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** You may obtain a copy of the License at 7e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** 8e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** http://www.apache.org/licenses/LICENSE-2.0 9e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** 10e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** Unless required by applicable law or agreed to in writing, software 11e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** distributed under the License is distributed on an "AS IS" BASIS, 12e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 13e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** See the License for the specific language governing permissions and 14e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard ** limitations under the License. 15e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard */ 16e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 17e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard/*********************************************************************** 18e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard File: convolve.c 19e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 205d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Description:Perform the convolution between two vectors x[] and h[] 215d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen and write the result in the vector y[] 22e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 23e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard************************************************************************/ 24e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 25e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#include "typedef.h" 26e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#include "basic_op.h" 27e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 2884333e0475bc911adc16417f4ca327c975cf6c36Andreas Huber#define UNUSED(x) (void)(x) 2984333e0475bc911adc16417f4ca327c975cf6c36Andreas Huber 30e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgardvoid Convolve ( 315d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word16 x[], /* (i) : input vector */ 325d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word16 h[], /* (i) : impulse response */ 335d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word16 y[], /* (o) : output vector */ 345d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word16 L /* (i) : vector size */ 355d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen ) 36e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard{ 375d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word32 i, n; 385d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word16 *tmpH,*tmpX; 395d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen Word32 s; 4084333e0475bc911adc16417f4ca327c975cf6c36Andreas Huber UNUSED(L); 4184333e0475bc911adc16417f4ca327c975cf6c36Andreas Huber 425d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen for (n = 0; n < 64;) 435d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen { 445d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpH = h+n; 455d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpX = x; 465d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i=n+1; 475d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen s = vo_mult32((*tmpX++), (*tmpH--));i--; 485d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen while(i>0) 495d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen { 505d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 515d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 525d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 535d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 545d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i -= 4; 555d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen } 565d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen y[n] = voround(L_shl(s, 1)); 575d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen n++; 58e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 595d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpH = h+n; 605d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpX = x; 615d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i=n+1; 625d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = vo_mult32((*tmpX++), (*tmpH--)); 635d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen i--; 645d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 655d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen i--; 66e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 675d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen while(i>0) 685d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen { 695d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 705d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 715d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 725d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 735d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i -= 4; 745d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen } 755d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen y[n] = voround(L_shl(s, 1)); 765d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen n++; 77e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 785d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpH = h+n; 795d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpX = x; 805d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i=n+1; 815d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = vo_mult32((*tmpX++), (*tmpH--)); 825d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen i--; 835d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 845d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen i--; 855d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 865d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen i--; 87e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 885d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen while(i>0) 895d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen { 905d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 915d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 925d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 935d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 945d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i -= 4; 955d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen } 965d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen y[n] = voround(L_shl(s, 1)); 975d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen n++; 98e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 995d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen s = 0; 1005d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpH = h+n; 1015d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen tmpX = x; 1025d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i=n+1; 1035d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen while(i>0) 1045d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen { 1055d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 1065d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 1075d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 1085d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen s = L_add(s, vo_mult32((*tmpX++), (*tmpH--))); 1095d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen i -= 4; 1105d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen } 1115d453222ae6dcc10efedb1e4805247d7c22a4168Marco Nelissen y[n] = voround(L_shl(s, 1)); 1125d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen n++; 1135d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen } 1145d5c3a132bb446ac78a37dfaac24a46cacf0dd73Marco Nelissen return; 115e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard} 116e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 117e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 118e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard 119