DGF.tpp

/*---------------------------------------------------------------------------*\
 *
 *  bitpit
 *
 *  Copyright (C) 2015-2021 OPTIMAD engineering Srl
 *
 *  -------------------------------------------------------------------------
 *  License
 *  This file is part of bitpit.
 *
 *  bitpit is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License v3 (LGPL)
 *  as published by the Free Software Foundation.
 *
 *  bitpit is distributed in the hope that it will be useful, but WITHOUT
 *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 *  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
 *  License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with bitpit. If not, see <http://www.gnu.org/licenses/>.
 *
\*---------------------------------------------------------------------------*/
 
// ========================================================================== //
//                              DGF IO FUNCTIONS                              //
//                                                                            //
// I/O functions for .dgf file formats                                        //
// ========================================================================== //
// INTFO                                                                      //
// ========================================================================== //
// Author  : Alessandro Alaia                                                 //
// Version : v3.0                                                             //
//                                                                            //
// All rights reserved.                                                       //
// ========================================================================== //
 
// ========================================================================== //
// TEMPLATE IMPLEMENTATIONS                                                   //
// ========================================================================== //
 
// class DGFObj templated methods ========================================== //
 
 
// -------------------------------------------------------------------------- // 
template< typename T, typename ... T2 >
void DGFObj::loadVData(
    const std::string           &data_name,
    int                         &n,
    std::vector< T >            &data,
    T2                          &... others
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
// none
 
// Counters
// none
 
// ========================================================================== //
// OPEN INPUT STREAM                                                          //
// ========================================================================== //
open("in");
 
// ========================================================================== //
// READ DATA SET                                                              //
// ========================================================================== //
err = dgf::readVertexData(ifile_handle, n, data, data_name);
 
// ========================================================================== //
// ITERATIVELY READ OTHER DATASET                                             //
// ========================================================================== //
loadVData(others ...);
 
return; };
 
// -------------------------------------------------------------------------- //
template< typename T, typename ... T2 >
void DGFObj::loadSData(
    const std::string           &data_name,
    int                         &n,
    std::vector< T >            &data,
    T2                          &... others
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
// none
 
// Counters
// none
 
// ========================================================================== //
// OPEN INPUT STREAM                                                          //
// ========================================================================== //
open("in");
 
// ========================================================================== //
// READ DATA SET                                                              //
// ========================================================================== //
err = dgf::readSimplexData(ifile_handle, n, data, data_name);
 
// ========================================================================== //
// ITERATIVELY READ OTHER DATASET                                             //
// ========================================================================== //
loadSData(others ...);
 
return; };
 
// -------------------------------------------------------------------------- // 
template< typename T, typename ... T2 >
void DGFObj::appendVData(
    const std::string           &data_name,
    int                         &n,
    std::vector< T >            &data,
    T2                          &... others
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
// none
 
// Counters
// none
 
// ========================================================================== //
// OPEN OUTPUT STREAM                                                         //
// ========================================================================== //
open("app");
 
// ========================================================================== //
// SAVE DATA SET                                                              //
// ========================================================================== //
err = dgf::writeVertexData(ofile_handle, n, data, data_name);
 
// ========================================================================== //
// RECURSIVELY SAVE OTHERS DATASETS                                           //
// ========================================================================== //
appendVData(others ...);
 
return; }
 
// -------------------------------------------------------------------------- // 
template< typename T, typename ... T2 >
void DGFObj::appendSData(
    const std::string           &data_name,
    int                         &n,
    std::vector< T >            &data,
    T2                          &... others
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
// none
 
// Counters
// none
 
// ========================================================================== //
// OPEN OUTPUT STREAM                                                         //
// ========================================================================== //
open("app");
 
// ========================================================================== //
// SAVE DATA SET                                                              //
// ========================================================================== //
err = dgf::writeSimplexData(ofile_handle, n, data, data_name);
 
// ========================================================================== //
// RECURSIVELY SAVE OTHERS DATASETS                                           //
// ========================================================================== //
appendSData(others ...);
 
return; }

 
// Input routines =========================================================== //
 
