Tetrahedron.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <iostream>
#include <tuple>
 
#include "mesh/Vertex.hpp"
#include "precice/impl/Types.hpp"
#include "utils/assertion.hpp"
 
// ----------------------------------------------------------- CLASS DEFINITION
 
namespace precice {
namespace mesh {
 
class Tetrahedron {
public:
  static constexpr int vertexCount{4};
 
  Tetrahedron(
      Vertex &vertexOne,
      Vertex &vertexTwo,
      Vertex &vertexThree,
      Vertex &vertexFour);
 
  int getDimensions() const;
 
  Vertex &vertex(int i);
 
  const Vertex &vertex(int i) const;
 
  double getVolume() const;
 
  const Eigen::VectorXd getCenter() const;
 
  double getEnclosingRadius() const;
 
  bool operator==(const Tetrahedron &other) const;
 
  bool operator!=(const Tetrahedron &other) const;
 
  bool operator<(const Tetrahedron &other) const
  {
    return std::make_tuple(_vertices[0]->getID(), _vertices[1]->getID(), _vertices[2]->getID(), _vertices[3]->getID()) <
           std::make_tuple(other._vertices[0]->getID(), other._vertices[1]->getID(), other._vertices[2]->getID(), other._vertices[3]->getID());
  }
 
private:
  std::array<Vertex *, 4> _vertices;
};
 
// --------------------------------------------------------- HEADER DEFINITIONS
 
inline Vertex &Tetrahedron::vertex(int i)
{
  PRECICE_ASSERT((i >= 0) && (i < 4), i);
  return *_vertices[i];
}
 
inline const Vertex &Tetrahedron::vertex(int i) const
{
  PRECICE_ASSERT((i >= 0) && (i < 4), i);
  return *_vertices[i];
}
 
std::ostream &operator<<(std::ostream &os, const Tetrahedron &t);
 
} // namespace mesh
} // namespace precice