Updated Blackfin version of compute_pitch_error()
[speexdsp.git] / ti / testenc-TI-C64x.c
1 /* Copyright (C) 2005 Psi Systems, Inc.
2    Author:  Jean-Marc Valin 
3    File: testenc-TI-C64x.c
4    Encoder/Decoder Loop Main file for TI TMS320C64xx processor\r
5    for use with TI Code Composer (TM) DSP development tools.
6    Modified from speexlib/testenc.c\r
7 \r
8
9    Redistribution and use in source and binary forms, with or without
10    modification, are permitted provided that the following conditions
11    are met:
12    
13    - Redistributions of source code must retain the above copyright
14    notice, this list of conditions and the following disclaimer.
15    
16    - Redistributions in binary form must reproduce the above copyright
17    notice, this list of conditions and the following disclaimer in the
18    documentation and/or other materials provided with the distribution.
19    
20    - Neither the name of the Xiph.org Foundation nor the names of its
21    contributors may be used to endorse or promote products derived from
22    this software without specific prior written permission.
23    
24    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
25    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
26    LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
27    A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR
28    CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
29    EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
30    PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
31    PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
32    LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
33    NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
34    SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
35 */
36 \r
37 #ifdef HAVE_CONFIG_H\r
38 #include "config.h"\r
39 #endif\r
40 \r
41 #include <speex/speex.h>\r
42 #include <stdio.h>\r
43 #include <stdlib.h>\r
44 #include <speex/speex_callbacks.h>\r
45 \r
46 #undef DECODE_ONLY\r
47 #define CHECK_RESULT  /* Compares original file with encoded/decoder file */\r
48 #define TESTENC_BYTES_PER_FRAME 20  /* 8kbps */\r
49 #define TESTENC_QUALITY 4 /* 8kbps */\r
50 //#define TESTENC_BYTES_PER_FRAME 28  /* 11kbps */\r
51 //#define TESTENC_QUALITY 5           /* 11 kbps */\r
52 \r
53 /* For narrowband, QUALITY maps to these bit rates (see modes.c, manual.pdf) \r
54  *   {1, 8, 2, 3, 3, 4, 4, 5, 5, 6, 7}\r
55  * 0 -> 2150\r
56  * 1 -> 3950\r
57  * 2 -> 5950\r
58  * 3 -> 8000\r
59  * 4 -> 8000\r
60  * 5 -> 11000\r
61  * 6 -> 11000\r
62  * 7 -> 15000\r
63  * 8 -> 15000\r
64  * 9 -> 18200\r
65  *10 -> 26400  */\r
66 \r
67 #ifdef FIXED_DEBUG\r
68 extern long long spx_mips;\r
69 #endif\r
70 \r
71 #ifdef MANUAL_ALLOC\r
72 /* Take all Speex space from this private heap */ \r
73 /* This is useful for multichannel applications */\r
74 #pragma DATA_SECTION(spxHeap, ".myheap"); \r
75 static char spxHeap[SPEEX_PERSIST_STACK_SIZE];\r
76 \r
77 #pragma DATA_SECTION(spxScratch, ".myheap"); \r
78 static char spxScratch[SPEEX_SCRATCH_STACK_SIZE];\r
79 \r
80 char *spxGlobalHeapPtr, *spxGlobalHeapEnd;\r
81 char *spxGlobalScratchPtr, *spxGlobalScratchEnd;\r
82 #endif    /* MANUAL_ALLOC */\r
83 \r
84 #include <math.h>\r
85 void main()\r
86 {\r
87    char *outFile, *bitsFile;\r
88    FILE *fout, *fbits=NULL;\r
89 #if !defined(DECODE_ONLY) || defined(CHECK_RESULT)\r
90    char *inFile;\r
91    FILE *fin;\r
92    short in_short[FRAME_SIZE];\r
93 #endif\r
94    short out_short[FRAME_SIZE];\r
95 #ifndef DECODE_ONLY\r
96    int nbBits;\r
97 #endif\r
98 #ifdef CHECK_RESULT\r
99    float sigpow,errpow,snr, seg_snr=0;\r
100    int snr_frames = 0;\r
101    int i;\r
102 #endif\r
103    char cbits[200];\r
104    void *st;\r
105    void *dec;\r
106    SpeexBits bits;\r
107    int tmp;\r
108    unsigned long bitCount=0;\r
109    int skip_group_delay;\r
110    SpeexCallback callback;\r
111 \r
112 #ifdef CHECK_RESULT\r
113    sigpow = 0;\r
114    errpow = 0;\r
115 #endif\r
116 \r
117 #ifdef MANUAL_ALLOC\r
118         spxGlobalHeapPtr = spxHeap;\r
119         spxGlobalHeapEnd = spxHeap + sizeof(spxHeap);\r
120 \r
121         spxGlobalScratchPtr = spxScratch;\r
122         spxGlobalScratchEnd = spxScratch + sizeof(spxScratch);\r
123 #endif\r
124    st = speex_encoder_init(&speex_nb_mode);\r
125 #ifdef MANUAL_ALLOC\r
126         spxGlobalScratchPtr = spxScratch;               /* Reuse scratch for decoder */\r
127 #endif\r
128    dec = speex_decoder_init(&speex_nb_mode);\r
129 \r
130    callback.callback_id = SPEEX_INBAND_CHAR;\r
131    callback.func = speex_std_char_handler;\r
132    callback.data = stderr;\r
133    speex_decoder_ctl(dec, SPEEX_SET_HANDLER, &callback);\r
134 \r
135    callback.callback_id = SPEEX_INBAND_MODE_REQUEST;\r
136    callback.func = speex_std_mode_request_handler;\r
137    callback.data = st;\r
138    speex_decoder_ctl(dec, SPEEX_SET_HANDLER, &callback);\r
139 \r
140    tmp=0;\r
141    speex_decoder_ctl(dec, SPEEX_SET_ENH, &tmp);\r
142    tmp=0;\r
143    speex_encoder_ctl(st, SPEEX_SET_VBR, &tmp);\r
144    tmp=TESTENC_QUALITY;\r
145    speex_encoder_ctl(st, SPEEX_SET_QUALITY, &tmp);\r
146    tmp=1;  /* Lowest */\r
147    speex_encoder_ctl(st, SPEEX_SET_COMPLEXITY, &tmp);\r
148 \r
149    speex_mode_query(&speex_nb_mode, SPEEX_MODE_FRAME_SIZE, &tmp);\r
150    fprintf (stderr, "frame size: %d\n", tmp);\r
151    skip_group_delay = tmp / 4;          /* 5ms algorithmic delay */\r
152 \r
153 #ifdef DECODE_ONLY\r
154    bitsFile = "c:\\speextrunktest\\samples\\malebitsin.dat";\r
155    fbits = fopen(bitsFile, "rb");\r
156 #else\r
157    bitsFile = "c:\\speextrunktest\\samples\\malebits.dat";\r
158    fbits = fopen(bitsFile, "wb");\r
159 #endif\r
160 #if !defined(DECODE_ONLY) || defined(CHECK_RESULT)\r
161    inFile = "c:\\speextrunktest\\samples\\male.snd";\r
162    fin = fopen(inFile, "rb");\r
163 #endif\r
164    outFile = "c:\\speextrunktest\\samples\\maleout.snd";\r
165    fout = fopen(outFile, "wb+");\r
166  \r
167    speex_bits_init(&bits);\r
168 #ifndef DECODE_ONLY\r
169    while (!feof(fin))\r
170    {\r
171       fread(in_short, sizeof(short), FRAME_SIZE, fin);\r
172       if (feof(fin))\r
173          break;\r
174       speex_bits_reset(&bits);\r
175 \r
176       speex_encode_int(st, in_short, &bits);\r
177       nbBits = speex_bits_write(&bits, cbits, 200);\r
178       bitCount+=bits.nbBits;\r
179 \r
180       fwrite(cbits, 1, nbBits, fbits);\r
181       speex_bits_rewind(&bits);\r
182 \r
183 #else /* DECODE_ONLY */\r
184    while (!feof(fbits))\r
185    {\r
186       fread(cbits, 1, TESTENC_BYTES_PER_FRAME, fbits);\r
187 \r
188       if (feof(fbits))\r
189          break;\r
190 \r
191       speex_bits_read_from(&bits, cbits, TESTENC_BYTES_PER_FRAME);\r
192 //      bitCount+=160;  /* only correct for 8kbps, but just for the printf */\r
193       bitCount+=bits.nbBits;\r
194 #endif\r
195       \r
196       speex_decode_int(dec, &bits, out_short);\r
197       speex_bits_reset(&bits);\r
198 \r
199       fwrite(&out_short[skip_group_delay], sizeof(short), FRAME_SIZE-skip_group_delay, fout);\r
200       skip_group_delay = 0;\r
201 #if 1\r
202    fprintf (stderr, "Bits so far: %lu \n", bitCount);\r
203 #endif\r
204    }\r
205    fprintf (stderr, "Total encoded size: %lu bits\n", bitCount);\r
206    speex_encoder_destroy(st);\r
207    speex_decoder_destroy(dec);\r
208 \r
209 #ifdef CHECK_RESULT \r
210    rewind(fin);\r
211    rewind(fout);\r
212 \r
213    while ( FRAME_SIZE == fread(in_short, sizeof(short), FRAME_SIZE, fin) \r
214            &&\r
215            FRAME_SIZE ==  fread(out_short, sizeof(short), FRAME_SIZE,fout) )\r
216    {\r
217         float s=0, e=0;\r
218         for (i=0;i<FRAME_SIZE;++i) {\r
219             s += (float)in_short[i] * in_short[i];\r
220             e += ((float)in_short[i]-out_short[i]) * ((float)in_short[i]-out_short[i]);\r
221         }\r
222         seg_snr += 10*log10((s+160)/(e+160));\r
223         sigpow += s;\r
224         errpow += e;\r
225         snr_frames++;\r
226    }\r
227 \r
228    snr = 10 * log10( sigpow / errpow );\r
229    seg_snr /= snr_frames;\r
230    fprintf(stderr,"SNR = %f\nsegmental SNR = %f\n",snr, seg_snr);\r
231 \r
232 #ifdef FIXED_DEBUG\r
233    printf ("Total: %f MIPS\n", (float)(1e-6*50*spx_mips/snr_frames));\r
234 #endif\r
235 #endif /* CHECK_RESULT */\r
236 #if !defined(DECODE_ONLY) || defined(CHECK_RESULT)\r
237    fclose(fin);\r
238 #endif\r
239    fclose(fout);\r
240    fclose(fbits);\r
241 }\r