CD_KMCSingleStateReaction.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_KMCSINGLESTATEREACTION_H
#define CD_KMCSINGLESTATEREACTION_H
 
// Std includes
#include <map>
#include <vector>
#include <list>
 
// Chombo includes
#include <REAL.H>
 
// Our includes
#include <CD_KMCSingleState.H>
#include <CD_NamespaceHeader.H>
 
template <typename State = KMCSingleState<long long>, typename T = long long>
class KMCSingleStateReaction
{
public:
  KMCSingleStateReaction() = default;
 
  KMCSingleStateReaction(const KMCSingleStateReaction&) = default;
 
  inline KMCSingleStateReaction(const std::list<size_t>& a_reactants, const std::list<size_t>& a_products) noexcept;
 
  inline virtual ~KMCSingleStateReaction();
 
  inline Real&
  rate() const noexcept;
 
  static inline T
  population(const size_t& a_reactant, const State& a_state) noexcept;
 
  inline Real
  propensity(const State& a_state) const noexcept;
 
  inline T
  computeCriticalNumberOfReactions(const State& a_state) const noexcept;
 
  inline std::list<size_t>
  getReactants() const noexcept;
 
  inline T
  getStateChange(const size_t a_particleReactant) const noexcept;
 
  inline void
  advanceState(State& a_state, const T& a_numReactions) const noexcept;
 
protected:
  mutable Real m_rate;
 
  Real m_propensityFactor;
 
  std::list<size_t> m_reactants;
 
  std::list<size_t> m_products;
 
  std::map<size_t, T> m_stateChange;
 
  inline void
  computeStateChanges() noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_KMCSingleStateReactionImplem.H>
 
#endif