gpclip.c revision e2e838afcf03e603a41a0455846eaf9614537c16
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: gpclip.c                                                     *
19e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*                                                                         *
20e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*      Description:To avoid unstable synthesis on frame erasure, the gain *
21e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*      need to be limited(gain pitch < 1.0) when the following            *
22e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*      case occurs                                                        *
23e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*      a resonance on LPC filter(lp_disp < 60Hz)                          *
24e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*      a good pitch prediction (lp_gp > 0.95)                             *
25e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard*                                                                         *
26e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard***************************************************************************/
27e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#include "typedef.h"
28e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#include "basic_op.h"
29e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
30e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#define DIST_ISF_MAX    307                /* 120 Hz (6400Hz=16384) */
31e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#define DIST_ISF_THRES  154                /* 60     (6400Hz=16384) */
32e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#define GAIN_PIT_THRES  14746              /* 0.9 in Q14 */
33e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#define GAIN_PIT_MIN    9830               /* 0.6 in Q14 */
34e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard#define M               16
35e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
36e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
37e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgardvoid Init_gp_clip(
38e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 mem[]                          /* (o) : memory of gain of pitch clipping algorithm */
39e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		)
40e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard{
41e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	mem[0] = DIST_ISF_MAX;
42e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	mem[1] = GAIN_PIT_MIN;
43e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard}
44e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
45e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
46e2e838afcf03e603a41a0455846eaf9614537c16Mans RullgardWord16 Gp_clip(
47e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 mem[]                          /* (i/o) : memory of gain of pitch clipping algorithm */
48e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	      )
49e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard{
50e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	Word16 clip = 0;
51e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	if ((mem[0] < DIST_ISF_THRES) && (mem[1] > GAIN_PIT_THRES))
52e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		clip = 1;
53e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
54e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	return (clip);
55e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard}
56e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
57e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
58e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgardvoid Gp_clip_test_isf(
59e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 isf[],                         /* (i)   : isf values (in frequency domain)           */
60e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 mem[]                          /* (i/o) : memory of gain of pitch clipping algorithm */
61e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		)
62e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard{
63e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	Word16 dist, dist_min;
64e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	Word32 i;
65e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
66e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	dist_min = vo_sub(isf[1], isf[0]);
67e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
68e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	for (i = 2; i < M - 1; i++)
69e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	{
70e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		dist = vo_sub(isf[i], isf[i - 1]);
71e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		if(dist < dist_min)
72e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		{
73e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard			dist_min = dist;
74e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		}
75e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	}
76e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
77e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	dist = extract_h(L_mac(vo_L_mult(26214, mem[0]), 6554, dist_min));
78e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
79e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	if (dist > DIST_ISF_MAX)
80e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	{
81e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		dist = DIST_ISF_MAX;
82e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	}
83e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	mem[0] = dist;
84e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
85e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	return;
86e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard}
87e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
88e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
89e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgardvoid Gp_clip_test_gain_pit(
90e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 gain_pit,                      /* (i) Q14 : gain of quantized pitch                    */
91e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		Word16 mem[]                          /* (i/o)   : memory of gain of pitch clipping algorithm */
92e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		)
93e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard{
94e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	Word16 gain;
95e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	Word32 L_tmp;
96e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	L_tmp = (29491 * mem[1])<<1;
97e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	L_tmp += (3277 * gain_pit)<<1;
98e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
99e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	gain = extract_h(L_tmp);
100e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
101e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	if(gain < GAIN_PIT_MIN)
102e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	{
103e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard		gain = GAIN_PIT_MIN;
104e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	}
105e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	mem[1] = gain;
106e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard	return;
107e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard}
108e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
109e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
110e2e838afcf03e603a41a0455846eaf9614537c16Mans Rullgard
111