CCZ4SBH_FV.h

Go to the documentation of this file.

// This file is part of the ExaHyPE2 project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include "AbstractCCZ4SBH_FV.h"
 
#include "MulticoreOrchestration.h"
 
#include "tarch/logging/Log.h"
 
namespace benchmarks {
  namespace exahype2 {
    namespace ccz4 {
      class CCZ4SBH_FV;
#ifdef PureFV
      extern TP::TwoPunctures* twoPunctures;
 
      void prepareTwoPunctures();
#endif
    } // namespace ccz4
  }   // namespace exahype2
} // namespace benchmarks
 
 
 
class benchmarks::exahype2::ccz4::CCZ4SBH_FV: public benchmarks::exahype2::ccz4::AbstractCCZ4SBH_FV {
private:
  static tarch::logging::Log                _log;
  static tarch::multicore::BooleanSemaphore _semaphore;
 
  int _numberOfPatches;
 
public:
  CCZ4SBH_FV();
 
  virtual void initialCondition(
    double* NOALIAS Q,
    const tarch::la::Vector<DIMENSIONS, double>& volumeCentre,
    const tarch::la::Vector<DIMENSIONS, double>& volumeH,
    bool                                         gridIsConstructed
  ) override;

  void reduceAdmissibleTimeStepSize(double timeStepSize);

  static bool isCellOverlappingWithBHImpactArea(const peano4::datamanagement::CellMarker& marker);
 
  static bool areAllFaceConnectedCellsOverlappingWithBHImpactArea(const peano4::datamanagement::CellMarker& marker);

  static bool areBothAdjacentCellsOverlappingWithBHImpactArea(const peano4::datamanagement::FaceMarker& marker);
 
  static bool isOneAdjacentCellOverlappingWithBHImpactArea(const peano4::datamanagement::FaceMarker& marker);

  static bool areBothAdjacentCellsOverlappingWithBHImpactArea(
    const peano4::datamanagement::CellMarker& marker, int faceNumber
  );

  virtual void startTimeStep(
    double globalMinTimeStamp, double globalMaxTimeStamp, double globalMinTimeStepSize, double globalMaxTimeStepSize
  ) override;

  virtual void finishTimeStep() override;
 
  void incNumberOfPatches();
 
private:
  static bool isCellOverlappingWithBHImpactArea(
    const tarch::la::Vector<DIMENSIONS, double>& cellCentre, const tarch::la::Vector<DIMENSIONS, double>& cellH
  );
};