PelePhysics.H

#ifndef PELEPHYSICS_H
#define PELEPHYSICS_H

#include <type_traits>
#include <iostream>
#include "PhysicsConstants.H"
#include "PelePhysicsConstraints.H"
#include "EOS.H"
#include "Transport.H"

namespace pele::physics {

template <typename T>
struct is_eos_type
{
  static constexpr bool value = std::is_base_of_v<typename T::eos_type, T>;
};

template <typename T>
struct is_transport_type
{
  static constexpr bool value =
    std::is_base_of_v<typename T::transport_type, T>;
};

template <typename EosModel, typename TransportModel>
struct PelePhysics
{
  static_assert(
    is_eos_type<EosModel>::value,
    "PelePhysics must have EOS model as its first argument");
  static_assert(
    is_transport_type<TransportModel>::value,
    "PelePhysics must have Transport model as its second argument");
  static_assert(
    is_valid_physics_combination<EosModel, TransportModel>::value,
    "Invalid physics combination attempted");

  using eos_type = typename EosModel::eos_type;
  using transport_type = typename TransportModel::transport_type;

  static std::string identifier()
  {
    if (std::is_same_v<EosModel, TransportModel>) {
      return EosModel::identifier();
    }
    return EosModel::identifier() + "-" + TransportModel::identifier();
  }

  template <class... Args>
  AMREX_GPU_HOST_DEVICE static eos_type eos(Args... args)
  {
    return eos_type(std::forward<Args>(args)...);
  }

  template <class... Args>
  AMREX_GPU_HOST_DEVICE static transport_type transport(Args... args)
  {
    return transport_type(std::forward<Args>(args)...);
  }
};

using PhysicsType = PelePhysics<EosType, TransportType>;

} // namespace pele::physics

#endif