A reaction network engine that uses a graph-based representation. More...

#include <engine_graph.h>

Inheritance diagram for gridfire::GraphEngine:

Classes
struct	constants

struct	PrecomputedReaction

Public Member Functions
	GraphEngine (const fourdst::composition::Composition &composition)
	Constructs a GraphEngine from a composition.

	GraphEngine (const reaction::LogicalReactionSet &reactions)
	Constructs a GraphEngine from a set of reactions.

StepDerivatives< double >	calculateRHSAndEnergy (const std::vector< double > &Y, const double T9, const double rho) const override
	Calculates the right-hand side (dY/dt) and energy generation rate.

void	generateJacobianMatrix (const std::vector< double > &Y, const double T9, const double rho) override
	Generates the Jacobian matrix for the current state.

void	generateStoichiometryMatrix () override
	Generates the stoichiometry matrix for the network.

double	calculateMolarReactionFlow (const reaction::Reaction &reaction, const std::vector< double > &Y, const double T9, const double rho) const override
	Calculates the molar reaction flow for a given reaction.

const std::vector< fourdst::atomic::Species > &	getNetworkSpecies () const override
	Gets the list of species in the network.

const reaction::LogicalReactionSet &	getNetworkReactions () const override
	Gets the set of logical reactions in the network.

double	getJacobianMatrixEntry (const int i, const int j) const override
	Gets an entry from the previously generated Jacobian matrix.

int	getStoichiometryMatrixEntry (const int speciesIndex, const int reactionIndex) const override
	Gets an entry from the stoichiometry matrix.

std::unordered_map< fourdst::atomic::Species, double >	getSpeciesTimescales (const std::vector< double > &Y, double T9, double rho) const override
	Computes timescales for all species in the network.

void	update (const NetIn &netIn) override
	Update the internal state of the engine.

bool	involvesSpecies (const fourdst::atomic::Species &species) const
	Checks if a given species is involved in the network.

void	exportToDot (const std::string &filename) const
	Exports the network to a DOT file for visualization.

void	exportToCSV (const std::string &filename) const
	Exports the network to a CSV file for analysis.

void	setScreeningModel (screening::ScreeningType) override
	Set the electron screening model.

screening::ScreeningType	getScreeningModel () const override
	Get the current electron screening model.

void	setPrecomputation (bool precompute)

bool	isPrecomputationEnabled () const

Public Member Functions inherited from gridfire::Engine
virtual	~Engine ()=default
	Virtual destructor.

Static Public Member Functions
static std::unordered_map< fourdst::atomic::Species, int >	getNetReactionStoichiometry (const reaction::Reaction &reaction)
	Gets the net stoichiometry for a given reaction.

Private Member Functions
void	syncInternalMaps ()
	Synchronizes the internal maps.

void	collectNetworkSpecies ()
	Collects the unique species in the network.

void	populateReactionIDMap ()
	Populates the reaction ID map.

void	populateSpeciesToIndexMap ()
	Populates the species-to-index map.

void	reserveJacobianMatrix ()
	Reserves space for the Jacobian matrix.

void	recordADTape ()
	Records the AD tape for the right-hand side of the ODE.

void	precomputeNetwork ()

bool	validateConservation () const
	Validates mass and charge conservation across all reactions.

void	validateComposition (const fourdst::composition::Composition &composition, double culling, double T9)
	Validates the composition against the current reaction set.

StepDerivatives< double >	calculateAllDerivativesUsingPrecomputation (const std::vector< double > &Y_in, const std::vector< double > &bare_rates, double T9, double rho) const

template<IsArithmeticOrAD T>
T	calculateMolarReactionFlow (const reaction::Reaction &reaction, const std::vector< T > &Y, const T T9, const T rho) const
	Calculates the molar reaction flow for a given reaction.

template<IsArithmeticOrAD T>
StepDerivatives< T >	calculateAllDerivatives (const std::vector< T > &Y_in, T T9, T rho) const
	Calculates all derivatives (dY/dt) and the energy generation rate.

