ReadWriteINT.cpp

Go to the documentation of this file.

//  ************************************************************************************************
//
//  BornAgain: simulate and fit reflection and scattering
//
//! @file      Device/IO/ReadWriteINT.cpp
//! @brief     Implements class ReadWriteINT.
//!
//! @homepage  http://www.bornagainproject.org
//! @license   GNU General Public License v3 or higher (see COPYING)
//! @copyright Forschungszentrum Jülich GmbH 2018
//! @authors   Scientific Computing Group at MLZ (see CITATION, AUTHORS)
//
//  ************************************************************************************************
 
#include "Device/IO/ReadWriteINT.h"
#include "Base/Axis/Bin.h"
#include "Base/Axis/IAxis.h"
#include "Base/Util/StringUtils.h"
#include "Device/Data/ArrayUtils.h"
#include "Device/IO/DataFormatUtils.h"
#include <iostream>
 
Datafield* ReadWriteINT::readDatafield(std::istream& input_stream)
{
    std::string line;
 
    std::vector<IAxis*> axes;
    while (std::getline(input_stream, line)) {
        line = BaseUtils::String::trim(line);
 
        if (line.find("axis") != std::string::npos)
            axes.emplace_back(DataUtils::Format::createAxis(input_stream));
 
        if (line.find("data") != std::string::npos)
            break;
    }
    // std::cout << "Read " << axes.size() << " axes" << std::endl;
    auto* result = new Datafield(axes);
    DataUtils::Format::fillDatafield(result, input_stream);
    return result;
}
 
void ReadWriteINT::writeDatafield(const Datafield& data, std::ostream& output_stream)
{
    output_stream << "# BornAgain Intensity Data\n\n";
 
    for (size_t i = 0; i < data.rank(); ++i) {
        std::string axis_name = std::string("axis") + std::to_string(i);
        std::unique_ptr<IAxis> axis(data.axis(i).clone());
        axis->setAxisName(axis_name);
        output_stream << std::endl;
        output_stream << "# axis-" << i << "\n";
        output_stream << (*axis) << "\n";
    }
    size_t n_columns = data.axis(data.rank() - 1).size();
 
    output_stream << "\n# data\n";
    writeDatafieldDoubles(data, output_stream, n_columns);
    output_stream << std::endl;
}
 
void ReadWriteINT::writeDatafieldDoubles(const Datafield& data, std::ostream& output_stream,
                                         size_t n_columns)
{
    output_stream.imbue(std::locale::classic());
    output_stream << std::scientific << std::setprecision(12);
    size_t ncol = 0;
    for (size_t i = 0; i < data.size(); ++i) {
        ncol++;
        output_stream << ignoreDenormalized(data[i]) << "    ";
        if (ncol == n_columns) {
            output_stream << std::endl;
            ncol = 0;
        }
    }
}
 
double ReadWriteINT::ignoreDenormalized(double value)
{
    return (std::fpclassify(value) == FP_SUBNORMAL) ? 0.0 : value;
}