SB-CELP work... more to go
[speexdsp.git] / libspeex / testenc_sb.c
1 #include "sb_celp.h"
2 #include <stdio.h>
3 #include <stdlib.h>
4
5 #define FRAME_SIZE 320
6 #include <math.h>
7 int main(int argc, char **argv)
8 {
9    char *inFile, *outFile, *bitsFile;
10    FILE *fin, *fout, *fbits=NULL;
11    short in[FRAME_SIZE];
12    float input[FRAME_SIZE], bak[FRAME_SIZE], bak2[FRAME_SIZE];
13    char cbits[200];
14    int nbBits;
15    int i;
16    SBEncState st;
17    SBDecState dec;
18    FrameBits bits;
19
20    for (i=0;i<FRAME_SIZE;i++)
21       bak2[i]=0;
22    sb_encoder_init(&st, &mp_sb_mode);
23    sb_decoder_init(&dec, &mp_sb_mode);
24    if (argc != 4 && argc != 3)
25    {
26       fprintf (stderr, "Usage: encode [in file] [out file] [bits file]\nargc = %d", argc);
27       exit(1);
28    }
29    inFile = argv[1];
30    fin = fopen(inFile, "r");
31    outFile = argv[2];
32    fout = fopen(outFile, "w");
33    if (argc==4)
34    {
35       bitsFile = argv[3];
36       fbits = fopen(bitsFile, "w");
37    }
38    frame_bits_init(&bits);
39    while (!feof(fin))
40    {
41       fread(in, sizeof(short), FRAME_SIZE, fin);
42       for (i=0;i<FRAME_SIZE;i++)
43          bak[i]=input[i]=in[i];
44       frame_bits_reset(&bits);
45       sb_encode(&st, input, &bits);
46       nbBits = frame_bits_write(&bits, cbits, 200);
47       printf ("Encoding frame in %d bits\n", nbBits*8);
48       if (argc==4)
49          fwrite(cbits, 1, nbBits, fbits);
50       {
51          float enoise=0, esig=0, snr;
52          for (i=0;i<FRAME_SIZE;i++)
53          {
54             enoise+=(bak2[i]-input[i])*(bak2[i]-input[i]);
55             esig += bak2[i]*bak2[i];
56          }
57          snr = 10*log10((esig+1)/(enoise+1));
58          printf ("real SNR = %f\n", snr);
59       }
60       frame_bits_rewind(&bits);
61       
62       sb_decode(&dec, &bits, input);
63
64       frame_bits_reset(&bits);
65       for (i=0;i<FRAME_SIZE;i++)
66       {
67          if (input[i]>32000)
68             input[i]=32000;
69          else if (input[i]<-32000)
70             input[i]=-32000;
71       }
72       for (i=0;i<FRAME_SIZE;i++)
73          in[i]=input[i];
74       for (i=0;i<FRAME_SIZE;i++)
75          bak2[i]=bak[i];
76       fwrite(in, sizeof(short), FRAME_SIZE, fout);
77    }
78    
79    sb_encoder_destroy(&st);
80    sb_decoder_destroy(&dec);
81    return 1;
82 }