StepDerivatives< double >	calculateAllDerivatives (const std::vector< double > &Y_in, const double T9, const double rho) const
	Calculates all derivatives (dY/dt) and the energy generation rate (double precision).

StepDerivatives< ADDouble >	calculateAllDerivatives (const std::vector< ADDouble > &Y_in, const ADDouble &T9, const ADDouble &rho) const
	Calculates all derivatives (dY/dt) and the energy generation rate (automatic differentiation).

Private Attributes
Config &	m_config = Config::getInstance()

quill::Logger *	m_logger = LogManager::getInstance().getLogger("log")

constants	m_constants

reaction::LogicalReactionSet	m_reactions
	Set of REACLIB reactions in the network.

std::unordered_map< std::string_view, reaction::Reaction * >	m_reactionIDMap
	Map from reaction ID to REACLIBReaction. //PERF: This makes copies of REACLIBReaction and could be a performance bottleneck.

std::vector< fourdst::atomic::Species >	m_networkSpecies
	Vector of unique species in the network.

std::unordered_map< std::string_view, fourdst::atomic::Species >	m_networkSpeciesMap
	Map from species name to Species object.

std::unordered_map< fourdst::atomic::Species, size_t >	m_speciesToIndexMap
	Map from species to their index in the stoichiometry matrix.

boost::numeric::ublas::compressed_matrix< int >	m_stoichiometryMatrix
	Stoichiometry matrix (species x reactions).

boost::numeric::ublas::compressed_matrix< double >	m_jacobianMatrix
	Jacobian matrix (species x species).

CppAD::ADFun< double >	m_rhsADFun
	CppAD function for the right-hand side of the ODE.

screening::ScreeningType	m_screeningType = screening::ScreeningType::BARE
	Screening type for the reaction network. Default to no screening.

std::unique_ptr< screening::ScreeningModel >	m_screeningModel = screening::selectScreeningModel(m_screeningType)

bool	m_usePrecomputation = true
	Flag to enable or disable using precomputed reactions for efficiency. Mathematically, this should not change the results. Generally end users should not need to change this.

std::vector< PrecomputedReaction >	m_precomputedReactions
	Precomputed reactions for efficiency.

Detailed Description

A reaction network engine that uses a graph-based representation.

The GraphEngine class implements the DynamicEngine interface using a graph-based representation of the reaction network. It uses sparse matrices for efficient storage and computation of the stoichiometry and Jacobian matrices. Automatic differentiation (AD) is used to calculate the Jacobian matrix.

The engine supports:

Calculation of the right-hand side (dY/dt) and energy generation rate.
Generation and access to the Jacobian matrix.
Generation and access to the stoichiometry matrix.
Calculation of molar reaction flows.
Access to the set of logical reactions in the network.
Computation of timescales for each species.
Exporting the network to DOT and CSV formats for visualization and analysis.

See also: engine_abstract.h

Definition at line 90 of file engine_graph.h.

Constructor & Destructor Documentation

◆ GraphEngine() [1/2]

gridfire::GraphEngine::GraphEngine ( const fourdst::composition::Composition & composition )

explicit

Constructs a GraphEngine from a composition.

Parameters

composition The composition of the material.

This constructor builds the reaction network from the given composition using the build_reaclib_nuclear_network function.

See also: build_reaclib_nuclear_network

Definition at line 26 of file engine_graph.cpp.

◆ GraphEngine() [2/2]

gridfire::GraphEngine::GraphEngine ( const reaction::LogicalReactionSet & reactions )

explicit

Constructs a GraphEngine from a set of reactions.

Parameters

reactions The set of reactions to use in the network.

This constructor uses the given set of reactions to construct the reaction network.

Definition at line 34 of file engine_graph.cpp.

Member Function Documentation

◆ calculateAllDerivatives() [1/3]

StepDerivatives< ADDouble > gridfire::GraphEngine::calculateAllDerivatives	(	const std::vector< ADDouble > &	Y_in,
		const ADDouble &	T9,
		const ADDouble &	rho ) const

