src/FDM/UIUCModel/uiuc_2DdataFileReader.cpp

   1 /**********************************************************************
   2
   3  FILENAME:     uiuc_2DdataFileReader.cpp
   4
   5 ----------------------------------------------------------------------
   6
   7  DESCRIPTION:  Reads name of data file to be opened and reads data
   8                into appropriate arrays or matrices
   9
  10 ----------------------------------------------------------------------
  11
  12  STATUS:       alpha version
  13
  14 ----------------------------------------------------------------------
  15
  16  REFERENCES:
  17
  18 ----------------------------------------------------------------------
  19
  20  HISTORY:      02/29/2000   initial release
  21                10/25/2001   (RD) Modified so that it recognizes a
  22                             blank line
  23                06/30/2003   (RD) replaced istrstream with istringstream
  24                             to get rid of the annoying warning about
  25                             using the strstream header
  26
  27 ----------------------------------------------------------------------
  28
  29  AUTHOR(S):    Jeff Scott         <jscott@mail.com>
  30                Robert Deters      <rdeters@uiuc.edu>
  31
  32 ----------------------------------------------------------------------
  33
  34  VARIABLES:
  35
  36 ----------------------------------------------------------------------
  37
  38  INPUTS:       -2D file name
  39                -conversion factor for x data
  40                -conversion factor for y data
  41                -conversion factor for z data
  42
  43 ----------------------------------------------------------------------
  44
  45  OUTPUTS:      -2D array of x data for each y case
  46                -1D array of y data
  47                -2D array of z data for each y case
  48                -1D array of max number of x-z data sets for each y case
  49                -max number of y data sets
  50
  51 ----------------------------------------------------------------------
  52
  53  CALLED BY:    uiuc_menu.cpp
  54
  55 ----------------------------------------------------------------------
  56
  57  CALLS TO:     specified 2D data file
  58
  59 ----------------------------------------------------------------------
  60
  61  COPYRIGHT:    (C) 2000 by Michael Selig
  62
  63  This program is free software; you can redistribute it and/or
  64  modify it under the terms of the GNU General Public License
  65  as published by the Free Software Foundation.
  66
  67  This program is distributed in the hope that it will be useful,
  68  but WITHOUT ANY WARRANTY; without even the implied warranty of
  69  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  70  GNU General Public License for more details.
  71
  72  You should have received a copy of the GNU General Public License
  73  along with this program; if not, write to the Free Software
  74  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  75
  76 **********************************************************************/
  77
  78 #include "uiuc_2DdataFileReader.h"
  79
  80
  81 void uiuc_2DdataFileReader( string file_name,
  82                             double x[100][100],
  83                             double y[100],
  84                             double z[100][100],
  85                             int xmax[100],
  86                             int &ymax)
  87 {
  88   ParseFile *matrix;
  89   double token_value1;
  90   double token_value2;
  91   int counter_y = 1, counter_x = 1;
  92   string linetoken1;
  93   string linetoken2;
  94
  95   stack command_list;
  96
  97   /* Read the file and get the list of commands */
  98   matrix = new ParseFile(file_name);
  99   command_list = matrix -> getCommands();
 100
 101   for (LIST command_line = command_list.begin(); command_line!=command_list.end(); ++command_line)
 102     {
 103       linetoken1 = matrix -> getToken(*command_line, 1); // gettoken(string,tokenNo);
 104       linetoken2 = matrix -> getToken(*command_line, 2); // 2 represents token No 2
 105
 106       istringstream token1(linetoken1.c_str());
 107       istringstream token2(linetoken2.c_str());
 108
 109       //reset token_value1 and token_value2 for first if statement
 110       token_value1 = -999;
 111       token_value2 = -999;
 112
 113       token1 >> token_value1;
 114       token2 >> token_value2;
 115
 116       //chenk to see if it is a blank line
 117       if (token_value1==-999 && token_value2==-999);
 118       //check to see if only one value on line (token2 blank)
 119       else if (token_value2 == -999)
 120         {
 121           y[counter_y] = token_value1 * convert_y;
 122
 123           ymax = counter_y;
 124           counter_y++;
 125           counter_x = 1;
 126         }
 127       else
 128         {
 129           x[counter_y-1][counter_x] = token_value1 * convert_x;
 130           z[counter_y-1][counter_x] = token_value2 * convert_z;
 131
 132           xmax[counter_y-1] = counter_x;
 133           counter_x++;
 134         }
 135     }
 136   delete matrix;
 137   return;
 138 }
 139
 140 // end uiuc_2DdataFileReader.cpp