lsPointData.hpp

Go to the documentation of this file.

#pragma once
 
#include <lsPreCompileMacros.hpp>
 
#include <array>
#include <cstdint>
#include <vector>
 
#include <lsConcepts.hpp>
 
#include <vcLogger.hpp>
 
namespace viennals {
 
using namespace viennacore;
 
template <class T = double,
          lsConcepts::IsFloatingPoint<T> = lsConcepts::assignable>
class PointData {
public:
  typedef std::vector<T> ScalarDataType;
  typedef std::vector<std::array<T, 3>> VectorDataType;
 
private:
  std::vector<ScalarDataType> scalarData;
  std::vector<std::string> scalarDataLabels;
  std::vector<VectorDataType> vectorData;
  std::vector<std::string> vectorDataLabels;
 
  template <class VectorType, class ReturnType = std::conditional_t<
                                  std::is_const<VectorType>::value,
                                  const typename VectorType::value_type *,
                                  typename VectorType::value_type *>>
  ReturnType indexPointerOrNull(VectorType &v, int index) const {
    if (index >= 0 && index < v.size())
      return &(v[index]);
    else
      Logger::getInstance()
          .addWarning("PointData: Tried to access out of bounds index! "
                      "Returned nullptr instead.")
          .print();
    return nullptr;
  }
 
  template <class DataType>
  void appendTranslateData(DataType &currentData, const DataType &source,
                           const std::vector<unsigned> &indices) {
    currentData.reserve(currentData.size() + indices.size());
    for (unsigned i = 0; i < indices.size(); ++i) {
      currentData.push_back(source[indices[i]]);
    }
  }
 
public:
  void insertNextScalarData(const ScalarDataType &scalars,
                            std::string label = "Scalars") {
    scalarData.push_back(scalars);
    scalarDataLabels.push_back(label);
  }
 
  void insertNextScalarData(ScalarDataType &&scalars,
                            std::string label = "Scalars") {
    scalarData.push_back(std::move(scalars));
    scalarDataLabels.push_back(label);
  }
 
  void insertNextVectorData(const VectorDataType &vectors,
                            std::string label = "Vectors") {
    vectorData.push_back(vectors);
    vectorDataLabels.push_back(label);
  }
 
  void insertNextVectorData(VectorDataType &&vectors,
                            std::string label = "Vectors") {
    vectorData.push_back(std::move(vectors));
    vectorDataLabels.push_back(label);
  }
 
  unsigned getScalarDataSize() const { return scalarData.size(); }
 
  unsigned getVectorDataSize() const { return vectorData.size(); }
 
  ScalarDataType *getScalarData(int index) {
    return indexPointerOrNull(scalarData, index);
  }
 
  const ScalarDataType *getScalarData(int index) const {
    return indexPointerOrNull(scalarData, index);
  }
 
  ScalarDataType *getScalarData(std::string searchLabel,
                                bool noWarning = false) {
    if (int i = getScalarDataIndex(searchLabel); i != -1) {
      return &(scalarData[i]);
    }
    if (!noWarning)
      Logger::getInstance()
          .addWarning("PointData attempted to access scalar data labeled '" +
                      searchLabel +
                      "', which does not exist. Returning nullptr instead.")
          .print();
    return nullptr;
  }
 
  const ScalarDataType *getScalarData(std::string searchLabel,
                                      bool noWarning = false) const {
    if (int i = getScalarDataIndex(searchLabel); i != -1) {
      return &(scalarData[i]);
    }
    if (!noWarning)
      Logger::getInstance()
          .addWarning("PointData attempted to access scalar data labeled '" +
                      searchLabel +
                      "', which does not exist. Returning nullptr instead.")
          .print();
    return nullptr;
  }
 
  int getScalarDataIndex(std::string searchLabel) const {
    for (int i = 0; i < scalarDataLabels.size(); ++i) {
      if (scalarDataLabels[i] == searchLabel) {
        return i;
      }
    }
    return -1;
  }
 
  std::string getScalarDataLabel(int index) const {
    return (index >= 0 && index < scalarDataLabels.size())
               ? scalarDataLabels[index]
               : "";
  }
 
  void setScalarDataLabel(int index, std::string newLabel) {
    scalarDataLabels[index] = newLabel;
  }
 
  void eraseScalarData(int index) {
    scalarData.erase(scalarData.begin() + index);
    scalarDataLabels.erase(scalarDataLabels.begin() + index);
  }
 