nodiscardprivate

Calculates all derivatives (dY/dt) and the energy generation rate (automatic differentiation).

Parameters

Y_in	Vector of current abundances for all species.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: StepDerivatives<ADDouble> containing dY/dt and energy generation rate.

This method calculates the time derivatives of all species and the specific nuclear energy generation rate for the current state using automatic differentiation.

Definition at line 352 of file engine_graph.cpp.

◆ calculateAllDerivatives() [2/3]

StepDerivatives< double > gridfire::GraphEngine::calculateAllDerivatives	(	const std::vector< double > &	Y_in,
		const double	T9,
		const double	rho ) const

nodiscardprivate

Calculates all derivatives (dY/dt) and the energy generation rate (double precision).

Parameters

Y_in	Vector of current abundances for all species.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: StepDerivatives<double> containing dY/dt and energy generation rate.

This method calculates the time derivatives of all species and the specific nuclear energy generation rate for the current state using double precision arithmetic.

Definition at line 344 of file engine_graph.cpp.

◆ calculateAllDerivatives() [3/3]

template<IsArithmeticOrAD T>

StepDerivatives< T > gridfire::GraphEngine::calculateAllDerivatives	(	const std::vector< T > &	Y_in,
		T	T9,
		T	rho ) const

nodiscardprivate

Calculates all derivatives (dY/dt) and the energy generation rate.

Template Parameters

T	The numeric type to use for the calculation.

Parameters

Y_in	Vector of current abundances for all species.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: StepDerivatives<T> containing dY/dt and energy generation rate.

This method calculates the time derivatives of all species and the specific nuclear energy generation rate for the current state.

Definition at line 521 of file engine_graph.h.

◆ calculateAllDerivativesUsingPrecomputation()

StepDerivatives< double > gridfire::GraphEngine::calculateAllDerivativesUsingPrecomputation	(	const std::vector< double > &	Y_in,
		const std::vector< double > &	bare_rates,
		double	T9,
		double	rho ) const

nodiscardprivate

Definition at line 223 of file engine_graph.cpp.

◆ calculateMolarReactionFlow() [1/2]

double gridfire::GraphEngine::calculateMolarReactionFlow	(	const reaction::Reaction &	reaction,
		const std::vector< double > &	Y,
		const double	T9,
		const double	rho ) const

nodiscardoverridevirtual

Calculates the molar reaction flow for a given reaction.

Parameters

reaction	The reaction for which to calculate the flow.
Y	Vector of current abundances.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: Molar flow rate for the reaction (e.g., mol/g/s).

This method computes the net rate at which the given reaction proceeds under the current state.

Implements gridfire::DynamicEngine.

Definition at line 377 of file engine_graph.cpp.

◆ calculateMolarReactionFlow() [2/2]

template<IsArithmeticOrAD T>

T gridfire::GraphEngine::calculateMolarReactionFlow	(	const reaction::Reaction &	reaction,
		const std::vector< T > &	Y,
		const T	T9,
		const T	rho ) const

private

Calculates the molar reaction flow for a given reaction.

Template Parameters

T	The numeric type to use for the calculation.

Parameters

reaction	The reaction for which to calculate the flow.
Y	Vector of current abundances.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: Molar flow rate for the reaction (e.g., mol/g/s).

This method computes the net rate at which the given reaction proceeds under the current state.

Definition at line 589 of file engine_graph.h.

◆ calculateRHSAndEnergy()

StepDerivatives< double > gridfire::GraphEngine::calculateRHSAndEnergy	(	const std::vector< double > &	Y,
		const double	T9,
		const double	rho ) const

nodiscardoverridevirtual

Calculates the right-hand side (dY/dt) and energy generation rate.

Parameters

Y	Vector of current abundances for all species.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: StepDerivatives<double> containing dY/dt and energy generation rate.

This method calculates the time derivatives of all species and the specific nuclear energy generation rate for the current state.

See also: StepDerivatives

