CellMarker.h

Go to the documentation of this file.

// This file is part of the Peano project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include "peano4/grid/GridTraversalEvent.h"
#include "peano4/utils/Globals.h"
#include "tarch/la/Vector.h"
 
#ifdef Parallel
#include <functional>
#include <mpi.h>
#endif
 
namespace peano4 {
  namespace datamanagement {
    struct CellMarker;
  } // namespace datamanagement
} // namespace peano4
 
std::ostream& operator<<(std::ostream& out, const peano4::datamanagement::CellMarker& marker);
 
struct peano4::datamanagement::CellMarker {
private:
  tarch::la::Vector<Dimensions, double> _centre;
 
  tarch::la::Vector<Dimensions, double> _h;
 
  bool _hasBeenRefined;
  bool _willBeRefined;
  bool _isLocal;
  bool _isParentLocal;
  bool _areAllVerticesRefined;
  bool _isOneVertexHanging;
  bool _isOneVertexCreatedOrDestroyed;
 
  bool _areAllVerticesInsideDomain;
  bool _invokingSpacetreeIsNotInvolvedInAnyDynamicLoadBalancing;
 
  bool _willBeEnclave;
  bool _hasBeenEnclave;
 
  tarch::la::Vector<Dimensions, int> _relativePositionOfCellWithinFatherCell;
 
public:
  CellMarker(const peano4::grid::GridTraversalEvent& event);
 
  bool hasBeenRefined() const;
 
  bool willBeRefined() const;
 
  tarch::la::Vector<Dimensions, double> x() const;
 
  tarch::la::Vector<Dimensions, double> getInvokingParentCellsCentre() const;
 
  bool isContained(const tarch::la::Vector<Dimensions, double>& x, double tolerance = tarch::la::NUMERICAL_ZERO_DIFFERENCE) const;
 
  static bool isContained(
    const tarch::la::Vector<Dimensions, double>&   x,
    const tarch::la::Vector<Dimensions, double>&   cellCentre,
    const tarch::la::Vector<Dimensions, double>&   h,
    double                                         tolerance
  );
 
  bool overlaps(const tarch::la::Vector<Dimensions, double>& offset, const tarch::la::Vector<Dimensions, double>& size) const;
 
  tarch::la::Vector<Dimensions, double> h() const;
 
  tarch::la::Vector<Dimensions, double> getOffset() const;
 
  std::string toString() const;
 
  bool isLocal() const;
 
  bool isParentLocal() const;
 
  bool hasBeenEnclaveCell() const;
 
  bool willBeEnclaveCell() const;
 
  bool hasBeenSkeletonCell() const;
 
  bool willBeSkeletonCell() const;
 
  tarch::la::Vector<Dimensions, int> getRelativePositionWithinFatherCell() const;
 
#if PeanoDebug > 0
  void setRelativePositionWithinFatherCell(int axis, int value);
#endif
 
#ifdef Parallel
  static void initDatatype();
  static void shutdownDatatype();
 
  static void send(const CellMarker& buffer, int destination, int tag, MPI_Comm communicator);
  static void receive(CellMarker& buffer, int source, int tag, MPI_Comm communicator);
 
  static void send(const CellMarker& buffer, int destination, int tag, std::function<void()> startCommunicationFunctor, std::function<void()> waitFunctor, MPI_Comm communicator);
  static void receive(CellMarker& buffer, int source, int tag, std::function<void()> startCommunicationFunctor, std::function<void()> waitFunctor, MPI_Comm communicator);
 
  static void sendAndPollDanglingMessages(const CellMarker& message, int destination, int tag, MPI_Comm communicator);
  static void receiveAndPollDanglingMessages(CellMarker& message, int source, int tag, MPI_Comm communicator);
 
  static MPI_Datatype Datatype;
#endif
};