  VectorDataType *getVectorData(int index) {
    return indexPointerOrNull(vectorData, index);
  }
 
  const VectorDataType *getVectorData(int index) const {
    return indexPointerOrNull(vectorData, index);
  }
 
  VectorDataType *getVectorData(std::string searchLabel,
                                bool noWarning = false) {
    if (int i = getVectorDataIndex(searchLabel); i != -1) {
      return &(vectorData[i]);
    }
    if (!noWarning)
      Logger::getInstance()
          .addWarning("PointData attempted to access vector data labeled '" +
                      searchLabel +
                      "', which does not exist. Returning nullptr instead.")
          .print();
    return nullptr;
  }
 
  const VectorDataType *getVectorData(std::string searchLabel,
                                      bool noWarning = false) const {
    if (int i = getVectorDataIndex(searchLabel); i != -1) {
      return &(vectorData[i]);
    }
    if (!noWarning)
      Logger::getInstance()
          .addWarning("PointData attempted to access vector data labeled '" +
                      searchLabel +
                      "', which does not exist. Returning nullptr instead.")
          .print();
    return nullptr;
  }
 
  int getVectorDataIndex(std::string searchLabel) const {
    for (int i = 0; i < vectorDataLabels.size(); ++i) {
      if (vectorDataLabels[i] == searchLabel) {
        return i;
      }
    }
    return -1;
  }
 
  std::string getVectorDataLabel(int index) const {
    return (index >= 0 && index < vectorDataLabels.size())
               ? vectorDataLabels[index]
               : "";
  }
 
  void setVectorDataLabel(int index, std::string newLabel) {
    vectorDataLabels[index] = newLabel;
  }
 
  void eraseVectorData(int index) {
    vectorData.erase(vectorData.begin() + index);
    vectorDataLabels.erase(vectorDataLabels.begin() + index);
  }
 
  void append(const PointData &passedData) {
    scalarData.insert(scalarData.end(), passedData.scalarData.begin(),
                      passedData.scalarData.end());
    scalarDataLabels.insert(scalarDataLabels.end(),
                            passedData.scalarDataLabels.begin(),
                            passedData.scalarDataLabels.end());
    vectorData.insert(vectorData.end(), passedData.vectorData.begin(),
                      passedData.vectorData.end());
    vectorDataLabels.insert(vectorDataLabels.end(),
                            passedData.vectorDataLabels.begin(),
                            passedData.vectorDataLabels.end());
  }
 
  void translateFromData(const PointData &source,
                         const std::vector<unsigned> &indices) {
    // scalars
    for (unsigned j = 0; j < source.getScalarDataSize(); ++j) {
      insertNextScalarData(ScalarDataType(), source.getScalarDataLabel(j));
      auto currentData = --scalarData.end();
      appendTranslateData(*currentData, source.scalarData[j], indices);
    }
 
    // vectors
    for (unsigned j = 0; j < source.getVectorDataSize(); ++j) {
      insertNextVectorData(VectorDataType(), source.getVectorDataLabel(j));
      auto currentData = --vectorData.end();
      appendTranslateData(*currentData, source.vectorData[j], indices);
    }
  }
 
  void translateFromMultiData(
      const PointData &source,
      const std::vector<std::vector<unsigned>> &indicesVector) {
    // scalars
    for (unsigned j = 0; j < source.getScalarDataSize(); ++j) {
      insertNextScalarData(ScalarDataType(), source.getScalarDataLabel(j));
      auto currentData = --scalarData.end();
      for (const auto &i : indicesVector) {
        appendTranslateData(*currentData, source.scalarData[j], i);
      }
    }
 
    // vectors
    for (unsigned j = 0; j < source.getVectorDataSize(); ++j) {
      insertNextVectorData(VectorDataType(), source.getVectorDataLabel(j));
      auto currentData = --vectorData.end();
      for (const auto &i : indicesVector) {
        appendTranslateData(*currentData, source.vectorData[j], i);
      }
    }
  }
 
  void clear() {
    scalarData.clear();
    scalarDataLabels.clear();
    vectorData.clear();
    vectorDataLabels.clear();
  }
 
  bool empty() { return scalarData.empty() && vectorData.empty(); }
 
