-/* One can do this the long way by generating the transfer function in
- the time domain and taking the forward FFT of the result. The
- results from direct calculation are cleaner and faster, however */
-
-static double vorbis_lpc_magnitude(double w,double *lpc, int m){
- int k;
- double real=1,imag=0;
- for(k=0;k<m;k++){
- real+=lpc[k]*cos((k+1)*w);
- imag+=lpc[k]*sin((k+1)*w);
- }
- return(1./sqrt(real*real+imag*imag));
-}
-
-/* generate the whole freq response curve on an LPC IIR filter */
+void vorbis_lpc_predict(float *coeff,float *prime,int m,
+ float *data,long n){