********************************************************************
function: LSP (also called LSF) conversion routines
- last mod: $Id$
The LSP generation code is taken (with minimal modification and a
few bugfixes) from "On the Computation of the LSP Frequencies" by
float *ftmp=lsp;
int c=m>>1;
- do{
+ while(c--){
q*=ftmp[0]-w;
p*=ftmp[1]-w;
ftmp+=2;
- }while(--c);
+ }
if(m&1){
/* odd order filter; slightly assymetric */
#define EPSILON 10e-7
static int Laguerre_With_Deflation(float *a,int ord,float *r){
int i,m;
- double lastdelta=0.f;
double *defl=alloca(sizeof(*defl)*(ord+1));
for(i=0;i<=ord;i++)defl[i]=a[i];
if(delta<0.f)delta*=-1;
if(fabs(delta/new)<10e-12)break;
- lastdelta=delta;
}
r[m-1]=new;