/src/speex/libspeex/vbr.c

Source
/* Copyright (C) 2002 Jean-Marc Valin
   File: vbr.c

   VBR-related routines

   Redistribution and use in source and binary forms, with or without
   modification, are permitted provided that the following conditions
   are met:

   - Redistributions of source code must retain the above copyright
   notice, this list of conditions and the following disclaimer.

   - Redistributions in binary form must reproduce the above copyright
   notice, this list of conditions and the following disclaimer in the
   documentation and/or other materials provided with the distribution.

   - Neither the name of the Xiph.org Foundation nor the names of its
   contributors may be used to endorse or promote products derived from
   this software without specific prior written permission.

   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
   ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
   A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR
   CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
   EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
   PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
   PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
   LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
   NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
   SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

*/

#ifdef HAVE_CONFIG_H
#include "config.h"
#endif

#include "vbr.h"
#include <math.h>


#define sqr(x) ((x)*(x))

#define MIN_ENERGY 6000
#define NOISE_POW .3

#ifndef DISABLE_VBR

const float vbr_nb_thresh[9][11]={
   {-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /*   CNG   */
   { 4.0f,  2.5f,  2.0f,  1.2f,  0.5f,-0.25f, -0.5f, -0.7f, -0.8f, -0.9f, -1.0f}, /*  2 kbps */
   {10.0f,  6.5f,  5.2f,  4.5f,  3.9f, 3.7f,  3.0f,  2.5f,  2.3f,  1.8f,  1.0f}, /*  6 kbps */
   {11.0f,  8.8f,  7.5f,  6.5f,  5.0f,  4.2f,  3.9f,  3.9f,  3.5f,  3.0f,  1.0f}, /*  8 kbps */
   {11.0f, 11.0f,  9.9f,  8.5f,  7.0f, 5.25f,  4.5f,  4.0f,  4.0f,  4.0f,  2.0f}, /* 11 kbps */
   {11.0f, 11.0f, 11.0f, 11.0f,  9.5f, 9.25f,  8.0f,  7.0f,  5.0f,  4.0f,  3.0f}, /* 15 kbps */
   {11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f,  9.5f,  8.5f,  6.2f,  5.2f,  5.0f}, /* 18 kbps */
   {11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 10.0f,  9.8f,  7.5f}, /* 24 kbps */
   { 7.0f,  4.5f,  3.7f,  3.0f,  2.5f,  1.0f,  1.8f,  1.5f,  1.0f,  0.0f,  0.0f}  /*  4 kbps */
};


const float vbr_hb_thresh[5][11]={
   {-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* silence */
   {-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /*  2 kbps */
   {11.0f, 11.0f,  9.5f,  8.5f,  7.5f,  6.0f,  5.0f,  3.9f,  3.0f,  2.0f,  1.0f}, /*  6 kbps */
   {11.0f, 11.0f, 11.0f, 11.0f, 11.0f,  9.5f,  8.7f,  7.8f,  7.0f,  6.5f,  4.0f}, /* 10 kbps */
   {11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f,  9.8f,  7.5f,  5.5f}  /* 18 kbps */
};

const float vbr_uhb_thresh[2][11]={
   {-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* silence */
   { 3.9f,  2.5f,  0.0f,  0.0f,  0.0f,  0.0f,  0.0f,  0.0f,  0.0f,  0.0f, -1.0f}  /*  2 kbps */
};

void vbr_init(VBRState *vbr)
{
   int i;

   vbr->average_energy=1600000;
   vbr->last_energy=1;
   vbr->accum_sum=0;
   vbr->soft_pitch=0;
   vbr->last_pitch_coef=0;
   vbr->last_quality=0;

   vbr->noise_accum = .05*pow(MIN_ENERGY, NOISE_POW);
   vbr->noise_accum_count=.05;
   vbr->noise_level=vbr->noise_accum/vbr->noise_accum_count;
   vbr->consec_noise=0;


   for (i=0;i<VBR_MEMORY_SIZE;i++)
      vbr->last_log_energy[i] = log(MIN_ENERGY);
}


/*
  This function should analyse the signal and decide how critical the
  coding error will be perceptually. The following factors should be
  taken into account:

  -Attacks (positive energy derivative) should be coded with more bits

  -Stationary voiced segments should receive more bits

  -Segments with (very) low absolute energy should receive less bits (maybe
  only shaped noise?)

  -DTX for near-zero energy?

  -Stationary fricative segments should have less bits

  -Temporal masking: when energy slope is decreasing, decrease the bit-rate

  -Decrease bit-rate for males (low pitch)?

  -(wideband only) less bits in the high-band when signal is very
  non-stationary (harder to notice high-frequency noise)???

*/

float vbr_analysis(VBRState *vbr, spx_word16_t *sig, int len, int pitch, float pitch_coef)
{
   int i;
   float ener=0, ener1=0, ener2=0;
   float qual=7;
   float log_energy;
   float non_st=0;
   float voicing;
   float pow_ener;

   for (i=0;i<len>>1;i++)
      ener1 += ((float)sig[i])*sig[i];

   for (i=len>>1;i<len;i++)
      ener2 += ((float)sig[i])*sig[i];
   ener=ener1+ener2;

   log_energy = log(ener+MIN_ENERGY);
   for (i=0;i<VBR_MEMORY_SIZE;i++)
      non_st += sqr(log_energy-vbr->last_log_energy[i]);
   non_st =  non_st/(30*VBR_MEMORY_SIZE);
   if (non_st>1)
      non_st=1;

   voicing = 3*(pitch_coef-.4)*fabs(pitch_coef-.4);
   vbr->average_energy = 0.9*vbr->average_energy + .1*ener;
   vbr->noise_level=vbr->noise_accum/vbr->noise_accum_count;
   pow_ener = pow(ener,NOISE_POW);
   if (vbr->noise_accum_count<.06 && ener>MIN_ENERGY)
      vbr->noise_accum = .05*pow_ener;

   if ((voicing<.3 && non_st < .2 && pow_ener < 1.2*vbr->noise_level)
       || (voicing<.3 && non_st < .05 && pow_ener < 1.5*vbr->noise_level)
       || (voicing<.4 && non_st < .05 && pow_ener < 1.2*vbr->noise_level)
       || (voicing<0 && non_st < .05))
   {
      float tmp;

      vbr->consec_noise++;
      if (pow_ener > 3*vbr->noise_level)
         tmp = 3*vbr->noise_level;
      else
         tmp = pow_ener;
      if (vbr->consec_noise>=4)
      {
         vbr->noise_accum = .95*vbr->noise_accum + .05*tmp;
         vbr->noise_accum_count = .95*vbr->noise_accum_count + .05;
      }
   } else {
      vbr->consec_noise=0;
   }

   if (pow_ener < vbr->noise_level && ener>MIN_ENERGY)
   {
      vbr->noise_accum = .95*vbr->noise_accum + .05*pow_ener;
      vbr->noise_accum_count = .95*vbr->noise_accum_count + .05;
   }

   /* Checking for very low absolute energy */
   if (ener < 30000)
   {
      qual -= .7;
      if (ener < 10000)
         qual-=.7;
      if (ener < 3000)
         qual-=.7;
   } else {
      float short_diff, long_diff;
      short_diff = log((ener+1)/(1+vbr->last_energy));
      long_diff = log((ener+1)/(1+vbr->average_energy));
      /*fprintf (stderr, "%f %f\n", short_diff, long_diff);*/

      if (long_diff<-5)
         long_diff=-5;
      if (long_diff>2)
         long_diff=2;

      if (long_diff>0)
         qual += .6*long_diff;
      if (long_diff<0)
         qual += .5*long_diff;
      if (short_diff>0)
      {
         if (short_diff>5)
            short_diff=5;
         qual += 1*short_diff;
      }
      /* Checking for energy increases */
      if (ener2 > 1.6*ener1)
         qual += .5;
   }
   vbr->last_energy = ener;
   vbr->soft_pitch = .8*vbr->soft_pitch + .2*pitch_coef;
   qual += 2.2*((pitch_coef-.4) + (vbr->soft_pitch-.4));

   if (qual < vbr->last_quality)
      qual = .5*qual + .5*vbr->last_quality;
   if (qual<4)
      qual=4;
   if (qual>10)
      qual=10;

   /*
   if (vbr->consec_noise>=2)
      qual-=1.3;
   if (vbr->consec_noise>=5)
      qual-=1.3;
   if (vbr->consec_noise>=12)
      qual-=1.3;
   */
   if (vbr->consec_noise>=3)
      qual=4;

   if (vbr->consec_noise)
      qual -= 1.0 * (log(3.0 + vbr->consec_noise)-log(3));
   if (qual<0)
      qual=0;

   if (ener<1600000)
   {
      if (vbr->consec_noise>2)
         qual-=0.5*(log(3.0 + vbr->consec_noise)-log(3));
      if (ener<10000&&vbr->consec_noise>2)
         qual-=0.5*(log(3.0 + vbr->consec_noise)-log(3));
      if (qual<0)
         qual=0;
      qual += .3*log(.0001+ener/1600000.0);
   }
   if (qual<-1)
      qual=-1;

   /*printf ("%f %f %f %f\n", qual, voicing, non_st, pow_ener/(.01+vbr->noise_level));*/

   vbr->last_pitch_coef = pitch_coef;
   vbr->last_quality = qual;

   for (i=VBR_MEMORY_SIZE-1;i>0;i--)
      vbr->last_log_energy[i] = vbr->last_log_energy[i-1];
   vbr->last_log_energy[0] = log_energy;

   /*printf ("VBR: %f %f %f %f\n", (float)(log_energy-log(vbr->average_energy+MIN_ENERGY)), non_st, voicing, vbr->noise_level);*/

   return qual;
}

void vbr_destroy(VBRState *vbr)
{
}

#endif /* #ifndef DISABLE_VBR */

Coverage Report

Created: 2026-10-17 07:18

Line	Count	Source
1		/* Copyright (C) 2002 Jean-Marc Valin
2		File: vbr.c
3
4		VBR-related routines
5
6		Redistribution and use in source and binary forms, with or without
7		modification, are permitted provided that the following conditions
8		are met:
9
10		- Redistributions of source code must retain the above copyright
11		notice, this list of conditions and the following disclaimer.
12
13		- Redistributions in binary form must reproduce the above copyright
14		notice, this list of conditions and the following disclaimer in the
15		documentation and/or other materials provided with the distribution.
16
17		- Neither the name of the Xiph.org Foundation nor the names of its
18		contributors may be used to endorse or promote products derived from
19		this software without specific prior written permission.
20
21		THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
22		``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
23		LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
24		A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR
25		CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
26		EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
27		PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
28		PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
29		LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
30		NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
31		SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
32
33		*/
34
35		#ifdef HAVE_CONFIG_H
36		#include "config.h"
37		#endif
38
39		#include "vbr.h"
40		#include <math.h>
41
42
43	72.1k	#define sqr(x) ((x)*(x))
44
45	45.2k	#define MIN_ENERGY 6000
46	16.9k	#define NOISE_POW .3
47
48		#ifndef DISABLE_VBR
49
50		const float vbr_nb_thresh[9][11]={
51		{-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* CNG */
52		{ 4.0f, 2.5f, 2.0f, 1.2f, 0.5f,-0.25f, -0.5f, -0.7f, -0.8f, -0.9f, -1.0f}, /* 2 kbps */
53		{10.0f, 6.5f, 5.2f, 4.5f, 3.9f, 3.7f, 3.0f, 2.5f, 2.3f, 1.8f, 1.0f}, /* 6 kbps */
54		{11.0f, 8.8f, 7.5f, 6.5f, 5.0f, 4.2f, 3.9f, 3.9f, 3.5f, 3.0f, 1.0f}, /* 8 kbps */
55		{11.0f, 11.0f, 9.9f, 8.5f, 7.0f, 5.25f, 4.5f, 4.0f, 4.0f, 4.0f, 2.0f}, /* 11 kbps */
56		{11.0f, 11.0f, 11.0f, 11.0f, 9.5f, 9.25f, 8.0f, 7.0f, 5.0f, 4.0f, 3.0f}, /* 15 kbps */
57		{11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 9.5f, 8.5f, 6.2f, 5.2f, 5.0f}, /* 18 kbps */
58		{11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 10.0f, 9.8f, 7.5f}, /* 24 kbps */
59		{ 7.0f, 4.5f, 3.7f, 3.0f, 2.5f, 1.0f, 1.8f, 1.5f, 1.0f, 0.0f, 0.0f} /* 4 kbps */
60		};
61
62
63		const float vbr_hb_thresh[5][11]={
64		{-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* silence */
65		{-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* 2 kbps */
66		{11.0f, 11.0f, 9.5f, 8.5f, 7.5f, 6.0f, 5.0f, 3.9f, 3.0f, 2.0f, 1.0f}, /* 6 kbps */
67		{11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 9.5f, 8.7f, 7.8f, 7.0f, 6.5f, 4.0f}, /* 10 kbps */
68		{11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 11.0f, 9.8f, 7.5f, 5.5f} /* 18 kbps */
69		};
70
71		const float vbr_uhb_thresh[2][11]={
72		{-1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f, -1.0f}, /* silence */
73		{ 3.9f, 2.5f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, -1.0f} /* 2 kbps */
74		};
75
76		void vbr_init(VBRState *vbr)
77	2.51k	{
78	2.51k	int i;
79
80	2.51k	vbr->average_energy=1600000;
81	2.51k	vbr->last_energy=1;
82	2.51k	vbr->accum_sum=0;
83	2.51k	vbr->soft_pitch=0;
84	2.51k	vbr->last_pitch_coef=0;
85	2.51k	vbr->last_quality=0;
86
87	2.51k	vbr->noise_accum = .05*pow(MIN_ENERGY, NOISE_POW);
88	2.51k	vbr->noise_accum_count=.05;
89	2.51k	vbr->noise_level=vbr->noise_accum/vbr->noise_accum_count;
90	2.51k	vbr->consec_noise=0;
91
92
93	15.1k	for (i=0;i<VBR_MEMORY_SIZE;i++)
94	12.5k	vbr->last_log_energy[i] = log(MIN_ENERGY);
95	2.51k	}
96
97
98		/*
99		This function should analyse the signal and decide how critical the
100		coding error will be perceptually. The following factors should be
101		taken into account:
102
103		-Attacks (positive energy derivative) should be coded with more bits
104
105		-Stationary voiced segments should receive more bits
106
107		-Segments with (very) low absolute energy should receive less bits (maybe
108		only shaped noise?)
109
110		-DTX for near-zero energy?
111
112		-Stationary fricative segments should have less bits
113
114		-Temporal masking: when energy slope is decreasing, decrease the bit-rate
115
116		-Decrease bit-rate for males (low pitch)?
117
118		-(wideband only) less bits in the high-band when signal is very
119		non-stationary (harder to notice high-frequency noise)???
120
121		*/
122
123		float vbr_analysis(VBRState vbr, spx_word16_t sig, int len, int pitch, float pitch_coef)
124	14.4k	{
125	14.4k	int i;
126	14.4k	float ener=0, ener1=0, ener2=0;
127	14.4k	float qual=7;
128	14.4k	float log_energy;
129	14.4k	float non_st=0;
130	14.4k	float voicing;
131	14.4k	float pow_ener;
132
133	1.16M	for (i=0;i<len>>1;i++)
134	1.15M	ener1 += ((float)sig[i])*sig[i];
135
136	1.16M	for (i=len>>1;i<len;i++)
137	1.15M	ener2 += ((float)sig[i])*sig[i];
138	14.4k	ener=ener1+ener2;
139
140	14.4k	log_energy = log(ener+MIN_ENERGY);
141	86.5k	for (i=0;i<VBR_MEMORY_SIZE;i++)
142	72.1k	non_st += sqr(log_energy-vbr->last_log_energy[i]);
143	14.4k	non_st = non_st/(30*VBR_MEMORY_SIZE);
144	14.4k	if (non_st>1)
145	7.04k	non_st=1;
146
147	14.4k	voicing = 3(pitch_coef-.4)fabs(pitch_coef-.4);
148	14.4k	vbr->average_energy = 0.9vbr->average_energy + .1ener;
149	14.4k	vbr->noise_level=vbr->noise_accum/vbr->noise_accum_count;
150	14.4k	pow_ener = pow(ener,NOISE_POW);
151	14.4k	if (vbr->noise_accum_count<.06 && ener>MIN_ENERGY)
152	3.73k	vbr->noise_accum = .05*pow_ener;
153
154	14.4k	if ((voicing<.3 && non_st < .2 && pow_ener < 1.2*vbr->noise_level)
155	11.2k	\|\| (voicing<.3 && non_st < .05 && pow_ener < 1.5*vbr->noise_level)
156	11.1k	\|\| (voicing<.4 && non_st < .05 && pow_ener < 1.2*vbr->noise_level)
157	11.1k	\|\| (voicing<0 && non_st < .05))
158	3.29k	{
159	3.29k	float tmp;
160
161	3.29k	vbr->consec_noise++;
162	3.29k	if (pow_ener > 3*vbr->noise_level)
163	5	tmp = 3*vbr->noise_level;
164	3.29k	else
165	3.29k	tmp = pow_ener;
166	3.29k	if (vbr->consec_noise>=4)
167	1.09k	{
168	1.09k	vbr->noise_accum = .95vbr->noise_accum + .05tmp;
169	1.09k	vbr->noise_accum_count = .95*vbr->noise_accum_count + .05;
170	1.09k	}
171	11.1k	} else {
172	11.1k	vbr->consec_noise=0;
173	11.1k	}
174
175	14.4k	if (pow_ener < vbr->noise_level && ener>MIN_ENERGY)
176	7.32k	{
177	7.32k	vbr->noise_accum = .95vbr->noise_accum + .05pow_ener;
178	7.32k	vbr->noise_accum_count = .95*vbr->noise_accum_count + .05;
179	7.32k	}
180
181		/* Checking for very low absolute energy */
182	14.4k	if (ener < 30000)
183	4.02k	{
184	4.02k	qual -= .7;
185	4.02k	if (ener < 10000)
186	3.52k	qual-=.7;
187	4.02k	if (ener < 3000)
188	3.18k	qual-=.7;
189	10.4k	} else {
190	10.4k	float short_diff, long_diff;
191	10.4k	short_diff = log((ener+1)/(1+vbr->last_energy));
192	10.4k	long_diff = log((ener+1)/(1+vbr->average_energy));
193		/fprintf (stderr, "%f %f\n", short_diff, long_diff);/
194
195	10.4k	if (long_diff<-5)
196	1.99k	long_diff=-5;
197	10.4k	if (long_diff>2)
198	1.61k	long_diff=2;
199
200	10.4k	if (long_diff>0)
201	4.04k	qual += .6*long_diff;
202	10.4k	if (long_diff<0)
203	6.35k	qual += .5*long_diff;
204	10.4k	if (short_diff>0)
205	4.17k	{
206	4.17k	if (short_diff>5)
207	1.80k	short_diff=5;
208	4.17k	qual += 1*short_diff;
209	4.17k	}
210		/* Checking for energy increases */
211	10.4k	if (ener2 > 1.6*ener1)
212	2.84k	qual += .5;
213	10.4k	}
214	14.4k	vbr->last_energy = ener;
215	14.4k	vbr->soft_pitch = .8vbr->soft_pitch + .2pitch_coef;
216	14.4k	qual += 2.2*((pitch_coef-.4) + (vbr->soft_pitch-.4));
217
218	14.4k	if (qual < vbr->last_quality)
219	7.55k	qual = .5qual + .5vbr->last_quality;
220	14.4k	if (qual<4)
221	3.79k	qual=4;
222	14.4k	if (qual>10)
223	2.15k	qual=10;
224
225		/*
226		if (vbr->consec_noise>=2)
227		qual-=1.3;
228		if (vbr->consec_noise>=5)
229		qual-=1.3;
230		if (vbr->consec_noise>=12)
231		qual-=1.3;
232		*/
233	14.4k	if (vbr->consec_noise>=3)
234	1.57k	qual=4;
235
236	14.4k	if (vbr->consec_noise)
237	3.29k	qual -= 1.0 * (log(3.0 + vbr->consec_noise)-log(3));
238	14.4k	if (qual<0)
239	0	qual=0;
240
241	14.4k	if (ener<1600000)
242	5.92k	{
243	5.92k	if (vbr->consec_noise>2)
244	1.42k	qual-=0.5*(log(3.0 + vbr->consec_noise)-log(3));
245	5.92k	if (ener<10000&&vbr->consec_noise>2)
246	1.38k	qual-=0.5*(log(3.0 + vbr->consec_noise)-log(3));
247	5.92k	if (qual<0)
248	0	qual=0;
249	5.92k	qual += .3*log(.0001+ener/1600000.0);
250	5.92k	}
251	14.4k	if (qual<-1)
252	413	qual=-1;
253
254		/printf ("%f %f %f %f\n", qual, voicing, non_st, pow_ener/(.01+vbr->noise_level));/
255
256	14.4k	vbr->last_pitch_coef = pitch_coef;
257	14.4k	vbr->last_quality = qual;
258
259	72.1k	for (i=VBR_MEMORY_SIZE-1;i>0;i--)
260	57.7k	vbr->last_log_energy[i] = vbr->last_log_energy[i-1];
261	14.4k	vbr->last_log_energy[0] = log_energy;
262
263		/printf ("VBR: %f %f %f %f\n", (float)(log_energy-log(vbr->average_energy+MIN_ENERGY)), non_st, voicing, vbr->noise_level);/
264
265	14.4k	return qual;
266	14.4k	}
267
268		void vbr_destroy(VBRState *vbr)
269	2.51k	{
270	2.51k	}
271
272		#endif /* #ifndef DISABLE_VBR */