/*-------------- Telecommunications & Signal Processing Lab ---------------
                             McGill University

Routine:
  void LPlpcPar (const char Fname[], float *pre, const float **Win, int *Lwin,
                 int *Woffs, int *Lframe, int *Np, float *bwexp)

Purpose:
  Read LPC analysis/sythesis parameters

Description:
  This routine reads LPC analysis parameters from a parameter file.

Parameters:
   -> const char Fname[]
      Parameter file
  <-  float *pre
      Preemphasis factor
  <-  const float **Win
      Array of window coefficients
  <-  int *Lwin
      Number of window coefficients
  <-  int *Woffs
      Window offset relative to frame
  <-  int *Lframe
      Frame length
  <-  int *Np
      Number of LPC coefficients
  <-  float *bwexp
      Bandwidth expansion factor

Author / revision:
  P. Kabal  Copyright (C) 1996
  $Revision: 1.8 $  $Date: 1996/07/16 14:25:02 $

-------------------------------------------------------------------------*/

static char rcsid[] = "$Id: LPlpcPar.c 1.8 1996/07/16 AFsp-V2R1 $";

#include <stdio.h>
#include <libtsp.h>

#ifndef LPSYN
#include "LPanal.h"
#else
#include "LPsyn.h"
#endif

enum {
  HAMMING = 0,
  RECT = 1
};

#define ERRSTOP(text,par)	UThalt ("%s: %s: \"%s\"", PROGRAM, text, par)

#define NENTRY	7
static const char *keytable[NENTRY+1]={
  "pre*emphasis_factor",
  "window_len*gth",
  "window_off*set",
  "window_t*ype",
  "frame_len*gth",
  "LPC_n*umber",
  "LPC_BW*_expansion",
  NULL
};
static const char *wkey[]={
  "H*amming",
  "r*ectangular",
  NULL
};

void
LPlpcPar (Fname, pre, Win, Lwin, Woffs, Lframe, Np, bwexp)

     const char Fname[];
     float *pre;
     const float **Win;
     int *Lwin;
     int *Woffs;
     int *Lframe;
     int *Np;
     float *bwexp;

{
  FILE *fp;
  char *Line;
  int ind;
  int Wtype;
  float prex, bwexpx;
  int Lwinx, Woffsx, Lframex, Npx;
  char line[80];
  static float Wdata[MAXWINDOW];

/* Default values */
  prex = 0.0;
  Lwinx = -1;
  Woffsx = 0;
  Wtype = HAMMING;
  Lframex = -1;
  Npx = -1;
  bwexpx = 1.0;

/* Open the parameter file */
  fp = fopen (Fname, "r");
  if (fp == NULL)
    UTerror ("%s: Cannot open file \"%s\"", PROGRAM, Fname);

/* Decode parameters */
  while (1) {
    Line = FLgetRec (fp, "!\0\\", 1);
    if (Line == NULL)
      break;
    ind = STkeyXpar (Line, "=", "\"\"", keytable, Line);
    if (ind < 0)
      break;
    switch (ind) {
    case 0:
      if (STdec1float (Line, &prex))
	ERRSTOP ("Invalid preemphasis value", Line);
      break;
    case 1:
      if (STdec1int (Line, &Lwinx) || Lwinx < 0 || Lwinx > MAXWINDOW)
	ERRSTOP ("Invalid window length", Line);
      break;
    case 2:
      if (STdec1int (Line, &Woffsx))
	ERRSTOP ("Invalid window offset value", Line);
      break;
    case 3:
      Wtype = STkeyMatch (Line, wkey);
      if (Wtype < 0)
	ERRSTOP ("Invalid window type", Line);
      break;
    case 4:
      if (STdec1int (Line, &Lframex) || Lframex < 0 || Lframex > MAXFRAME)
	ERRSTOP ("Invalid frame length", Line);
      break;
    case 5:
      if (STdec1int (Line, &Npx) || Npx < 0 || Npx > MAXNP)
	ERRSTOP ("Invalid number of LPC coefficients", Line);
      break;
    case 6:
      if (STdec1float (Line, &bwexpx))
	ERRSTOP ("Invalid LPC bandwidth expansion factor", Line);
      break;
    }
  }
  fclose (fp);

/* Set up the analysis window */
  if (Wtype == HAMMING)
    FIwinHamm (Wdata, Lwinx, 0.46);
  else
    FIwinHamm (Wdata, Lwinx, 0.0);

/* Error checks */
  if (Lframex < 0)
    UThalt ("%s: Frame length not set", PROGRAM);
  if (Npx < 0)
    UThalt ("%s: No. LPC coefficients not set", PROGRAM);
  if (Lwinx < 0)
    UThalt ("%s: Window length not set", PROGRAM);

/* Write the parameters to the header information string */
  AFsetHinfo ("\nparameters:\\");
  sprintf (line, "\n  preemphasis_factor = %g\\", prex);
  AFsetHinfo (line);
  sprintf (line, "\n  window_length = %d\\", Lwinx);
  AFsetHinfo (line);
  sprintf (line, "\n  window_offset = %d\\", Woffsx);
  AFsetHinfo (line);
  if (Wtype == HAMMING)
    AFsetHinfo ("\n  window_type = Hamming\\");
  else
    AFsetHinfo ("\n  window_type = rectangular\\");
  sprintf (line, "\n  frame_length = %d\\", Lframex);
  AFsetHinfo (line);
  sprintf (line, "\n  LPC_number = %d\\", Npx);
  AFsetHinfo (line);
  sprintf (line, "\n  LPC_BW_expansion = %g", bwexpx);
  AFsetHinfo (line);

/* Return the values */
  *pre = prex;
  *Win = Wdata;
  *Lwin = Lwinx;
  *Woffs = Woffsx;
  *Lframe = Lframex;
  *Np = Npx;
  *bwexp = bwexpx;

  return;
}