  std::ostream &serialize(std::ostream &stream) {
    // HEADER
    // identifier: "PointData"
    // 4 byte: number of scalar data sets
    // 4 byte: number of vector data sets
    stream << "lsPointData";
    uint32_t numberOfScalarData = scalarData.size();
    uint32_t numberOfVectorData = vectorData.size();
    stream.write(reinterpret_cast<const char *>(&numberOfScalarData),
                 sizeof(uint32_t));
    stream.write(reinterpret_cast<const char *>(&numberOfVectorData),
                 sizeof(uint32_t));
 
    // Scalar Data
    {
      auto labelIt = scalarDataLabels.begin();
      // iterate over all scalar data sets
      for (auto data : scalarData) {
        // write name of scalar data and size of set
        uint32_t sizeOfName = labelIt->length();
        stream.write(reinterpret_cast<char *>(&sizeOfName), sizeof(uint32_t));
        stream << *labelIt;
        uint32_t numberOfValues = data.size();
        stream.write(reinterpret_cast<const char *>(&numberOfValues),
                     sizeof(uint32_t));
        // iterate over scalars in data set
        for (auto value : data) {
          stream.write(reinterpret_cast<const char *>(&value),
                       sizeof(typename ScalarDataType::value_type));
        }
        ++labelIt;
      }
    }
 
    // Vector Data
    {
      auto labelIt = vectorDataLabels.begin();
      // iterate over all vector data sets
      for (auto data : vectorData) {
        // write name of vector data and size of set
        uint32_t sizeOfName = labelIt->length();
        stream.write(reinterpret_cast<char *>(&sizeOfName), sizeof(uint32_t));
        stream << *labelIt;
        uint32_t numberOfVectors = data.size();
        stream.write(reinterpret_cast<const char *>(&numberOfVectors),
                     sizeof(uint32_t));
        // iterate over vectors in data set
        for (auto vector : data) {
          // over values in vector
          for (auto value : vector) {
            stream.write(
                reinterpret_cast<const char *>(&value),
                sizeof(typename VectorDataType::value_type::value_type));
          }
        }
        ++labelIt;
      }
    }
 
    return stream;
  }
 
  std::istream &deserialize(std::istream &stream) {
    char identifier[11];
    stream.read(identifier, 11);
    if (std::string(identifier).compare(0, 11, "lsPointData")) {
      Logger::getInstance()
          .addWarning("Reading PointData from stream failed. Header could "
                      "not be found.")
          .print();
      return stream;
    }
 
    // read number of data sets
    unsigned numberOfScalarData = 0;
    unsigned numberOfVectorData = 0;
    stream.read(reinterpret_cast<char *>(&numberOfScalarData),
                sizeof(uint32_t));
    stream.read(reinterpret_cast<char *>(&numberOfVectorData),
                sizeof(uint32_t));
 
    // read scalar data
    for (unsigned i = 0; i < numberOfScalarData; ++i) {
      uint32_t sizeOfName;
      stream.read(reinterpret_cast<char *>(&sizeOfName), sizeof(uint32_t));
      std::vector<char> dataLabel(sizeOfName);
      stream.read(dataLabel.data(), sizeOfName);
      uint32_t numberOfValues;
      stream.read(reinterpret_cast<char *>(&numberOfValues), sizeof(uint32_t));
      ScalarDataType scalarData;
      scalarData.resize(numberOfValues);
      // read all scalar values into the data
      for (auto &value : scalarData) {
        stream.read(reinterpret_cast<char *>(&value),
                    sizeof(typename ScalarDataType::value_type));
      }
      // now add this scalar data to current PointData
      insertNextScalarData(scalarData,
                           std::string(dataLabel.begin(), dataLabel.end()));
    }
 
    // read vector data
    for (unsigned i = 0; i < numberOfVectorData; ++i) {
      uint32_t sizeOfName;
      stream.read(reinterpret_cast<char *>(&sizeOfName), sizeof(uint32_t));
      std::vector<char> dataLabel(sizeOfName);
      stream.read(dataLabel.data(), sizeOfName);
      uint32_t numberOfValues;
      stream.read(reinterpret_cast<char *>(&numberOfValues), sizeof(uint32_t));
      VectorDataType vectorData;
      vectorData.resize(numberOfValues);
      // read all vector values into the vector data
      for (auto &vector : vectorData) {
        for (auto &value : vector) {
          stream.read(reinterpret_cast<char *>(&value),
                      sizeof(typename VectorDataType::value_type::value_type));
        }
      }
      // now add this scalar data to current PointData
      insertNextVectorData(vectorData,
                           std::string(dataLabel.begin(), dataLabel.end()));
    }
 
    return stream;
  }
};
 
// add all template specialisations for this class
PRECOMPILE_PRECISION(PointData);
 
} // namespace viennals