Big changes in both narrowband and wideband. Retrained LSP codebook,
[speexdsp.git] / libspeex / testenc.c
1 #include "speex.h"
2 #include <stdio.h>
3 #include <stdlib.h>
4
5 #define FRAME_SIZE 160
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    void *st;
17    void *dec;
18    FrameBits bits;
19
20    for (i=0;i<FRAME_SIZE;i++)
21       bak2[i]=0;
22    st = encoder_init(&speex_nb_mode);
23    dec = decoder_init(&speex_nb_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    speex_bits_init(&bits);
39    while (!feof(fin))
40    {
41       fread(in, sizeof(short), FRAME_SIZE, fin);
42       if (feof(fin))
43          break;
44       for (i=0;i<FRAME_SIZE;i++)
45          bak[i]=input[i]=in[i];
46       speex_bits_reset(&bits);
47       encode(st, input, &bits);
48       nbBits = speex_bits_write(&bits, cbits, 200);
49       printf ("Encoding frame in %d bits\n", nbBits*8);
50       if (argc==4)
51          fwrite(cbits, 1, nbBits, fbits);
52       {
53          float enoise=0, esig=0, snr;
54          for (i=0;i<FRAME_SIZE;i++)
55          {
56             enoise+=(bak2[i]-input[i])*(bak2[i]-input[i]);
57             esig += bak2[i]*bak2[i];
58          }
59          snr = 10*log10((esig+1)/(enoise+1));
60          printf ("real SNR = %f\n", snr);
61       }
62       speex_bits_rewind(&bits);
63       
64       decode(dec, &bits, input);
65
66       /* Save the bits here */
67       for (i=0;i<FRAME_SIZE;i++)
68       {
69          if (input[i]>32000)
70             input[i]=32000;
71          else if (input[i]<-32000)
72             input[i]=-32000;
73       }
74       speex_bits_reset(&bits);
75       for (i=0;i<FRAME_SIZE;i++)
76          in[i]=input[i];
77       for (i=0;i<FRAME_SIZE;i++)
78          bak2[i]=bak[i];
79       fwrite(in, sizeof(short), FRAME_SIZE, fout);
80    }
81    
82    encoder_destroy(st);
83    decoder_destroy(dec);
84    return 1;
85 }