Implements gridfire::Engine.

Definition at line 42 of file engine_graph.cpp.

◆ collectNetworkSpecies()

void gridfire::GraphEngine::collectNetworkSpecies ( )

private

Collects the unique species in the network.

This method collects the unique species in the network from the reactants and products of all reactions.

Definition at line 72 of file engine_graph.cpp.

◆ exportToCSV()

void gridfire::GraphEngine::exportToCSV ( const std::string & filename ) const

Exports the network to a CSV file for analysis.

Parameters

filename The name of the CSV file to create.

This method generates a CSV file containing information about the reactions in the network, including the reactants, products, Q-value, and reaction rate coefficients.

Exceptions

std::runtime_error If the file cannot be opened for writing.

Example usage:

engine.exportToCSV("network.csv");

Definition at line 484 of file engine_graph.cpp.

◆ exportToDot()

void gridfire::GraphEngine::exportToDot ( const std::string & filename ) const

Exports the network to a DOT file for visualization.

Parameters

filename The name of the DOT file to create.

This method generates a DOT file that can be used to visualize the reaction network as a graph. The DOT file can be converted to a graphical image using Graphviz.

Exceptions

std::runtime_error If the file cannot be opened for writing.

Example usage:

engine.exportToDot("network.dot");

Definition at line 436 of file engine_graph.cpp.

◆ generateJacobianMatrix()

void gridfire::GraphEngine::generateJacobianMatrix	(	const std::vector< double > &	Y,
		const double	T9,
		const double	rho )

overridevirtual

Generates the Jacobian matrix for the current state.

Parameters

Y	Vector of current abundances.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

This method computes and stores the Jacobian matrix (∂(dY/dt)_i/∂Y_j) for the current state using automatic differentiation. The matrix can then be accessed via getJacobianMatrixEntry().

See also: getJacobianMatrixEntry()

Implements gridfire::DynamicEngine.

Definition at line 386 of file engine_graph.cpp.

◆ generateStoichiometryMatrix()

void gridfire::GraphEngine::generateStoichiometryMatrix ( )

overridevirtual

Generates the stoichiometry matrix for the network.

This method computes and stores the stoichiometry matrix, which encodes the net change of each species in each reaction.

Implements gridfire::DynamicEngine.

Definition at line 303 of file engine_graph.cpp.

◆ getJacobianMatrixEntry()

double gridfire::GraphEngine::getJacobianMatrixEntry	(	const int	i,
		const int	j ) const

nodiscardoverridevirtual

Gets an entry from the previously generated Jacobian matrix.

Parameters

i	Row index (species index).
j	Column index (species index).

Returns: Value of the Jacobian matrix at (i, j).

The Jacobian must have been generated by generateJacobianMatrix() before calling this.

See also: generateJacobianMatrix()

Implements gridfire::DynamicEngine.

Definition at line 419 of file engine_graph.cpp.

◆ getNetReactionStoichiometry()

std::unordered_map< fourdst::atomic::Species, int > gridfire::GraphEngine::getNetReactionStoichiometry ( const reaction::Reaction & reaction )

staticnodiscard

Gets the net stoichiometry for a given reaction.

Parameters

reaction The reaction for which to get the stoichiometry.

Returns: Map of species to their stoichiometric coefficients.

Definition at line 423 of file engine_graph.cpp.

◆ getNetworkReactions()

const reaction::LogicalReactionSet & gridfire::GraphEngine::getNetworkReactions ( ) const

nodiscardoverridevirtual

Gets the set of logical reactions in the network.

Returns: Reference to the LogicalReactionSet containing all reactions.

Implements gridfire::DynamicEngine.

Definition at line 134 of file engine_graph.cpp.

◆ getNetworkSpecies()

const std::vector< fourdst::atomic::Species > & gridfire::GraphEngine::getNetworkSpecies ( ) const

nodiscardoverridevirtual

Gets the list of species in the network.

Returns: Vector of Species objects representing all network species.

Implements gridfire::Engine.

