/*
 *  ===================================================================
 *  prog_dataset.c : This program reads and stores an unknown
 *                   one-dimensional stream of data from an input file.
 * 
 *  Copyright (C) 1997 by Mark Austin.
 *                                                                     
 *  This software is provided "as is" without express or implied warranty.
 *  Permission is granted to use this software on any computer system,
 *  and to redistribute it freely, subject to the following restrictions:
 * 
 *  1. The author is not responsible for the consequences of use of
 *     this software, even if they arise from defects in the software.
 *  2. The origin of this software must not be misrepresented, either
 *     by explicit claim or by omission.
 *  3. Altered versions must be plainly marked as such, and must not
 *     be misrepresented as being the original software.
 *  4. This notice is to remain intact.
 * 
 *  Written By : Mark Austin                            September 1997 
 *  ===================================================================
 */

#include <stdio.h>
#include <stdlib.h> 
#include "miscellaneous.h"

void vecPrint( char *, float *, int );   /* declaration for vector function */

int main( void ) {
enum { InitialLength = 3, Increment = 2 };
float     *fpData;
float  fDataPoint;
int iNoPoints = 0;
int iVectorLength;

    /* [a] : Allocate and print initial data vector */

    fpData = (float *) safeCalloc( InitialLength, sizeof (float),
                                   __FILE__, __LINE__ );

    iVectorLength = InitialLength;

    vecPrint( "Initial data vector", fpData , iVectorLength );

    /* [b] : Read data from standard input until an EOF is obtained */

    while ( scanf( "%f%*c", &fDataPoint ) != EOF ) {

        iNoPoints = iNoPoints + 1;
        if( iNoPoints > iVectorLength ) {

            iVectorLength += Increment;

            fpData = (float *) realloc( fpData ,
                                        iVectorLength*sizeof (float) );
            if(fpData == NULL) {
               printf("ERROR >> Memory Reallocation Failure\n");
               exit(1);
            }

            vecPrint( "Reallocated Data Vector", fpData , iVectorLength );
        }

        fpData[ iNoPoints - 1 ] = fDataPoint;
    }

    vecPrint( "Final Data Vector", fpData , iVectorLength );

    /* [c] : Free vector memory */

    free( fpData );
}

/*
 *  ================================================================
 *  vecPrint() : Naive Implementation of Vector Print
 *  
 *  Input :   char *cpMessage -- pointer to vector message.
 *            float *fpVector -- pointer to float array
 *            int     iLength -- length of float array
 *  Output :  void.
 *  ================================================================
 */

void vecPrint( char *cpMessage , float *fpVector , int iLength ) {
int iRow;

    if(cpMessage != (char *) NULL) 
       printf ("\nVECTOR \"%s\" \n\n", cpMessage );
    else 
       printf("\nVECTOR : \"UNTITLED\" \n\n");

    for(iRow = 1; iRow <= iLength; iRow++) 
       printf(" %3d   %16.5e\n", iRow, fpVector[ iRow-1 ]);

}