Adding sub-band CELP (SB-CELP) -like encoding. Still incomplete.
[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    EncState st;
17    DecState dec;
18    FrameBits bits;
19
20    for (i=0;i<FRAME_SIZE;i++)
21       bak2[i]=0;
22    encoder_init(&st, &mp_nb_mode);
23    decoder_init(&dec, &mp_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    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       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       /*decode(&dec, &bits, input);*/
63
64       /* Save the bits here */
65       frame_bits_reset(&bits);
66       for (i=0;i<FRAME_SIZE;i++)
67          in[i]=input[i];
68       for (i=0;i<FRAME_SIZE;i++)
69          bak2[i]=bak[i];
70       fwrite(in, sizeof(short), FRAME_SIZE, fout);
71    }
72    
73    encoder_destroy(&st);
74    decoder_destroy(&dec);
75    return 1;
76 }