CreateChemicalReactionData.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateChemicalReactionData.h"
 
#include <boost/algorithm/string/predicate.hpp>
 
#include "BaseLib/ConfigTree.h"
#include "ChemicalReaction.h"
 
namespace ChemistryLib
{
namespace SelfContainedSolverData
{
std::vector<std::unique_ptr<ChemicalReaction>> createChemicalReactionData(
    BaseLib::ConfigTree const& config)
{
    std::vector<std::unique_ptr<ChemicalReaction>> chemical_reactions;
 
    for (
        auto const& reaction_config :
        config.getConfigSubtreeList("chemical_reaction"))
    {
        auto const stoichiometric_vector =
            reaction_config.getConfigParameter<std::vector<double>>(
                "stoichiometric_coefficients");
 
        auto const reaction_type =
            reaction_config.getConfigParameter<std::string>("reaction_type");
        if (boost::iequals(reaction_type, "FirstOrderReaction"))
        {
            chemical_reactions.emplace_back(std::make_unique<
                                            FirstOrderReaction>(
                stoichiometric_vector,
                reaction_config.getConfigParameter<double>(
                    "first_order_rate_constant")));
        }
    }
 
    return chemical_reactions;
}
}  // namespace SelfContainedSolverData
}  // namespace ChemistryLib