system_base.hpp

Go to the documentation of this file.

// Copyright (c) Lawrence Livermore National Security, LLC and
// other Smith Project Developers. See the top-level LICENSE file for
// details.
//
// SPDX-License-Identifier: (BSD-3-Clause)
 
#pragma once
 
#include <string>
#include <vector>
#include <memory>
#include <utility>
#include "field_state.hpp"
#include "field_store.hpp"
#include "coupled_system_solver.hpp"
#include "state_advancer.hpp"
#include "smith/physics/common.hpp"
#include "mfem.hpp"
 
namespace smith {
 
template <typename... Args>
struct Parameters;
 
struct ReactionInfo {
  std::string name;                                    
  const mfem::ParFiniteElementSpace* space = nullptr;  
};
 
namespace detail {
 
template <std::size_t, typename T>
using always_t = T;
 
template <typename Space, typename... Tail, std::size_t... Is>
auto time_rule_params_impl(std::index_sequence<Is...>) -> Parameters<always_t<Is, Space>..., Tail...>;
 
}  // namespace detail
 
template <typename Rule, typename Space, typename... Tail>
using TimeRuleParams =
    decltype(detail::time_rule_params_impl<Space, Tail...>(std::make_index_sequence<Rule::num_states>{}));
 
struct SystemBase {
  std::shared_ptr<FieldStore> field_store;      
  std::shared_ptr<CoupledSystemSolver> solver;  
  std::shared_ptr<StateAdvancer> advancer;      
  std::vector<FieldState> parameter_fields;     
  std::string prepend_name;                     
 
  const std::vector<FieldState>& getParameterFields() const { return parameter_fields; }
 
  std::string prefix(const std::string& name) const
  {
    if (prepend_name.empty()) {
      return name;
    }
    return prepend_name + "_" + name;
  }
 
  std::vector<ReactionInfo> getReactionInfos() const { return {}; }
};
 
}  // namespace smith