Line data    Source code

       1              : // This file is part of the Mamico project. For conditions of distribution
       2              : // and use, please see the copyright notice in Mamico's main folder, or at
       3              : // www5.in.tum.de/mamico
       4              : 
       5              : /** Implementation of ReadFromFile.h
       6              :  *  @author Felix Maurer
       7              :  */
       8              : 
       9              : // Member functions of ReadFromFile.h:
      10            0 : template <unsigned int dim> void coupling::filtering::ReadFromFile<dim>::operator()() {
      11              : #ifdef DEBUG_READ_FROM_FILE
      12              :   std::cout << "             RFF: Now reading input data from: " << _location;
      13              : #endif
      14              : 
      15            0 :   _file.open(_location);
      16              : 
      17              :   // Skip all lines of previous iterations
      18            0 :   std::string skipstr;
      19            0 :   for (unsigned int s = 0; s < _iteration * ((coupling::filtering::FilterInterface<dim>::_inputCells.size()) + 1); s++)
      20            0 :     std::getline(_file, skipstr);
      21              : 
      22              :   // TODO: check if number of lines match number of cells
      23              : 
      24            0 :   for (unsigned int ic_index = 0; ic_index < coupling::filtering::FilterInterface<dim>::_inputCells.size(); ic_index++) {
      25              :     // get next input line
      26            0 :     std::string linestr;
      27            0 :     std::getline(_file, linestr);
      28            0 :     std::istringstream lineiss(linestr);
      29              : 
      30              :     // store all values from line
      31            0 :     std::vector<double> doublesFromLine;
      32              : 
      33              :     // interpret substrings as doubles
      34            0 :     std::string d;
      35            0 :     while (std::getline(lineiss, d, ';')) {
      36            0 :       if (d != " ") { // ignore what comes after last ';'
      37            0 :         doublesFromLine.push_back(std::stod(d));
      38              : 
      39              : #ifdef DEBUG_READ_FROM_FILE
      40              :         if (ic_index == 0)
      41              :           std::cout << "             RFF: New double read for cell 0: " << doublesFromLine.back() << std::endl;
      42              : #endif
      43              :       }
      44              :     }
      45              : 
      46              :     // TODO: check if number of doubles match dim*_vectorSetters.size() +
      47              :     // _scalarSetters.size()
      48              : 
      49              :     // iteration, indexes. Adjust this if write-to-file's output changes!
      50            0 :     const unsigned int const_offset = 1;
      51              : 
      52              :     // increases for each double read from line
      53            0 :     unsigned int incr_offset = 0;
      54              : 
      55              :     // Apply scalar setters
      56            0 :     for (const auto scalarProperty : coupling::filtering::FilterInterface<dim>::_scalarAccessFunctionPairs) {
      57            0 :       (coupling::filtering::FilterInterface<dim>::_outputCells[ic_index]->*scalarProperty.set)(doublesFromLine[const_offset + incr_offset]);
      58              : 
      59              : #ifdef DEBUG_READ_FROM_FILE
      60              :       if (ic_index == 0)
      61              :         std::cout << "               RFF: Scalar sent to cell 0: " << doublesFromLine[0 /*const_offset + incr_offset*/] << std::endl;
      62              : #endif
      63              : 
      64            0 :       incr_offset++;
      65              :     }
      66              :     // Apply vector setters
      67            0 :     for (const auto vectorProperty : coupling::filtering::FilterInterface<dim>::_vectorAccessFunctionPairs) {
      68            0 :       tarch::la::Vector<dim, double> tempVec;
      69            0 :       for (unsigned int d = 0; d < dim; d++) {
      70            0 :         tempVec[d] = doublesFromLine[const_offset + incr_offset + d];
      71              :       }
      72              : 
      73              : #ifdef DEBUG_READ_FROM_FILE
      74              :       if (ic_index == 0)
      75              :         std::cout << "               RFF: Vector sent to cell 0: " << tempVec << std::endl;
      76              : #endif
      77              : 
      78            0 :       (coupling::filtering::FilterInterface<dim>::_outputCells[ic_index]->*vectorProperty.set)(tempVec);
      79            0 :       incr_offset += dim;
      80              :     }
      81              :   }
      82              : 
      83            0 :   _file.close();
      84            0 :   _iteration++;
      85              : 
      86              : #ifdef DEBUG_READ_FROM_FILE
      87              :   std::cout << " ... done" << std::endl;
      88              : #endif
      89            0 : }