Definition at line 128 of file engine_graph.cpp.

◆ getScreeningModel()

screening::ScreeningType gridfire::GraphEngine::getScreeningModel ( ) const

nodiscardoverridevirtual

Get the current electron screening model.

Returns: The currently active screening model type.

Usage Example:: screening::ScreeningType currentModel = myEngine.getScreeningModel();

gridfire::screening::ScreeningType
ScreeningType
Enumerates the available plasma screening models.
Definition screening_types.h:15

Implements gridfire::DynamicEngine.

Definition at line 365 of file engine_graph.cpp.

◆ getSpeciesTimescales()

std::unordered_map< fourdst::atomic::Species, double > gridfire::GraphEngine::getSpeciesTimescales	(	const std::vector< double > &	Y,
		double	T9,
		double	rho ) const

nodiscardoverridevirtual

Computes timescales for all species in the network.

Parameters

Y	Vector of current abundances.
T9	Temperature in units of 10^9 K.
rho	Density in g/cm^3.

Returns: Map from Species to their characteristic timescales (s).

This method estimates the timescale for abundance change of each species, which can be used for timestep control or diagnostics.

Implements gridfire::DynamicEngine.

Definition at line 538 of file engine_graph.cpp.

◆ getStoichiometryMatrixEntry()

int gridfire::GraphEngine::getStoichiometryMatrixEntry	(	const int	speciesIndex,
		const int	reactionIndex ) const

nodiscardoverridevirtual

Gets an entry from the stoichiometry matrix.

Parameters

speciesIndex	Index of the species.
reactionIndex	Index of the reaction.

Returns: Stoichiometric coefficient for the species in the reaction.

The stoichiometry matrix must have been generated by generateStoichiometryMatrix().

See also: generateStoichiometryMatrix()

Implements gridfire::DynamicEngine.

Definition at line 429 of file engine_graph.cpp.

◆ involvesSpecies()

bool gridfire::GraphEngine::involvesSpecies ( const fourdst::atomic::Species & species ) const

nodiscard

Checks if a given species is involved in the network.

Parameters

species The species to check.

Returns: True if the species is involved in the network, false otherwise.

Definition at line 140 of file engine_graph.cpp.

◆ isPrecomputationEnabled()

bool gridfire::GraphEngine::isPrecomputationEnabled ( ) const

nodiscard

Definition at line 373 of file engine_graph.cpp.

◆ populateReactionIDMap()

void gridfire::GraphEngine::populateReactionIDMap ( )

private

Populates the reaction ID map.

This method populates the reaction ID map, which maps reaction IDs to REACLIBReaction objects.

Definition at line 101 of file engine_graph.cpp.

◆ populateSpeciesToIndexMap()

void gridfire::GraphEngine::populateSpeciesToIndexMap ( )

private

Populates the species-to-index map.

This method populates the species-to-index map, which maps species to their index in the stoichiometry matrix.

Definition at line 110 of file engine_graph.cpp.

◆ precomputeNetwork()

void gridfire::GraphEngine::precomputeNetwork ( )

private

Definition at line 603 of file engine_graph.cpp.

◆ recordADTape()

void gridfire::GraphEngine::recordADTape ( )

private

Records the AD tape for the right-hand side of the ODE.

This method records the AD tape for the right-hand side of the ODE, which is used to calculate the Jacobian matrix using automatic differentiation.

Exceptions

std::runtime_error If there are no species in the network.

Definition at line 558 of file engine_graph.cpp.

◆ reserveJacobianMatrix()

void gridfire::GraphEngine::reserveJacobianMatrix ( )

private

Reserves space for the Jacobian matrix.

This method reserves space for the Jacobian matrix, which is used to store the partial derivatives of the right-hand side of the ODE with respect to the species abundances.

Definition at line 117 of file engine_graph.cpp.

◆ setPrecomputation()

void gridfire::GraphEngine::setPrecomputation ( bool precompute )

Definition at line 369 of file engine_graph.cpp.

◆ setScreeningModel()

