ReadWriteNumpyTXT Class Reference

Description

Class for reading and writing Datafield from simple ASCII file with the layout as in numpy.savetxt.

Definition at line 30 of file ReadWriteNumpyTXT.h.

Public Member Functions
Datafield *	readDatafield (std::istream &input_stream)
void	writeDatafield (const Datafield &data, std::ostream &output_stream)

Static Private Member Functions
static double	ignoreDenormalized (double value)
static void	write1DRepresentation (const Datafield &data, std::ostream &output_stream)
static void	write2DRepresentation (const Datafield &data, std::ostream &output_stream)

Member Function Documentation

ignoreDenormalized()

double ReadWriteNumpyTXT::ignoreDenormalized ( double  value )

staticprivate

Definition at line 134 of file ReadWriteNumpyTXT.cpp.

{
    return (std::fpclassify(value) == FP_SUBNORMAL) ? 0.0 : value;
}

Referenced by write1DRepresentation(), and write2DRepresentation().

readDatafield()

Datafield * ReadWriteNumpyTXT::readDatafield

(

std::istream &

input_stream

)

Definition at line 33 of file ReadWriteNumpyTXT.cpp.

{
    std::string line;
    std::vector<std::vector<double>> data;
 
    // Read numbers from input stream:
    while (std::getline(input_stream, line)) {
        line = BaseUtils::String::trim(line);
        if (line.empty() || !isDoubleStartChar(line[0]))
            continue;
 
        try {
            std::vector<double> dataInRow = DataUtils::Format::parse_doubles(line);
            data.push_back(dataInRow);
        } catch (...) {
            continue;
        }
    }
 
    // validating
    size_t nrows = data.size();
    size_t ncols(0);
    if (nrows)
        ncols = data[0].size();
 
    if (ncols == 0)
        throw std::runtime_error("ReadNumpyTXTStrategy::readDatafield() -> Error. "
                                 "Can't parse file");
 
    for (size_t row = 0; row < nrows; row++) {
        if (data[row].size() != ncols)
            throw std::runtime_error("ReadNumpyTXTStrategy::readDatafield() -> Error. "
                                     "Number of elements is different from row to row.");
    }
 
    if (nrows < 2)
        return DataUtils::Array::createPField1D(data[0]).release();
 
    if (ncols < 2) {
        const size_t size = data.size();
        std::vector<double> vector1d(size);
        for (size_t i = 0; i < size; ++i)
            vector1d[i] = data[i][0];
        return DataUtils::Array::createPField1D(vector1d).release();
    }
 
    return DataUtils::Array::createPField2D(data).release();
}

References DataUtils::Array::createPField1D(), DataUtils::Array::createPField2D(), DataUtils::Format::parse_doubles(), and BaseUtils::String::trim().

Referenced by IOFactory::readDatafield().

Here is the call graph for this function:

write1DRepresentation()

void ReadWriteNumpyTXT::write1DRepresentation ( const Datafield &  data, std::ostream &  output_stream  )

staticprivate

Definition at line 101 of file ReadWriteNumpyTXT.cpp.

{
    output_stream << "# coordinates         intensities" << std::endl;
    output_stream.imbue(std::locale::classic());
    output_stream << std::scientific << std::setprecision(12);
 
    const std::vector<double> axis_values = data.axis(0).binCenters();
 
    // printing coordinate and associated intensity
    for (size_t i = 0, nrows = axis_values.size(); i < nrows; ++i)
        output_stream << axis_values[i] << "    " << ignoreDenormalized(data[i]) << std::endl;
}

References Datafield::axis(), IAxis::binCenters(), ignoreDenormalized(), and BaseUtils::String::scientific().

Referenced by writeDatafield().

Here is the call graph for this function:

write2DRepresentation()

void ReadWriteNumpyTXT::write2DRepresentation ( const Datafield &  data, std::ostream &  output_stream  )

staticprivate

Definition at line 114 of file ReadWriteNumpyTXT.cpp.

{
    const size_t nrows = data.axis(1).size();
    const size_t ncols = data.axis(0).size();
 
    output_stream << "# [nrows=" << nrows << ", ncols=" << ncols << "]" << std::endl;
 
    std::vector<std::vector<double>> dataArray = DataUtils::Array::createVector2D(data);
    output_stream.imbue(std::locale::classic());
    output_stream << std::scientific << std::setprecision(12);
 
    for (size_t i = 0; i < nrows; i++) {
        for (size_t j = 0; j < ncols; j++) {
            double z_value = dataArray[i][j];
            output_stream << ignoreDenormalized(z_value) << "    ";
        }
        output_stream << std::endl;
    }
}

References Datafield::axis(), DataUtils::Array::createVector2D(), ignoreDenormalized(), BaseUtils::String::scientific(), and IAxis::size().

Referenced by writeDatafield().

Here is the call graph for this function:

writeDatafield()

void ReadWriteNumpyTXT::writeDatafield	(	const Datafield &	data,
		std::ostream &	output_stream
	)

Definition at line 82 of file ReadWriteNumpyTXT.cpp.

{
    output_stream << "# BornAgain Intensity Data" << std::endl;
    output_stream << "# Simple array suitable for numpy, matlab etc." << std::endl;
 
    const size_t dim = data.rank();
    switch (dim) {
    case 1:
        write1DRepresentation(data, output_stream);
        break;
    case 2:
        write2DRepresentation(data, output_stream);
        break;
    default:
        throw std::runtime_error("Error in DatafieldWriteNumpyTXTStrategy::writeDatafield: data "
                                 "of unsupported dimensions");
    }
}

References Datafield::rank(), write1DRepresentation(), and write2DRepresentation().

Referenced by IOFactory::writeDatafield().

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• /home/build/builds/o5h8MZZm/0/mlz/bornagain/Device/IO/ReadWriteNumpyTXT.h
• /home/build/builds/o5h8MZZm/0/mlz/bornagain/Device/IO/ReadWriteNumpyTXT.cpp