Not a permanent file, but I don;t want to lose it yet.

Monty 19990807

svn path=/trunk/vorbis/; revision=22

diff --git a/lib/dsptest.c b/lib/dsptest.c
new file mode 100644 (file)
index 0000000..33f6fae
--- /dev/null
+++ b/lib/dsptest.c
@@ -0,0 +1,236 @@
+#include <math.h>
+#include <stdio.h>
+#include "codec.h"
+#include "envelope.h"
+#include "mdct.h"
+
+#define READ 4096
+
+int main(){
+   vorbis_dsp_state encode,decode;
+   vorbis_info vi;
+   vorbis_block vb;
+   long counterin=0;
+   long counterout=0;
+   int done=0;
+   char *temp[]={ "Test" ,"the Test band", "test records",NULL };
+   int mtemp[]={0,1};
+   int frame=0;
+
+   signed char buffer[READ*4+44];
+   
+  
+   vi.channels=2;
+   vi.rate=44100;
+   vi.version=0;
+   vi.mode=0;
+   vi.user_comments=temp;
+   vi.vendor="Xiphophorus";
+   vi.smallblock=512;
+   vi.largeblock=2048;
+   vi.envelopesa=64;
+   vi.envelopech=2;
+   vi.envelopemap=mtemp;
+   vi.channelmap=NULL;
+   vi.preecho_thresh=10.;
+   vi.preecho_clamp=4.;
+
+   vorbis_analysis_init(&encode,&vi);
+   vorbis_synthesis_init(&decode,&vi);
+   vorbis_block_init(&encode,&vb);
+ 
+   fread(buffer,1,44,stdin);
+   fwrite(buffer,1,44,stdout);
+   
+   
+   while(!done){
+     long bread=fread(buffer,1,READ*4,stdin);
+     double **buf=vorbis_analysis_buffer(&encode,READ);
+     long i;
+     
+     /* uninterleave samples */
+     
+     for(i=0;i<bread/4;i++){
+       buf[0][i]=(int)((buffer[i*4+1]<<8)|(0x00ff&(int)buffer[i*4]));
+       buf[1][i]=(int)((buffer[i*4+3]<<8)|(0x00ff&(int)buffer[i*4+2]));
+     }
+  
+     vorbis_analysis_wrote(&encode,i);
+ 
+     while(vorbis_analysis_blockout(&encode,&vb)){
+       double **pcm;
+       int avail;
+       double *window=encode.window[vb.W][vb.lW][vb.nW];
+ 
+       /* analysis */
+
+       /* apply envelope */
+       _ve_envelope_sparsify(&vb);
+       _ve_envelope_apply(&vb,0);
+ 
+       for(i=0;i<vb.pcm_channels;i++)
+        mdct_forward(&vb.vd->vi[vb.W],vb.pcm[i],vb.pcm[i],window);
+ 
+
+       /* synthesis */
+
+       for(i=0;i<vb.pcm_channels;i++)
+        mdct_backward(&vb.vd->vi[vb.W],vb.pcm[i],vb.pcm[i],window);
+ 
+ 
+       /*{
+        FILE *out;
+        char path[80];
+        int i;
+        
+        int avail=encode.block_size[vb.W];
+        int beginW=countermid-avail/2;
+        
+        sprintf(path,"ana%d",vb.frameno);
+        out=fopen(path,"w");
+        
+        for(i=0;i<avail;i++)
+          fprintf(out,"%d %g\n",i+beginW,vb.pcm[0][i]);
+        fprintf(out,"\n");
+        for(i=0;i<avail;i++)
+          fprintf(out,"%d %g\n",i+beginW,window[i]);
+        
+        fclose(out);
+       }*/
+ 
+       _ve_envelope_apply(&vb,1);
+
+       counterin+=bread/4;
+       vorbis_synthesis_blockin(&decode,&vb);
+        
+       while((avail=vorbis_synthesis_pcmout(&decode,&pcm))){
+        if(avail>READ)avail=READ;
+
+        for(i=0;i<avail;i++){
+          int l=rint(pcm[0][i]*32768.);
+          int r=rint(pcm[1][i]*32768.);
+          if(l>32767)l=32767;
+          if(r>32767)r=32767;
+          if(l<-32768)l=-32768;
+          if(r<-32768)r=-32768;
+          buffer[i*4]=l&0xff;
+          buffer[i*4+1]=(l>>8)&0xff;
+          buffer[i*4+2]=r&0xff;
+          buffer[i*4+3]=(r>>8)&0xff;
+        }    
+        
+        fwrite(buffer,1,avail*4,stdout);
+        
+        /*{
+          FILE *out;
+          char path[80];
+          int i;
+          
+          sprintf(path,"syn%d",frame);
+          out=fopen(path,"w");
+          
+          for(i=0;i<avail;i++)
+            fprintf(out,"%ld %g\n",i+counterout,pcm[0][i]);
+          fprintf(out,"\n");
+          for(i=0;i<avail;i++)
+            fprintf(out,"%ld %g\n",i+counterout,pcm[1][i]);
+          
+          fclose(out);
+          }*/
+
+        vorbis_synthesis_read(&decode,avail);
+        
+        counterout+=avail;
+        frame++;
+       }
+       
+ 
+       if(vb.eofflag){
+        done=1;
+        break;
+       }
+     }
+   }
+   return 0;
+}
+
+
+
+
+ #endif
+
+  signed char buffer[BANDS*4+44];
+  long count=0;
+  double max=0;
+  int imax=0;
+  
+  double left[BANDS+POINTS],right[BANDS+POINTS];
+
+  fread(buffer,1,44,stdin);
+  fwrite(buffer,1,44,stdout);
+
+
+  while(!feof(stdin)){
+    long bread=fread(buffer,1,BANDS*4,stdin);
+    long i;
+
+    /* uninterleave samples */
+    
+    for(i=0;i<BANDS;i++){
+      left[i]=(int)((buffer[i*4+1]<<8)|(0x00ff&(int)buffer[i*4]));
+      right[i]=(int)((buffer[i*4+3]<<8)|(0x00ff&(int)buffer[i*4+2]));
+    }
+  
+    /*{
+      FILE *out=fopen("in.m","w");
+      for(i=0;i<BANDS;i++)
+       fprintf(out,"%g\n",left[i]);
+      fclose(out);
+    }*/
+
+    /* analysis bank */
+      
+    input_subbander(filterIL,left);
+    input_subbander(filterIR,right);
+  
+    ppQMF_subband_b32_w512(filterIL,left);
+    ppQMF_subband_b32_w512(filterIR,right);
+
+    for(i=0;i<BANDS;i++){
+      double c=left[i]+right[i];
+      double d=(left[i]-right[i])*4.; /* Stereo expansion */
+      
+      left[i]=(c+d)/2.;
+      right[i]=(c-d)/2.;
+
+      left[i]=alog(left[i]); 
+      right[i]=alog(right[i]);
+
+    }
+
+    ppQMF_synth_b32_w512(filterOL,left);
+    ppQMF_synth_b32_w512(filterOR,right);
+    
+    output_subbander(filterOL,left);
+    output_subbander(filterOR,right);
+
+    for(i=0;i<BANDS;i++){
+      int l=rint(left[i]/80.);
+      int r=rint(right[i]/80.);
+      if(abs(l)>32767 || abs(r)>32768){
+       fprintf(stderr,"\nClipping!\n");
+       exit(0);
+      }
+      buffer[i*4]=l&0xff;
+      buffer[i*4+1]=(l>>8)&0xff;
+      buffer[i*4+2]=r&0xff;
+      buffer[i*4+3]=(r>>8)&0xff;
+    }    
+
+    fwrite(buffer,1,bread,stdout);
+  }
+  return(0);
+}
+
+
+