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

Routine:
  int FLfReadTF (FILE *fp, int MaxNval, int cc, float x[])

Purpose:
  Read numeric data in text form from a stream

Description:
  This procedure reads floating point numeric data from a stream.  The stream
  contains data in text form.  Data fields are free format, with data values
  separated by white-space (as defined by isspace).  Zero or more data values
  can appear in each line of input.  Commas can also be used to separate data
  values, but only within a line, i.e. a comma should not appear at the end of
  a line.  The cc argument specifies a character which marks the beginning of
  comments.

  This routine prints an error message and halts execution on detection of an
  error.

Parameters:
  <-  int FLfReadTF
      Number of values returned
   -> FILE *fp
      File pointer to the stream
   -> int MaxNval
      Maximum number of values to be returned
   -> int cc
      Comment delimiter
  <-  float x[]
      Array of output values

Author / revision:
  P. Kabal  Copyright (C) 1996
  $Revision: 1.8 $  $Date: 1996/06/28 18:43:11 $

-------------------------------------------------------------------------*/

static char rcsid[] = "$Id: FLfReadTF.c 1.8 1996/06/28 AFsp-V2R1 $";

#include <stdio.h>
#include <string.h>
#include <libtsp.h>
#include <libtsp/nucleus.h>

int
FLfReadTF (fp, MaxNval, cc, x)

     FILE *fp;
     int MaxNval;
     int cc;
     float x[];

{
  char *line;
  char *c;
  int N;
  int Nvr, Nv;

/* Read the values until end-of-file */
  N = 0;
  while (1) {

    /* Read a line of input */
    line = FLgetLine (fp);
    if (line == NULL)
      break;

    /* Process comments */
    c = strchr (line, cc);
    if (c != NULL)
      *c = '\0';		/* Force an end-of-line at the comment */

    /* Decode data records */
    Nvr = MaxNval - N;
    if (Nvr > 0) {
      if (STdecNfloat (line, 0, Nvr, &x[N], &Nv))
	UThalt ("FLfReadTF: Data format error in data file");
      N = N + Nv;
    }
    else if (*(STtrimIws (line)) != '\0')
      UThalt ("FLfReadTF: Too many values or extra data at end of data file");
  }

/* Return the number of values */
  return N;
}