+++ /dev/null
-/*
-** Copyright (C) 1999-2011 Erik de Castro Lopo <erikd@mega-nerd.com>
-**
-** This program is free software; you can redistribute it and/or modify
-** it under the terms of the GNU General Public License as published by
-** the Free Software Foundation; either version 2 of the License, or
-** (at your option) any later version.
-**
-** This program is distributed in the hope that it will be useful,
-** but WITHOUT ANY WARRANTY; without even the implied warranty of
-** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-** GNU General Public License for more details.
-**
-** You should have received a copy of the GNU General Public License
-** along with this program; if not, write to the Free Software
-** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
-*/
-
-#include <stdio.h>
-#include <unistd.h>
-#include <stdlib.h>
-#include <string.h>
-#include <math.h>
-
-#include "g72x.h"
-#include "g72x_priv.h"
-
-#ifndef M_PI
-#define M_PI 3.14159265358979323846264338
-#endif
-
-#define BUFFER_SIZE (1<<14) /* Should be (1<<14) */
-#define SAMPLE_RATE 11025
-
-
-static void g721_test (void) ;
-static void g723_test (double margin) ;
-
-static void gen_signal_double (double *data, double scale, int datalen) ;
-static int error_function (double data, double orig, double margin) ;
-
-static int oct_save_short (short *a, short *b, int len) ;
-
-int
-main (int argc, char *argv [])
-{ int bDoAll = 0 ;
- int nTests = 0 ;
-
- if (argc != 2)
- { printf ("Usage : %s <test>\n", argv [0]) ;
- printf (" Where <test> is one of the following:\n") ;
- printf (" g721 - test G721 encoder and decoder\n") ;
- printf (" g723 - test G721 encoder and decoder\n") ;
- printf (" all - perform all tests\n") ;
- exit (1) ;
- } ;
-
- bDoAll=!strcmp (argv [1], "all");
-
- if (bDoAll || ! strcmp (argv [1], "g721"))
- { g721_test () ;
- nTests++ ;
- } ;
-
- if (bDoAll || ! strcmp (argv [1], "g723"))
- { g723_test (0.53) ;
- nTests++ ;
- } ;
-
- if (nTests == 0)
- { printf ("Mono : ************************************\n") ;
- printf ("Mono : * No '%s' test defined.\n", argv [1]) ;
- printf ("Mono : ************************************\n") ;
- return 1 ;
- } ;
-
- return 0 ;
-} /* main */
-
-static void
-g721_test (void)
-{
- return ;
-} /* g721_test */
-
-static void
-g723_test (double margin)
-{ static double orig_buffer [BUFFER_SIZE] ;
- static short orig [BUFFER_SIZE] ;
- static short data [BUFFER_SIZE] ;
-
- G72x_STATE encoder_state, decoder_state ;
-
- long k ;
- int code, position, max_err ;
-
- private_init_state (&encoder_state) ;
- encoder_state.encoder = g723_24_encoder ;
- encoder_state.codec_bits = 3 ;
-
- private_init_state (&decoder_state) ;
- decoder_state.decoder = g723_24_decoder ;
- decoder_state.codec_bits = 3 ;
-
- memset (data, 0, BUFFER_SIZE * sizeof (short)) ;
- memset (orig, 0, BUFFER_SIZE * sizeof (short)) ;
-
- printf (" g723_test : ") ;
- fflush (stdout) ;
-
- gen_signal_double (orig_buffer, 32000.0, BUFFER_SIZE) ;
- for (k = 0 ; k < BUFFER_SIZE ; k++)
- orig [k] = (short) orig_buffer [k] ;
-
- /* Write and read data here. */
- position = 0 ;
- max_err = 0 ;
- for (k = 0 ; k < BUFFER_SIZE ; k++)
- { code = encoder_state.encoder (orig [k], &encoder_state) ;
- data [k] = decoder_state.decoder (code, &decoder_state) ;
- if (abs (orig [k] - data [k]) > max_err)
- { position = k ;
- max_err = abs (orig [k] - data [k]) ;
- } ;
- } ;
-
- printf ("\n\nMax error of %d at postion %d.\n", max_err, position) ;
-
- for (k = 0 ; k < BUFFER_SIZE ; k++)
- { if (error_function (data [k], orig [k], margin))
- { printf ("Line %d: Incorrect sample A (#%ld : %d should be %d).\n", __LINE__, k, data [k], orig [k]) ;
- oct_save_short (orig, data, BUFFER_SIZE) ;
- exit (1) ;
- } ;
- } ;
-
-
- printf ("ok\n") ;
-
- return ;
-} /* g723_test */
-
-
-#define SIGNAL_MAXVAL 30000.0
-#define DECAY_COUNT 1000
-
-static void
-gen_signal_double (double *gendata, double scale, int gendatalen)
-{ int k, ramplen ;
- double amp = 0.0 ;
-
- ramplen = DECAY_COUNT ;
-
- for (k = 0 ; k < gendatalen ; k++)
- { if (k <= ramplen)
- amp = scale * k / ((double) ramplen) ;
- else if (k > gendatalen - ramplen)
- amp = scale * (gendatalen - k) / ((double) ramplen) ;
-
- gendata [k] = amp * (0.4 * sin (33.3 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE))
- + 0.3 * cos (201.1 * 2.0 * M_PI * ((double) (k+1)) / ((double) SAMPLE_RATE))) ;
- } ;
-
- return ;
-} /* gen_signal_double */
-
-static int
-error_function (double data, double orig, double margin)
-{ double error ;
-
- if (fabs (orig) <= 500.0)
- error = fabs (fabs (data) - fabs(orig)) / 2000.0 ;
- else if (fabs (orig) <= 1000.0)
- error = fabs (data - orig) / 3000.0 ;
- else
- error = fabs (data - orig) / fabs (orig) ;
-
- if (error > margin)
- { printf ("\n\n*******************\nError : %f\n", error) ;
- return 1 ;
- } ;
- return 0 ;
-} /* error_function */
-
-static int
-oct_save_short (short *a, short *b, int len)
-{ FILE *file ;
- int k ;
-
- if (! (file = fopen ("error.dat", "w")))
- return 1 ;
-
- fprintf (file, "# Not created by Octave\n") ;
-
- fprintf (file, "# name: a\n") ;
- fprintf (file, "# type: matrix\n") ;
- fprintf (file, "# rows: %d\n", len) ;
- fprintf (file, "# columns: 1\n") ;
-
- for (k = 0 ; k < len ; k++)
- fprintf (file, "% d\n", a [k]) ;
-
- fprintf (file, "# name: b\n") ;
- fprintf (file, "# type: matrix\n") ;
- fprintf (file, "# rows: %d\n", len) ;
- fprintf (file, "# columns: 1\n") ;
-
- for (k = 0 ; k < len ; k++)
- fprintf (file, "% d\n", b [k]) ;
-
- fclose (file) ;
- return 0 ;
-} /* oct_save_short */
-
projects / Faustine.git / blobdiff