ArrayUtils Namespace Reference

Classes
class	CreateDataImpl

Functions
template<class T >
std::pair< size_t, size_t >	getShape (const T &data)
template<class T >
CreateDataImpl::ReturnType< T >	createData (const T &vec)
PyObject *	createNumpyArray (const std::vector< double > &data)
template<class T >
decltype(auto)	createVector1D (const T &data)
template<class T >
decltype(auto)	createVector2D (const T &data)

Detailed Description

Array and Numpy utility functions getShape, createNumpyArray.

Function Documentation

getShape()

template<class T >

std::pair< size_t, size_t > ArrayUtils::getShape

(

const T &

data

)

Returns shape nrows, ncols of 2D array.

Definition at line 123 of file ArrayUtils.h.

{
    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("Util::getShape() -> Error. "
                                     "Number of elements is different from row to row.");
    return std::make_pair(nrows, ncols);
}

Referenced by Histogram2D::addContent(), ArrayUtils::CreateDataImpl::createDataImpl(), and Histogram2D::initFromShape().

createData()

template<class T >

CreateDataImpl::ReturnType<T> ArrayUtils::createData

(

const T &

vec

)

Creates OutputData array from input vector.

Definition at line 61 of file ArrayUtils.h.

{
    constexpr const int size = CreateDataImpl::nDim<T>::value;
    static_assert(
        size == 1 || size == 2,
        "Error in ArrayUtils::createData: invalid dimensionality or type of the input argument");
    static_assert(std::is_same<CreateDataImpl::ReturnType<T>,
                               decltype(CreateDataImpl::createDataImpl(vec))>::value,
                  "Error in ArrayUtils::createData: invalid return type.");
    return CreateDataImpl::createDataImpl(vec);
}

Referenced by FitObjective::addSimulationAndData(), PyArrayImport::importArrayToOutputData(), and OutputDataReadNumpyTXTStrategy::readOutputData().

createNumpyArray()

PyObject * ArrayUtils::createNumpyArray

(

const std::vector< double > &

data

)

Definition at line 20 of file ArrayUtils.cpp.

{
    const size_t ndim(1);
    npy_int ndim_numpy = ndim;
    npy_intp* ndimsizes_numpy = new npy_intp[ndim];
    ndimsizes_numpy[0] = data.size();
 
    // creating standalone numpy array
    PyObject* pyarray = PyArray_SimpleNew(ndim_numpy, ndimsizes_numpy, NPY_DOUBLE);
    delete[] ndimsizes_numpy;
    if (pyarray == nullptr)
        throw Exceptions::RuntimeErrorException("ExportOutputData() -> Panic in PyArray_SimpleNew");
 
    // getting pointer to data buffer of numpy array
    double* array_buffer = (double*)PyArray_DATA((PyArrayObject*)pyarray);
 
    for (size_t index = 0; index < data.size(); ++index)
        *array_buffer++ = data[index];
 
    return pyarray;
}

Referenced by Histogram1D::getBinCentersNumpy(), Histogram1D::getBinErrorsNumpy(), and Histogram1D::getBinValuesNumpy().

createVector1D()

template<class T >

decltype(auto) ArrayUtils::createVector1D

(

const T &

data

)

Creates 1D vector from OutputData.

createVector2D()

template<class T >

decltype(auto) ArrayUtils::createVector2D

(

const T &

data

)

Creates 2D vector from OutputData.

Referenced by SpectrumUtils::FindPeaks(), and anonymous_namespace{OutputDataWriteStrategy.cpp}::Write2DRepresentation().