VectorColumn.h

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2013 ISIS Rutherford Appleton Laboratory UKRI,
//   NScD Oak Ridge National Laboratory, European Spallation Source,
//   Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
// SPDX - License - Identifier: GPL - 3.0 +
#pragma once
 
#include "MantidAPI/Column.h"
#include "MantidDataObjects/DllConfig.h"
#include "MantidKernel/StringTokenizer.h"
 
#include <boost/algorithm/string/join.hpp>
#include <boost/algorithm/string/trim.hpp>
#include <boost/lexical_cast.hpp>
 
#include <cmath>
namespace Mantid {
namespace DataObjects {
 
template <class Type> class MANTID_DATAOBJECTS_DLL VectorColumn : public API::Column {
public:
  VectorColumn() { m_type = typeName(); }
 
  size_t size() const override { return m_data.size(); }
 
  const std::type_info &get_type_info() const override { return typeid(std::vector<Type>); }
 
  const std::type_info &get_pointer_type_info() const override { return typeid(std::vector<Type> *); }
 
  void print(size_t index, std::ostream &s) const override {
    const std::vector<Type> &values = m_data.at(index);
 
    auto it = values.begin();
 
    if (it != values.end()) {
      s << *it;
      ++it;
    }
 
    for (; it != values.end(); ++it) {
      s << ',';
      s << *it;
    }
  }
 
  void read(size_t index, const std::string &text) override {
    std::vector<Type> newValues;
    Mantid::Kernel::StringTokenizer elements(text, ",", Mantid::Kernel::StringTokenizer::TOK_TRIM);
 
    for (const auto &element : elements) {
      try {
        newValues.emplace_back(boost::lexical_cast<Type>(element));
      } catch (boost::bad_lexical_cast &) {
        throw std::invalid_argument("Unable to convert one of the elements: " + element);
      }
    }
 
    m_data.at(index) = newValues;
  }
 
  void read(const size_t index, std::istringstream &in) override {
    std::string s;
    in >> s;
    read(index, s);
  }
 
  bool isBool() const override { return false; }
 
  bool isNumber() const override { return false; }
 
  long int sizeOfData() const override {
    long int dataSize(0);
 
    for (auto elemIt = m_data.begin(); elemIt != m_data.end(); ++elemIt) {
      dataSize += static_cast<long int>(elemIt->size() * sizeof(Type));
    }
 
    return dataSize;
  }
 
  VectorColumn *clone() const override {
    auto newColumn = new VectorColumn<Type>();
    newColumn->m_data = m_data;
    newColumn->setName(m_name);
    newColumn->setPlotType(m_plotType);
    return newColumn;
  }
 
  double toDouble(size_t i) const override {
    UNUSED_ARG(i);
    throw std::runtime_error("VectorColumn is not convertible to double.");
  }
 
  void fromDouble(size_t i, double value) override {
    UNUSED_ARG(i);
    UNUSED_ARG(value);
    throw std::runtime_error("VectorColumn is not assignable from double.");
  }
 
  const std::vector<std::vector<Type>> &data() const { return m_data; }
 
  bool equals(const Column &otherColumn, double tolerance) const override {
    if (!possibleToCompare(otherColumn)) {
      return false;
    }
    const auto &otherColumnTyped = static_cast<const VectorColumn<Type> &>(otherColumn);
    const auto &otherData = otherColumnTyped.data();
    for (size_t i = 0; i < m_data.size(); i++) {
      if (m_data[i].size() != otherData[i].size()) {
        return false;
      }
      for (size_t j = 0; j < m_data[i].size(); j++) {
        if (fabs((double)m_data[i][j] - (double)otherData[i][j]) > tolerance) {
          return false;
        }
      }
    }
    return true;
  }
 
  bool equalsRelErr(const Column &otherColumn, double tolerance) const override {
    if (!possibleToCompare(otherColumn)) {
      return false;
    }
    const auto &otherColumnTyped = static_cast<const VectorColumn<Type> &>(otherColumn);
    const auto &otherData = otherColumnTyped.data();
    for (size_t i = 0; i < m_data.size(); i++) {
      if (m_data[i].size() != otherData[i].size()) {
        return false;
      }
      for (size_t j = 0; j < m_data[i].size(); j++) {
        double num = fabs((double)m_data[i][j] - (double)otherData[i][j]);
        double den = (fabs((double)m_data[i][j]) + fabs((double)otherData[i][j])) / 2;
        if (den < tolerance && num > tolerance) {
          return false;
        } else if (num / den > tolerance) {
          return false;
        }
      }
    }
    return true;
  }
 
protected:
  void resize(size_t count) override { m_data.resize(count); }
 
  void insert(size_t index) override {
    // Insert empty vector at the given position
    m_data.insert(m_data.begin() + index, std::vector<Type>());
  }
 
  void remove(size_t index) override { m_data.erase(m_data.begin() + index); }
 
  void *void_pointer(size_t index) override { return &m_data.at(index); }
 
  const void *void_pointer(size_t index) const override { return &m_data.at(index); }
 
private:
  std::vector<std::vector<Type>> m_data;
 
  std::string typeName();
};
 
} // namespace DataObjects
} // namespace Mantid
 
#define DECLARE_VECTORCOLUMN(Type, TypeName)                                                                           \
  template <> std::string VectorColumn<Type>::typeName() { return #TypeName; }                                         \
  Kernel::RegistrationHelper register_column_##TypeName(                                                               \
      (API::ColumnFactory::Instance().subscribe<VectorColumn<Type>>(#TypeName), 0));