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., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
  75  USA or view http://www.gnu.org/copyleft/gpl.html.
  76
  77 **********************************************************************/
  78
  79 #include "uiuc_2DdataFileReader.h"
  80
  81
  82 void uiuc_2DdataFileReader( string file_name,
  83                             double x[100][100],
  84                             double y[100],
  85                             double z[100][100],
  86                             int xmax[100],
  87                             int &ymax)
  88 {
  89   ParseFile *matrix;
  90   double token_value1;
  91   double token_value2;
  92   int counter_y = 1, counter_x = 1;
  93   string linetoken1;
  94   string linetoken2;
  95
  96   stack command_list;
  97
  98   /* Read the file and get the list of commands */
  99   matrix = new ParseFile(file_name);
 100   command_list = matrix -> getCommands();
 101
 102   for (LIST command_line = command_list.begin(); command_line!=command_list.end(); ++command_line)
 103     {
 104       linetoken1 = matrix -> getToken(*command_line, 1); // gettoken(string,tokenNo);
 105       linetoken2 = matrix -> getToken(*command_line, 2); // 2 represents token No 2
 106
 107       istringstream token1(linetoken1.c_str());
 108       istringstream token2(linetoken2.c_str());
 109
 110       //reset token_value1 and token_value2 for first if statement
 111       token_value1 = -999;
 112       token_value2 = -999;
 113
 114       token1 >> token_value1;
 115       token2 >> token_value2;
 116
 117       //chenk to see if it is a blank line
 118       if (token_value1==-999 && token_value2==-999);
 119       //check to see if only one value on line (token2 blank)
 120       else if (token_value2 == -999)
 121         {
 122           y[counter_y] = token_value1 * convert_y;
 123
 124           ymax = counter_y;
 125           counter_y++;
 126           counter_x = 1;
 127         }
 128       else
 129         {
 130           x[counter_y-1][counter_x] = token_value1 * convert_x;
 131           z[counter_y-1][counter_x] = token_value2 * convert_z;
 132
 133           xmax[counter_y-1] = counter_x;
 134           counter_x++;
 135         }
 136     }
 137   return;
 138 }
 139
 140 // end uiuc_2DdataFileReader.cpp