void gridfire::GraphEngine::setScreeningModel ( screening::ScreeningType model )

overridevirtual

Set the electron screening model.

Parameters

model The type of screening model to use for reaction rate calculations.

This method allows changing the screening model at runtime. Screening corrections account for the electrostatic shielding of nuclei by electrons, which affects reaction rates in dense stellar plasmas.

Usage Example:: myEngine.setScreeningModel(screening::ScreeningType::WEAK);

gridfire::screening::ScreeningType::WEAK
@ WEAK
Weak screening model (Salpeter, 1954).
Definition screening_types.h:35

Postcondition: The engine will use the specified screening model for subsequent rate calculations.

Implements gridfire::DynamicEngine.

Definition at line 360 of file engine_graph.cpp.

◆ syncInternalMaps()

void gridfire::GraphEngine::syncInternalMaps ( )

private

Synchronizes the internal maps.

This method synchronizes the internal maps used by the engine, including the species map, reaction ID map, and species-to-index map. It also generates the stoichiometry matrix and records the AD tape.

Definition at line 62 of file engine_graph.cpp.

◆ update()

void gridfire::GraphEngine::update ( const NetIn & netIn )

overridevirtual

Update the internal state of the engine.

Parameters

netIn A struct containing the current network input, such as temperature, density, and composition.

This method is intended to be implemented by derived classes to update their internal state based on the provided network conditions. For example, an adaptive engine might use this to re-evaluate which reactions and species are active. For other engines that do not support manually updating, this method might do nothing.

Usage Example:: NetIn input = { ... };

myEngine.update(input);

gridfire::NetIn
Definition network.h:53

Postcondition: The internal state of the engine is updated to reflect the new conditions.

Implements gridfire::DynamicEngine.

Definition at line 554 of file engine_graph.cpp.

◆ validateComposition()

void gridfire::GraphEngine::validateComposition	(	const fourdst::composition::Composition &	composition,
		double	culling,
		double	T9 )

private

Validates the composition against the current reaction set.

Parameters

composition	The composition to validate.
culling	The culling threshold to use.
T9	The temperature to use.

This method validates the composition against the current reaction set. If the composition is not compatible with the reaction set, the reaction set is rebuilt from the composition.

Definition at line 203 of file engine_graph.cpp.

◆ validateConservation()

bool gridfire::GraphEngine::validateConservation ( ) const

nodiscardprivate

Validates mass and charge conservation across all reactions.

Returns: True if all reactions conserve mass and charge, false otherwise.

This method checks that all reactions in the network conserve mass and charge. If any reaction does not conserve mass or charge, an error message is logged and false is returned.

Definition at line 148 of file engine_graph.cpp.

Member Data Documentation

◆ m_config

Config& gridfire::GraphEngine::m_config = Config::getInstance()

private

Definition at line 327 of file engine_graph.h.

◆ m_constants

constants gridfire::GraphEngine::m_constants

private

Definition at line 330 of file engine_graph.h.

◆ m_jacobianMatrix

boost::numeric::ublas::compressed_matrix<double> gridfire::GraphEngine::m_jacobianMatrix

private

Jacobian matrix (species x species).

Definition at line 340 of file engine_graph.h.

◆ m_logger

quill::Logger* gridfire::GraphEngine::m_logger = LogManager::getInstance().getLogger("log")

private

Definition at line 328 of file engine_graph.h.

◆ m_networkSpecies

std::vector<fourdst::atomic::Species> gridfire::GraphEngine::m_networkSpecies

private

Vector of unique species in the network.

Definition at line 335 of file engine_graph.h.

◆ m_networkSpeciesMap

std::unordered_map<std::string_view, fourdst::atomic::Species> gridfire::GraphEngine::m_networkSpeciesMap

private

Map from species name to Species object.

Definition at line 336 of file engine_graph.h.

◆ m_precomputedReactions

std::vector<PrecomputedReaction> gridfire::GraphEngine::m_precomputedReactions

private

