view frtest/max-out.c @ 187:6ad363f7ea30

libgsmefr/d1035pf.c: perf opt
author Mychaela Falconia <falcon@freecalypso.org>
date Tue, 03 Jan 2023 03:22:04 +0000
parents a5ffec18e4cd
children
line wrap: on
line source

/*
 * This program reads an extended-libgsm binary file (which must be in FR
 * codec format, not EFR) and performs the following processing, in this
 * order:
 *
 * 1) dumps every frame in human-readable form just like gsmrec-dump;
 * 2) passes the frame stream through our DTX Rx preprocessor;
 * 3) dumps the output from the preprocessor;
 * 4) feeds this output to the standard 06.10 decoder (libgsm);
 * 5) finds the highest-magnitude PCM sample in the decoder output and
 *    reports this maximum output magnitude.
 *
 * This program is intended to serve as a debug aid, for troubleshooting
 * of our libgsmfrp+libgsm stack.
 */

#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>
#include <gsm.h>
#include "../libgsmfrp/gsm_fr_preproc.h"
#include "../libtest/binreader.h"

main(argc, argv)
	char **argv;
{
	FILE *binf;
	gsm dec_state;
	struct gsmfr_preproc_state *pp_state;
	unsigned frame_index;
	uint8_t frame[BINFILE_MAX_FRAME];
	int16_t pcm[160];
	gsm_signal fr_params[76];
	int rc, bfi, taf;
	int i, j, n;
	unsigned samp_abs, samp_max;

	if (argc != 2) {
		fprintf(stderr, "usage: %s bin-stream-file\n", argv[0]);
		exit(1);
	}
	binf = fopen(argv[1], "r");
	if (!binf) {
		perror(argv[1]);
		exit(1);
	}
	dec_state = gsm_create();
	if (!dec_state) {
		fprintf(stderr, "gsm_create() failed!\n");
		exit(1);
	}
	pp_state = gsmfr_preproc_create();
	if (!pp_state) {
		fprintf(stderr, "gsmfr_preproc_create() failed!\n");
		exit(1);
	}
	for (frame_index = 0; ; frame_index++) {
		rc = binfile_read_frame(binf, frame);
		if (rc < 0) {
			fprintf(stderr, "error: garbage in %s\n", argv[1]);
			exit(1);
		}
		if (!rc)
			break;
		printf("#%u: ", frame_index);
		switch (frame[0] & 0xF0) {
		case 0xB0:
			bfi = 1;
			taf = frame[1] & 1;
			printf("BFI TAF=%u\n", taf);
			break;
		case 0xD0:
			bfi = 0;
			fputs("FR", stdout);
			gsm_explode(dec_state, frame, fr_params);
			n = 0;
			for (i = 0; i < 8; i++)
				printf(" %d", fr_params[n++]);
			putchar('\n');
			for (i = 0; i < 4; i++) {
				putchar(' ');
				for (j = 0; j < 17; j++)
					printf(" %d", fr_params[n++]);
				putchar('\n');
			}
			break;
		default:
			fprintf(stderr, "error: %s is not in FR codec format\n",
				argv[1]);
			exit(1);
		}
		if (bfi)
			gsmfr_preproc_bfi(pp_state, taf, frame);
		else
			gsmfr_preproc_good_frame(pp_state, frame);
		fputs("Xform:", stdout);
		gsm_explode(dec_state, frame, fr_params);
		n = 0;
		for (i = 0; i < 8; i++)
			printf(" %d", fr_params[n++]);
		putchar('\n');
		for (i = 0; i < 4; i++) {
			putchar(' ');
			for (j = 0; j < 17; j++)
				printf(" %d", fr_params[n++]);
			putchar('\n');
		}
		gsm_decode(dec_state, frame, pcm);
		samp_max = 0;
		for (i = 0; i < 160; i++) {
			if (pcm[i] >= 0)
				samp_abs = pcm[i];
			else
				samp_abs = -pcm[i];
			if (samp_abs > samp_max)
				samp_max = samp_abs;
		}
		printf("  MAX=%u\n", samp_max);
	}
	exit(0);
}