namespace dgf {
 
// -------------------------------------------------------------------------- //
template< typename T >
unsigned int readData(
    std::ifstream               &file_handle,
    int                         &N,
    std::vector< T >            &Data
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
bool                     check;
std::streamoff           start_pos;
std::string              word, line;
std::stringstream        sline;
 
// Counters
int                 n = 0;
 
// ========================================================================== //
// CHECK INPUT STREAM                                                         //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// SCAN DATA BLOCK                                                            //
// ========================================================================== //
file_handle.clear();
start_pos = file_handle.tellg();
dgf::scanData(file_handle, n);
file_handle.clear();
file_handle.seekg(start_pos);
 
// ========================================================================== //
// RESIZE INPUT DATA                                                          //
// ========================================================================== //
Data.resize(n+N);
 
// ========================================================================== //
// READ DATA                                                                  //
// ========================================================================== //
n = 0;
check = true;
while (!file_handle.eof() && check) {
 
    // Get current line
    start_pos = file_handle.tellg();
    std::getline(file_handle, line);
    line = utils::string::trim(line);
    sline.clear();
    sline.str(line);
 
    // Read data
    if (sline >> word) {
         check =  ((word.compare("#") != 0)
              &&  (word.compare("VERTEX") != 0)
              &&  (word.compare("SIMPLEX") != 0)
              &&  (word.compare("VERTEXDATA") != 0)
              &&  (word.compare("SIMPLEXDATA") != 0));
        if (check) {
            sline.seekg(0);
            sline >> Data[N+n];
            n++;
        }
    }
} //next line
if (word.compare("#") != 0) {
    file_handle.clear();
    file_handle.seekg(start_pos);
}
 
// Update counters
N+=n;
 
return(0); }
 
// -------------------------------------------------------------------------- // 
template <typename T>
unsigned int readVertexData(
    std::ifstream               &file_handle,
    int                         &n,
    std::vector< T >            &data,
    const std::string           &data_name
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
bool            check = false;
long int        current_pos, start_pos;
std::string          header, line, word;
std::stringstream    sline;
 
// Counters
 
// ========================================================================== //
// CHECK STREAM STATUS                                                        //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// PARAMETERS                                                                 //
// ========================================================================== //
sline << "VERTEXDATA " << data_name << std::endl;
header = sline.str();
header = utils::string::trim(header);
 
// ========================================================================== //
// SCAN DGF FILE LOOKING FOR DATASET WITH SPECIFIED NAME                      //
// ========================================================================== //
file_handle.clear();
start_pos = file_handle.tellg();
current_pos = -1;
while (start_pos != current_pos) {
 
    // Get current line
    std::getline(file_handle, line);
    line = utils::string::trim(line);
    sline.clear();
    sline.str(line);
 
    // Check eof
    if (file_handle.eof()) {
        file_handle.clear();
        file_handle.seekg(0);
    }
    current_pos = file_handle.tellg();
 
    // Look for keyword
    if ((sline >> word) && (word.compare("VERTEXDATA") == 0)) {
        if (data_name.compare("") == 0) {
            start_pos = current_pos;
            check = true;
        }
        else {
            if (line.compare(header) == 0) {
                start_pos = current_pos;
                check = true;
            }
        }
    }
 
    
} //next line
 
// ========================================================================== //
// LOAD DATA                                                                  //
// ========================================================================== //
if (check) {
    dgf::readData(file_handle, n, data);
}
 
return(0); }
 
// -------------------------------------------------------------------------- //
template <typename T>
unsigned int readSimplexData(
    std::ifstream               &file_handle,
    int                         &n,
    std::vector< T >            &data,
    const std::string           &data_name
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
bool                 check = false;
std::streamoff       current_pos, start_pos;
std::string          header, line, word;
std::stringstream    sline;
 
// Counters
 
// ========================================================================== //
// CHECK STREAM STATUS                                                        //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// PARAMETERS                                                                 //
// ========================================================================== //
sline << "SIMPLEXDATA " << data_name << std::endl;
header = sline.str();
header = utils::string::trim(header);
 
// ========================================================================== //
// SCAN DGF FILE LOOKING FOR DATASET WITH SPECIFIED NAME                      //
// ========================================================================== //
file_handle.clear();
start_pos = file_handle.tellg();
current_pos = -1;
while (start_pos != current_pos) {
 
    // Get current line
    std::getline(file_handle, line);
    line = utils::string::trim(line);
    sline.clear();
    sline.str(line);
 
    // Check eof
    if (file_handle.eof()) {
        file_handle.clear();
        file_handle.seekg(0);
    }
    current_pos = file_handle.tellg();
 
    // Look for keyword
    if ((sline >> word) && (word.compare("SIMPLEXDATA") == 0)) {
        if (data_name.compare("") == 0) {
            start_pos = current_pos;
            check = true;
        }
        else {
            if (line.compare(header) == 0) {
                start_pos = current_pos;
                check = true;
            }
        }
    }
 
    
} //next line
 
// ========================================================================== //
// LOAD DATA                                                                  //
// ========================================================================== //
if (check) {
    dgf::readData(file_handle, n, data);
}
 
return(0); }
 
// Output routines ========================================================== //
 
// -------------------------------------------------------------------------- // 
template < typename T >
unsigned int writeData(
    std::ofstream               &file_handle,
    int                         &N,
    std::vector< T >            &Data
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
// none
 
// Counters
int             i;
 
// ========================================================================== //
// CHECK STREAM STATUS                                                        //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// EXPORT DATA                                                                //
// ========================================================================== //
for (i = 0; i < N; i++) {
    file_handle << Data[i] << std::endl;
} //next i
file_handle << "#" << std::endl << std::endl;
 
return(0); }
 
// -------------------------------------------------------------------------- // 
template < typename T >
unsigned int writeVertexData(
    std::ofstream               &file_handle,
    int                         &N,
    std::vector< T >            &Data,
    const std::string           &Data_name
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
unsigned int        err = 0;
std::stringstream        sheader;
std::string              dataset;
std::string              header;
 
// Counters
// None
 
// ========================================================================== //
// CHECK STREAM STATUS                                                        //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// EXPORT DATA                                                                //
// ========================================================================== //
 
// Data header -------------------------------------------------------------- //
dataset = Data_name;
dataset = utils::string::trim(dataset);
sheader << "VERTEXDATA " << dataset << std::endl;
header = sheader.str();
header = utils::string::trim(header);
file_handle << header << std::endl;
 
// Export data -------------------------------------------------------------- //
err = dgf::writeData(file_handle, N, Data);
 
return(err); };
 
// -------------------------------------------------------------------------- // 
template < typename T >
unsigned int writeSimplexData(
    std::ofstream               &file_handle,
    int                         &N,
    std::vector< T >            &Data,
    const std::string           &Data_name
) {
 
// ========================================================================== //
// VARIABLES DECLARATION                                                      //
// ========================================================================== //
 
// Local variables
unsigned int        err = 0;
std::stringstream        sheader;
std::string              dataset;
std::string              header;
 
// Counters
// None
 
// ========================================================================== //
// CHECK STREAM STATUS                                                        //
// ========================================================================== //
if (!file_handle.good()) { return(1); }
 
// ========================================================================== //
// EXPORT DATA                                                                //
// ========================================================================== //
 
// Data header -------------------------------------------------------------- //
dataset = Data_name;
dataset = utils::string::trim(dataset);
sheader << "SIMPLEXDATA " << dataset << std::endl;
header = sheader.str();
header = utils::string::trim(header);
file_handle << header << std::endl;
 
// Export data -------------------------------------------------------------- //
err = dgf::writeData(file_handle, N, Data);
 
return(err); };
 
}