EigenFortranVector.h

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2016 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 <utility>
 
namespace Mantid {
namespace CurveFitting {
 
template <class VectorClass> class FortranVector : public VectorClass {
  int m_base;
  using ElementConstType = decltype(std::declval<const VectorClass>()[0]);
  using ElementRefType = decltype(std::declval<VectorClass>()[0]);
 
public:
  FortranVector();
  explicit FortranVector(const int n);
  FortranVector(const int iFrom, const int iTo);
  void allocate(int firstIndex, int lastIndex);
  void allocate(int newSize);
  int len() const;
 
  ElementConstType operator()(int i) const;
  ElementRefType operator()(int i);
  ElementConstType operator[](int i) const;
  ElementRefType operator[](int i);
 
  VectorClass moveToBaseVector();
 
private:
  static size_t makeSize(int firstIndex, int lastIndex);
};
 
template <class VectorClass> size_t FortranVector<VectorClass>::makeSize(int firstIndex, int lastIndex) {
  if (lastIndex < firstIndex) {
    throw std::invalid_argument("Vector defined with invalid index range.");
  }
  return static_cast<size_t>(lastIndex - firstIndex + 1);
}
 
template <class VectorClass> FortranVector<VectorClass>::FortranVector() : VectorClass(makeSize(1, 1)), m_base(1) {}
 
template <class VectorClass>
FortranVector<VectorClass>::FortranVector(const int n) : VectorClass(makeSize(1, n)), m_base(1) {}
 
template <class VectorClass>
FortranVector<VectorClass>::FortranVector(const int iFirst, const int iLast)
    : VectorClass(makeSize(iFirst, iLast)), m_base(iFirst) {}
 
template <class VectorClass> void FortranVector<VectorClass>::allocate(int iFirst, int iLast) {
  m_base = iFirst;
  this->resize(makeSize(iFirst, iLast));
}
 
template <class VectorClass> void FortranVector<VectorClass>::allocate(int newSize) {
  m_base = 1;
  this->resize(makeSize(1, newSize));
}
 
template <class VectorClass>
typename FortranVector<VectorClass>::ElementConstType FortranVector<VectorClass>::operator()(int i) const {
  return this->VectorClass::operator[](static_cast<size_t>(i - m_base));
}
 
template <class VectorClass>
typename FortranVector<VectorClass>::ElementRefType FortranVector<VectorClass>::operator()(int i) {
  return this->VectorClass::operator[](static_cast<size_t>(i - m_base));
}
 
template <class VectorClass>
typename FortranVector<VectorClass>::ElementConstType FortranVector<VectorClass>::operator[](int i) const {
  return this->VectorClass::operator[](static_cast<size_t>(i - m_base));
}
 
template <class VectorClass>
typename FortranVector<VectorClass>::ElementRefType FortranVector<VectorClass>::operator[](int i) {
  return this->VectorClass::operator[](static_cast<size_t>(i - m_base));
}
 
template <class VectorClass> VectorClass FortranVector<VectorClass>::moveToBaseVector() {
  return VectorClass(std::move(*this));
}
 
template <class VectorClass> int FortranVector<VectorClass>::len() const { return static_cast<int>(this->size()); }
 
} // namespace CurveFitting
} // namespace Mantid