MultiBlockGriddingAlgorithm.hpp

Go to the documentation of this file.

// Copyright (c) 2019 Maikel Nadolski
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of this software and associated documentation files (the "Software"), to deal
// in the Software without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of the Software, and to permit persons to whom the Software is
// furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included in
// all copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
// SOFTWARE.
 
#ifndef FUB_AMREX_COUPLED_GRIDDING_ALGORITHM_HPP
#define FUB_AMREX_COUPLED_GRIDDING_ALGORITHM_HPP
 
#include "fub/AMReX/multi_block/MultiBlockBoundary.hpp"
 
#include "fub/AMReX/GriddingAlgorithm.hpp"
#include "fub/AMReX/cutcell/GriddingAlgorithm.hpp"
 
#include "fub/output/CounterOutput.hpp"
 
#include "fub/core/span.hpp"
 
#include "fub/equations/ideal_gas_mix/FlameMasterReactor.hpp"
 
namespace fub {
namespace amrex {
 
/// \ingroup GriddingAlgorithm
class MultiBlockGriddingAlgorithm {
public:
  MultiBlockGriddingAlgorithm(
      FlameMasterReactor reactor,
      std::vector<std::shared_ptr<GriddingAlgorithm>> tubes,
      std::vector<std::shared_ptr<cutcell::GriddingAlgorithm>> plena,
      std::vector<BlockConnection> connectivity);
 
  MultiBlockGriddingAlgorithm(
      FlameMasterReactor reactor,
      std::vector<std::shared_ptr<GriddingAlgorithm>> tubes,
      std::vector<std::shared_ptr<cutcell::GriddingAlgorithm>> plena,
      std::vector<BlockConnection> connectivity,
      std::vector<std::shared_ptr<PressureValve>> valves);
 
  MultiBlockGriddingAlgorithm(const MultiBlockGriddingAlgorithm& other);
  MultiBlockGriddingAlgorithm&
  operator=(const MultiBlockGriddingAlgorithm& other);
 
  MultiBlockGriddingAlgorithm(MultiBlockGriddingAlgorithm&& other) noexcept =
      default;
 
  MultiBlockGriddingAlgorithm&
  operator=(MultiBlockGriddingAlgorithm&& other) noexcept = default;
 
  std::ptrdiff_t GetCycles() const noexcept { return plena_[0]->GetCycles(); }
 
  Duration GetTimePoint() const noexcept { return plena_[0]->GetTimePoint(); }
 
  [[nodiscard]] span<const std::shared_ptr<GriddingAlgorithm>> GetTubes() const
      noexcept;
  [[nodiscard]] span<const std::shared_ptr<cutcell::GriddingAlgorithm>>
  GetPlena() const noexcept;
 
  [[nodiscard]] span<const BlockConnection> GetConnectivity() const noexcept;
  [[nodiscard]] span<MultiBlockBoundary> GetBoundaries(int level = 0) noexcept {
    return boundaries_[static_cast<std::size_t>(level)];
  }
 
  void RegridAllFinerLevels(int which_level);
 
private:
  FlameMasterReactor reactor_;
  std::vector<std::shared_ptr<GriddingAlgorithm>> tubes_;
  std::vector<std::shared_ptr<cutcell::GriddingAlgorithm>> plena_;
  std::vector<BlockConnection> connectivity_;
  std::vector<std::vector<MultiBlockBoundary>> boundaries_;
};
 
} // namespace amrex
 
template <typename PrintDuration>
class CounterOutput<amrex::MultiBlockGriddingAlgorithm, PrintDuration>
    : public OutputAtFrequencyOrInterval<amrex::MultiBlockGriddingAlgorithm> {
public:
  CounterOutput(const ProgramOptions& po,
                std::chrono::steady_clock::time_point ref)
      : OutputAtFrequencyOrInterval<amrex::MultiBlockGriddingAlgorithm>(po),
        reference_{ref} {}
 
  CounterOutput(const ProgramOptions& po)
      : OutputAtFrequencyOrInterval<amrex::MultiBlockGriddingAlgorithm>(po),
        reference_{std::chrono::steady_clock::now()} {}
 
  CounterOutput(std::chrono::steady_clock::time_point reference,
                std::vector<std::ptrdiff_t> frequencies,
                std::vector<Duration> intervals)
      : OutputAtFrequencyOrInterval<amrex::MultiBlockGriddingAlgorithm>(
            std::move(frequencies), std::move(intervals)),
        reference_(reference) {}
 
  void operator()(const amrex::MultiBlockGriddingAlgorithm& grid) override {
    std::chrono::steady_clock::time_point now =
        std::chrono::steady_clock::now();
    auto diff =
        std::chrono::duration_cast<std::chrono::nanoseconds>(now - reference_);
    std::vector<CounterResult> statistics = grid.GetPlena()[0]
                                                ->GetPatchHierarchy()
                                                .GetCounterRegistry()
                                                ->gather_statistics();
    if (statistics.size()) {
      print_statistics<PrintDuration>(statistics, diff.count());
    }
  }
 
private:
  std::chrono::steady_clock::time_point reference_;
};
 
} // namespace fub
 
#endif // FUB_AMREX_COUPLED_GRIDDING_ALGORITHM_HPP