Sample.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <Eigen/Core>
 
#include "utils/assertion.hpp"
 
namespace precice::time {

struct Sample {
  explicit Sample(int dims) noexcept
      : dataDims(dims) {}

  Sample(int dims, int dataCount)
      : dataDims(dims), values(dims * dataCount) {}

  Sample(int dataDims, int nVertices, int meshDims)
      : dataDims(dataDims), values(nVertices * dataDims), gradients(meshDims, nVertices * dataDims)
  {
  }

  Sample(int dims, Eigen::VectorXd inValues)
      : dataDims(dims), values(std::move(inValues))
  {
    PRECICE_ASSERT(dataDims > 0);
  }

  Sample(int dims, Eigen::VectorXd inValues, Eigen::MatrixXd inGradients)
      : dataDims(dims), values(std::move(inValues)), gradients(std::move(inGradients))
  {
    PRECICE_ASSERT(dataDims > 0);
    PRECICE_ASSERT(gradients.size() == 0 || gradients.cols() == values.size(), gradients.size(), gradients.cols(), values.size());
  }
 
  Sample(const Sample &) = default;
  Sample(Sample &&)      = default;
 
  Sample &operator=(const Sample &) = default;
  Sample &operator=(Sample &&)      = default;

  Sample &setZero()
  {
    values.setZero();
    gradients.setZero();
    return *this;
  }

  int dataDims;

  Eigen::VectorXd values;

  Eigen::MatrixXd gradients;
};
 
} // namespace precice::time