DataUtils::Array Namespace Reference

Description

Array utility functions getShape.

Functions
std::unique_ptr< Datafield >	createPField1D (const std::vector< double > &vec)
std::unique_ptr< Datafield >	createPField2D (const std::vector< std::vector< double >> &vec)
std::vector< double >	createVector1D (const Datafield &data)
	Creates 1D vector from Datafield. More...
std::vector< std::vector< double > >	createVector2D (const Datafield &data)
	Creates 2D vector from Datafield. More...
std::pair< size_t, size_t >	getShape (const std::vector< std::vector< double >> &data)
	Returns shape nrows, ncols of 2D array. More...

Function Documentation

createPField1D()

std::unique_ptr< Datafield > DataUtils::Array::createPField1D

(

const std::vector< double > &

vec

)

Definition at line 32 of file ArrayUtils.cpp.

{
    const size_t N = vec.size();
    return std::unique_ptr<Datafield>(new Datafield{{
                new FixedBinAxis("axis0", N, 0.0, (double)N)}, vec});
}

References N.

Referenced by FitObjective::addSimulationAndData(), DataUtils::Data::importArrayToDatafield(), and ReadWriteNumpyTXT::readDatafield().

createPField2D()

std::unique_ptr< Datafield > DataUtils::Array::createPField2D

(

const std::vector< std::vector< double >> &

vec

)

Definition at line 40 of file ArrayUtils.cpp.

{
    auto shape = DataUtils::Array::getShape(vec);
    const size_t nrows = shape.first;
    const size_t ncols = shape.second;
 
    if (nrows == 0 || ncols == 0)
        throw std::runtime_error(
            "Error in DataUtils::Array::createDataImpl: input argument contains empty dimensions");
 
    auto* frame = new Frame({new FixedBinAxis("axis0", ncols, 0.0, (double)ncols),
                             new FixedBinAxis("axis1", nrows, 0.0, (double)nrows)});
 
    std::vector<double> out(frame->size());
    for (size_t row = 0; row < nrows; ++row) {
        for (size_t col = 0; col < ncols; ++col) {
            size_t iout = nrows - row - 1 + col * nrows;
            out[iout] = vec[row][col];
        }
    }
 
    return std::make_unique<Datafield>(frame, out);
}

References getShape().

Referenced by FitObjective::addSimulationAndData(), DataUtils::Data::importArrayToDatafield(), and ReadWriteNumpyTXT::readDatafield().

Here is the call graph for this function:

createVector1D()

std::vector< double > DataUtils::Array::createVector1D

(

const Datafield &

data

)

Creates 1D vector from Datafield.

Definition at line 64 of file ArrayUtils.cpp.

{
    if (data.rank() != 1)
        throw std::runtime_error("DataUtils::Array::createVector1D() -> Error. Not 1D data.");
 
    std::vector<double> result = data.flatVector();
    return result;
}

References Datafield::flatVector(), and Datafield::rank().

Here is the call graph for this function:

createVector2D()

std::vector< std::vector< double > > DataUtils::Array::createVector2D

(

const Datafield &

data

)

Creates 2D vector from Datafield.

Definition at line 73 of file ArrayUtils.cpp.

{
    std::vector<std::vector<double>> result;
 
    const size_t nrows = data.axis(1).size();
    const size_t ncols = data.axis(0).size();
 
    result.resize(nrows);
 
    for (size_t row = 0; row < nrows; ++row) {
        result[row].resize(ncols, 0.0);
        for (size_t col = 0; col < ncols; ++col) {
            size_t globalbin = nrows - row - 1 + col * nrows;
            result[row][col] = data[globalbin];
        }
    }
 
    return result;
}

References Datafield::axis(), and IAxis::size().

Referenced by DataUtils::Data::FindPeaks(), and ReadWriteNumpyTXT::write2DRepresentation().

Here is the call graph for this function:

getShape()

std::pair< size_t, size_t > DataUtils::Array::getShape

(

const std::vector< std::vector< double >> &

data

)

Returns shape nrows, ncols of 2D array.

Definition at line 20 of file ArrayUtils.cpp.

{
    size_t nrows = data.size();
    size_t ncols(0);
    if (nrows)
        ncols = data[0].size();
    for (size_t row = 0; row < nrows; row++)
        if (data[row].size() != ncols)
            throw std::runtime_error("Invalid data: not rectangular");
    return std::make_pair(nrows, ncols);
}

Referenced by createPField2D().