CreateReactionRate.cpp

Go to the documentation of this file.

 
#include <optional>
 
#include "BaseLib/ConfigTree.h"
#include "ChemistryLib/PhreeqcIOData/ReactionRate.h"
#include "ChemistryLib/PhreeqcKernelData/ReactionRate.h"
 
namespace ChemistryLib
{
template <typename ReactionRate>
std::vector<ReactionRate> createReactionRates(
    std::optional<BaseLib::ConfigTree> const& config)
{
    if (!config)
    {
        return {};
    }
 
    std::vector<ReactionRate> reaction_rates;
    for (auto const& rate_config :
         config->getConfigSubtreeList("rate"))
    {
        auto kinetic_reactant =
            rate_config.getConfigParameter<std::string>("kinetic_reactant");
 
        auto const expression_config =
            rate_config.getConfigSubtree("expression");
        auto const& statements =
            expression_config.getConfigParameterList<std::string>("statement");
 
        std::vector<std::string> expression_statements;
        expression_statements.reserve(statements.size());
        std::copy(begin(statements),
                  end(statements),
                  back_inserter(expression_statements));
 
        reaction_rates.emplace_back(std::move(kinetic_reactant),
                                    std::move(expression_statements));
    }
 
    return reaction_rates;
}
 
template std::vector<PhreeqcIOData::ReactionRate>
createReactionRates<PhreeqcIOData::ReactionRate>(
    std::optional<BaseLib::ConfigTree> const& config);
 
template std::vector<PhreeqcKernelData::ReactionRate>
createReactionRates<PhreeqcKernelData::ReactionRate>(
    std::optional<BaseLib::ConfigTree> const& config);
}  // namespace ChemistryLib