Polation.hpp

Go to the documentation of this file.

#pragma once
 
#include <boost/container/static_vector.hpp>
#include <iosfwd>
#include "Eigen/Core"
#include "mesh/Edge.hpp"
#include "mesh/Tetrahedron.hpp"
#include "mesh/Triangle.hpp"
#include "mesh/Vertex.hpp"
 
namespace precice::mapping {

struct WeightedElement {
  int    vertexID;
  double weight;
};

class Polation {
public:
  Polation(const Eigen::VectorXd &location, const mesh::Vertex &element);

  Polation(const Eigen::VectorXd &location, const mesh::Edge &element);

  Polation(const Eigen::VectorXd &location, const mesh::Triangle &element);

  Polation(const Eigen::VectorXd &location, const mesh::Tetrahedron &element);

  std::size_t nElements() const;

  const boost::container::static_vector<WeightedElement, 4> &getWeightedElements() const;

  bool isInterpolation() const;

  double distance() const;
 
private:
  boost::container::static_vector<WeightedElement, 4> _weightedElements;
  double                                              _distance;
};

std::ostream &operator<<(std::ostream &os, const WeightedElement &w);

std::ostream &operator<<(std::ostream &os, const Polation &p);
 
} // namespace precice::mapping