- /* extract the spectral envelope and residue */
-
- {
- double floor1[n/2];
- double floor2[n/2];
- double curve[n/2];
- double work[n/2];
- double lpc1[40];
- double lpc2[40];
- double amp1,amp2;
-
- memset(floor1,0,sizeof(floor1));
- memset(floor2,0,sizeof(floor2));
- memset(work,0,sizeof(work));
-
- _vp_noise_floor(C,floor1,n/2);
- /* _vp_noise_floor(D,floor1,n/2);*/
- _vp_mask_floor(C,floor1,n/2);
- /* _vp_mask_floor(D,floor1,n/2);*/
-
- _vp_psy_sparsify(C,floor1,n/2);
- _vp_psy_sparsify(D,floor1,n/2);
-
-
- memmove(floor1+n/4,floor1,(n/4)*sizeof(double));
-
- amp2=sqrt(vorbis_curve_to_lpc(floor1,n/2,lpc2,20));
- vorbis_lpc_to_curve(work,n/2,lpc2,amp2,20);
-
-
- /* amp2=sqrt(vorbis_curve_to_lpc(floor1,n/2/8,lpc2,12));
- vorbis_lpc_to_curve(work,n/2/8,lpc2,amp2,12);
- amp2=sqrt(vorbis_curve_to_lpc(floor1+n/2/8,n/2/8,lpc2,10));
- vorbis_lpc_to_curve(work+n/2/8,n/2/8,lpc2,amp2,10);
- amp2=sqrt(vorbis_curve_to_lpc(floor1+n/2/4,n/2/4,lpc2,6));
- vorbis_lpc_to_curve(work+n/2/4,n/2/4,lpc2,amp2,6);
- amp2=sqrt(vorbis_curve_to_lpc(floor1+n/2/2,n/2/2,lpc2,6));
- vorbis_lpc_to_curve(work+n/2/2,n/2/2,lpc2,amp2,6);*/
-
- _vp_psy_quantize(C,work,n/2);
- _vp_psy_quantize(D,work,n/2);
-
-
- {
- FILE *out;
- char buffer[80];
-
-
- sprintf(buffer,"mdctC%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++)
- fprintf(out,"%g\n",fabs(C[i]));
- fclose(out);
-
- sprintf(buffer,"mdctD%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++)
- fprintf(out,"%g\n",fabs(D[i]));
- fclose(out);
-
- sprintf(buffer,"floor%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++)
- fprintf(out,"%g\n",floor1[i]);
- fclose(out);
-
-
- sprintf(buffer,"lpc%d.m",frameno);
- out=fopen(buffer,"w+");
- for(i=0;i<n/2;i++)
- fprintf(out,"%g\n",work[i]);
- fclose(out);
-
- frameno++;
- }
+#ifdef ANALYSIS
+int analysis_noisy=1;
+
+/* there was no great place to put this.... */
+void _analysis_output_always(char *base,int i,float *v,int n,int bark,int dB,ogg_int64_t off){
+ int j;
+ FILE *of;
+ char buffer[80];
+
+ sprintf(buffer,"%s_%d.m",base,i);
+ of=fopen(buffer,"w");
+
+ if(!of)perror("failed to open data dump file");
+
+ for(j=0;j<n;j++){
+ if(bark){
+ float b=toBARK((4000.f*j/n)+.25);
+ fprintf(of,"%f ",b);
+ }else
+ if(off!=0)
+ fprintf(of,"%f ",(double)(j+off)/8000.);
+ else
+ fprintf(of,"%f ",(double)j);
+
+ if(dB){
+ float val;
+ if(v[j]==0.)
+ val=-140.;
+ else
+ val=todB(v+j);
+ fprintf(of,"%f\n",val);
+ }else{
+ fprintf(of,"%f\n",v[j]);