Precomputed reactions for efficiency.

Definition at line 349 of file engine_graph.h.

◆ m_reactionIDMap

std::unordered_map<std::string_view, reaction::Reaction*> gridfire::GraphEngine::m_reactionIDMap

private

Map from reaction ID to REACLIBReaction. //PERF: This makes copies of REACLIBReaction and could be a performance bottleneck.

Definition at line 333 of file engine_graph.h.

◆ m_reactions

reaction::LogicalReactionSet gridfire::GraphEngine::m_reactions

private

Set of REACLIB reactions in the network.

Definition at line 332 of file engine_graph.h.

◆ m_rhsADFun

CppAD::ADFun<double> gridfire::GraphEngine::m_rhsADFun

private

CppAD function for the right-hand side of the ODE.

Definition at line 342 of file engine_graph.h.

◆ m_screeningModel

std::unique_ptr<screening::ScreeningModel> gridfire::GraphEngine::m_screeningModel = screening::selectScreeningModel(m_screeningType)

private

Definition at line 345 of file engine_graph.h.

◆ m_screeningType

screening::ScreeningType gridfire::GraphEngine::m_screeningType = screening::ScreeningType::BARE

private

Screening type for the reaction network. Default to no screening.

Definition at line 344 of file engine_graph.h.

◆ m_speciesToIndexMap

std::unordered_map<fourdst::atomic::Species, size_t> gridfire::GraphEngine::m_speciesToIndexMap

private

Map from species to their index in the stoichiometry matrix.

Definition at line 337 of file engine_graph.h.

◆ m_stoichiometryMatrix

boost::numeric::ublas::compressed_matrix<int> gridfire::GraphEngine::m_stoichiometryMatrix

private

Stoichiometry matrix (species x reactions).

Definition at line 339 of file engine_graph.h.

◆ m_usePrecomputation

bool gridfire::GraphEngine::m_usePrecomputation = true

private

Flag to enable or disable using precomputed reactions for efficiency. Mathematically, this should not change the results. Generally end users should not need to change this.

Definition at line 347 of file engine_graph.h.

The documentation for this class was generated from the following files:

src/network/include/gridfire/engine/engine_graph.h
src/network/lib/engine/engine_graph.cpp

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GraphEngine() [1/2]

◆ GraphEngine() [2/2]

Member Function Documentation

◆ calculateAllDerivatives() [1/3]

◆ calculateAllDerivatives() [2/3]

◆ calculateAllDerivatives() [3/3]

◆ calculateAllDerivativesUsingPrecomputation()

◆ calculateMolarReactionFlow() [1/2]

◆ calculateMolarReactionFlow() [2/2]

◆ calculateRHSAndEnergy()

◆ collectNetworkSpecies()

◆ exportToCSV()

◆ exportToDot()

◆ generateJacobianMatrix()

◆ generateStoichiometryMatrix()

◆ getJacobianMatrixEntry()

◆ getNetReactionStoichiometry()

◆ getNetworkReactions()

◆ getNetworkSpecies()

◆ getScreeningModel()

◆ getSpeciesTimescales()

◆ getStoichiometryMatrixEntry()

◆ involvesSpecies()

◆ isPrecomputationEnabled()

◆ populateReactionIDMap()

◆ populateSpeciesToIndexMap()

◆ precomputeNetwork()

◆ recordADTape()

◆ reserveJacobianMatrix()

◆ setPrecomputation()

◆ setScreeningModel()

◆ syncInternalMaps()

◆ update()

◆ validateComposition()

◆ validateConservation()

Member Data Documentation

◆ m_config

◆ m_constants

◆ m_jacobianMatrix

◆ m_logger

◆ m_networkSpecies

◆ m_networkSpeciesMap

◆ m_precomputedReactions

◆ m_reactionIDMap

◆ m_reactions

◆ m_rhsADFun

◆ m_screeningModel

◆ m_screeningType

◆ m_speciesToIndexMap

◆ m_stoichiometryMatrix

◆ m_usePrecomputation