//! @file Units.cpp

// This file is part of Cantera. See License.txt in the top-level directory or
// at https://cantera.org/license.txt for license and copyright information.

#include "cantera/base/Units.h"
#include "cantera/base/ctexceptions.h"
#include "cantera/base/global.h"
#include "cantera/base/stringUtils.h"
#include "cantera/base/AnyMap.h"
#include "cantera/base/utilities.h"
#include <regex>

namespace {
using namespace Cantera;

const map<string, Units> knownUnits{
    {"", Units(1.0)},
    {"1", Units(1.0)},

    // Mass [M]
    {"kg", Units(1.0, 1, 0, 0)},
    {"g", Units(1e-3, 1, 0, 0)},

    // Length [L]
    {"m", Units(1.0, 0, 1, 0)},
    {"micron", Units(1e-6, 0, 1, 0)},
    {"angstrom", Units(1e-10, 0, 1, 0)},
    {"Å", Units(1e-10, 0, 1, 0)},

    // Time [T]
    {"s", Units(1.0, 0, 0, 1)},
    {"min", Units(60, 0, 0, 1)},
    {"hr", Units(3600, 0, 0, 1)},

    // Temperature [K]
    {"K", Units(1.0, 0, 0, 0, 1)},
    {"C", Units(1.0, 0, 0, 0, 1)},

    // Current [A]
    {"A", Units(1.0, 0, 0, 0, 0, 1)},

    // Quantity [Q]
    {"mol", Units(1e-3, 0, 0, 0, 0, 0, 1)},
    {"gmol", Units(1e-3, 0, 0, 0, 0, 0, 1)},
    {"mole", Units(1e-3, 0, 0, 0, 0, 0, 1)},
    {"kmol", Units(1.0, 0, 0, 0, 0, 0, 1)},
    {"kgmol", Units(1.0, 0, 0, 0, 0, 0, 1)},
    {"molec", Units(1.0/Avogadro, 0, 0, 0, 0, 0, 1)},

    // Energy [M*L^2/T^2]
    {"J", Units(1.0, 1, 2, -2)},
    {"cal", Units(4.184, 1, 2, -2)},
    {"erg", Units(1e-7, 1, 2, -2)},
    {"eV", Units(ElectronCharge, 1, 2, -2)},

    // Force [M*L/T^2]
    {"N", Units(1.0, 1, 1, -2)},
    {"dyn", Units(1e-5, 1, 1, -2)},

    // Pressure [M/L/T^2]
    {"Pa", Units(1.0, 1, -1, -2)},
    {"atm", Units(OneAtm, 1, -1, -2)},
    {"bar", Units(1.0e5, 1, -1, -2)},
    {"dyn/cm^2", Units(0.1, 1, -1, -2)},

    // Volume [L^3]
    {"m^3", Units(1.0, 0, 3, 0)},
    {"m³", Units(1.0, 0, 3, 0)},
    {"liter", Units(0.001, 0, 3, 0)},
    {"L", Units(0.001, 0, 3, 0)},
    {"l", Units(0.001, 0, 3, 0)},
    {"cc", Units(1.0e-6, 0, 3, 0)},

    // Other electrical units
    {"ohm", Units(1.0, 1, 2, -3, 0, -2)}, // kg*m^2/s^3/A^2
    {"V", Units(1.0, 1, 2, -3, 0, -1)}, // kg*m^2/s^3/A
    {"coulomb", Units(1.0, 0, 0, 1, 0, 1)}, // A*s

    //! Activation energy units [M*L^2/T^2/Q]
    {"J/kmol", Units(1.0, 1, 2, -2, 0, 0, -1)},
};

const map<string, double> prefixes{
    {"Y", 1e24},
    {"Z", 1e21},
    {"E", 1e18},
    {"P", 1e15},
    {"T", 1e12},
    {"G", 1e9},
    {"M", 1e6},
    {"k", 1e3},
    {"h", 1e2},
    {"d", 1e-1},
    {"c", 1e-2},
    {"m", 1e-3},
    {"u", 1e-6},
    {"n", 1e-9},
    {"p", 1e-12},
    {"f", 1e-15},
    {"a", 1e-18},
    {"z", 1e-21},
    {"y", 1e-24}
};
}

namespace Cantera
{

Units::Units(double factor, double mass, double length, double time,
             double temperature, double current, double quantity)
    : m_factor(factor)
    , m_mass_dim(mass)
    , m_length_dim(length)
    , m_time_dim(time)
    , m_temperature_dim(temperature)
    , m_current_dim(current)
    , m_quantity_dim(quantity)
{
    if (mass != 0 && length == -mass && time == -2 * mass
        && temperature == 0 && current == 0 && quantity == 0) {
        // Dimension looks like Pa^n
        m_pressure_dim = mass;
    } else if (mass != 0 && length == 2 * mass && time == -2 * mass
               && temperature == 0 && current == 0 && quantity == 0)
    {
        // Dimension looks like J^n
        m_energy_dim = mass;
    }
}

Units::Units(const string& name, bool force_unity)
{
    size_t start = 0;

    // Determine factor
    static std::regex regexp("[+-]?([0-9]+([.][0-9]*)?|[.][0-9]+)");
    std::smatch matched;
    std::regex_search(name, matched, regexp);
    if (matched.size()) {
        string factor = *matched.begin();
        if (name.find(factor) == 0) {
            m_factor = fpValueCheck(factor);
            start = factor.size();
        }
    }

    while (true) {
        // Split into groups of the form 'unit^exponent'
        size_t stop = name.find_first_of("*/", start);
        size_t caret = name.find('^', start);
        if (caret > stop) {
            // No caret in this group
            caret = npos;
        }
        string unit = trimCopy(
            name.substr(start, std::min(caret, stop) - start));

        double exponent = 1.0;
        if (caret != npos) {
            exponent = fpValueCheck(name.substr(caret+1, stop-caret-1));
        }
        if (start != 0 && name[start-1] == '/') {
            // This unit is in the denominator
            exponent = -exponent;
        }

        if (knownUnits.find(unit) != knownUnits.end()) {
            // Incorporate the unit defined by the current group
            *this *= knownUnits.at(unit).pow(exponent);
        } else {
            // See if the unit looks like a prefix + base unit
            string prefix = unit.substr(0, 1);
            string suffix = unit.substr(1);
            if (prefixes.find(prefix) != prefixes.end() &&
                knownUnits.find(suffix) != knownUnits.end()) {
                Units u = knownUnits.at(suffix);
                u.scale(prefixes.at(prefix));
                *this *= u.pow(exponent);
            } else {
                throw CanteraError("Units::Units(string)",
                    "Unknown unit '{}' in unit string '{}'", unit, name);
            }
        }

        start = stop+1;
        if (stop == npos) {
            break;
        }
    }

    if (force_unity && (std::abs(m_factor - 1.) > SmallNumber)) {
        throw CanteraError("Units::Units(string)",
            "Detected non-unity conversion factor:\n"
            "input '{}' is equivalent to '{}' where the conversion factor is '{}'",
            name, str(), m_factor);
    }
}

bool Units::convertible(const Units& other) const
{
    return (m_mass_dim == other.m_mass_dim &&
            m_length_dim == other.m_length_dim &&
            m_time_dim == other.m_time_dim &&
            m_temperature_dim == other.m_temperature_dim &&
            m_current_dim == other.m_current_dim &&
            m_quantity_dim == other.m_quantity_dim);
}

Units& Units::operator*=(const Units& other)
{
    m_factor *= other.m_factor;
    m_mass_dim += other.m_mass_dim;
    m_length_dim += other.m_length_dim;
    m_time_dim += other.m_time_dim;
    m_temperature_dim += other.m_temperature_dim;
    m_current_dim += other.m_current_dim;
    m_quantity_dim += other.m_quantity_dim;
    m_pressure_dim += other.m_pressure_dim;
    m_energy_dim += other.m_energy_dim;
    return *this;
}

Units Units::pow(double exponent) const
{
    return Units(std::pow(m_factor, exponent),
                 m_mass_dim * exponent,
                 m_length_dim * exponent,
                 m_time_dim * exponent,
                 m_temperature_dim * exponent,
                 m_current_dim * exponent,
                 m_quantity_dim * exponent);
}

string Units::str(bool skip_unity) const
{
    map<string, double> dims{
        {"kg", m_mass_dim},
        {"m", m_length_dim},
        {"s", m_time_dim},
        {"K", m_temperature_dim},
        {"A", m_current_dim},
        {"kmol", m_quantity_dim},
    };

    string num = "";
    string den = "";
    for (auto const& [dimension, exponent] : dims) {
        int rounded = (int)round(exponent);
        if (exponent == 0.) {
            // skip
        } else if (exponent == 1.) {
            num.append(fmt::format(" * {}", dimension));
        } else if (exponent == -1.) {
            den.append(fmt::format(" / {}", dimension));
        } else if (exponent == rounded && rounded > 0) {
            num.append(fmt::format(" * {}^{}", dimension, rounded));
        } else if (exponent == rounded) {
            den.append(fmt::format(" / {}^{}", dimension, -rounded));
        } else if (exponent > 0) {
            num.append(fmt::format(" * {}^{}", dimension, exponent));
        } else {
            den.append(fmt::format(" / {}^{}", dimension, -exponent));
        }
    }

    if (skip_unity && (std::abs(m_factor - 1.) < SmallNumber)) {
        if (num.size()) {
            return fmt::format("{}{}", num.substr(3), den);
        }
        // print '1' as the numerator is empty
        return fmt::format("1{}", den);
    }

    string factor;
    if (m_factor == round(m_factor)) {
        // ensure that fmt::format does not round to integer
        factor = fmt::format("{:.1f}", m_factor);
    } else {
        factor = fmt::format("{}", m_factor);
    }

    if (num.size()) {
        // concatenate factor, numerator and denominator (skipping leading '*')
        return fmt::format("{} {}{}", factor, num.substr(3), den);
    }

    return fmt::format("{}{}", factor, den);
}

bool Units::operator==(const Units& other) const
{
    return m_factor == other.m_factor
           && m_mass_dim == other.m_mass_dim
           && m_length_dim == other.m_length_dim
           && m_time_dim == other.m_time_dim
           && m_temperature_dim == other.m_temperature_dim
           && m_current_dim == other.m_current_dim
           && m_quantity_dim == other.m_quantity_dim
           && m_pressure_dim == other.m_pressure_dim
           && m_energy_dim == other.m_energy_dim;
}

double Units::dimension(const string& primary) const
{
    if (primary == "mass") {
        return m_mass_dim;
    } else if (primary == "length") {
        return m_length_dim;
    } else if (primary == "time") {
        return m_time_dim;
    } else if (primary == "temperature") {
        return m_temperature_dim;
    } else if (primary == "current") {
        return m_current_dim;
    } else if (primary == "quantity") {
        return m_quantity_dim;
    } else {
        throw CanteraError("Units::dimension",
            "Unknown primary unit '{}'.", primary);
    }
}

Units UnitStack::standardUnits() const
{
    if (!stack.size()) {
        return Units(0);
    }
    return stack[0].first;
}

void UnitStack::setStandardUnits(Units& standardUnits)
{
    if (!stack.size()) {
        stack.emplace_back(standardUnits, 0.);
        return;
    }
    if (stack[0].second != 0.) {
        throw CanteraError("UnitStack::setStandardUnit",
            "Standard unit is already defined.");
    }
    stack[0].first = standardUnits;
}

double UnitStack::standardExponent() const
{
    if (!stack.size()) {
        return NAN;
    }
    return stack[0].second;
}

void UnitStack::join(double exponent)
{
    if (!stack.size()) {
        throw CanteraError("UnitStack::join",
            "Standard unit is not defined.");
    }
    stack[0].second += exponent;

}

void UnitStack::update(const Units& units, double exponent)
{
    bool found = false;
    for (auto& [current_unit, current_exp] : stack) {
        if (current_unit == units) {
            current_exp += exponent;
            found = true;
            break;
        }
    }
    if (!found) {
        stack.emplace_back(units, exponent);
    }
}

Units UnitStack::product() const
{
    if (!stack.size()) {
        return Units(0.);
    }
    Units out = Units(1.);
    for (auto& [units, exponent] : stack) {
        if (exponent == 1) {
            out *= units;
        } else {
            out *= units.pow(exponent);
        }
    }
    return out;
}

UnitSystem::UnitSystem(std::initializer_list<string> units)
{
    setDefaults(units);
}

map<string, string> UnitSystem::defaults() const
{
    // Unit system defaults
    map<string, string> units{
        {"mass", "kg"},
        {"length", "m"},
        {"time", "s"},
        {"quantity", "kmol"},
        {"pressure", "Pa"},
        {"energy", "J"},
        {"temperature", "K"},
        {"current", "A"},
        {"activation-energy", "J / kmol"},
    };

    // Overwrite entries that have conversion factors
    for (const auto& [dimension, default_unit] : m_defaults) {
        units[dimension] = default_unit;
    }

    // Activation energy follows specified energy and quantity units
    // unless given explicitly
    if (!m_defaults.count("activation-energy")) {
        units["activation-energy"] = fmt::format(
            "{} / {}", units["energy"], units["quantity"]);
    }

    return units;
}

void UnitSystem::setDefaults(std::initializer_list<string> units)
{
    for (const auto& name : units) {
        auto unit = Units(name);
        if (unit.convertible(knownUnits.at("kg"))) {
            m_mass_factor = unit.factor();
            m_defaults["mass"] = name;
        } else if (unit.convertible(knownUnits.at("m"))) {
            m_length_factor = unit.factor();
            m_defaults["length"] = name;
        } else if (unit.convertible(knownUnits.at("s"))) {
            m_time_factor = unit.factor();
            m_defaults["time"] = name;
        } else if (unit.convertible(knownUnits.at("kmol"))) {
            m_quantity_factor = unit.factor();
            m_defaults["quantity"] = name;
        } else if (unit.convertible(knownUnits.at("Pa"))) {
            m_pressure_factor = unit.factor();
            m_defaults["pressure"] = name;
        } else if (unit.convertible(knownUnits.at("J"))) {
            m_energy_factor = unit.factor();
            m_defaults["energy"] = name;
        } else if (unit.convertible(knownUnits.at("K"))) {
            // Do nothing -- no other scales are supported for temperature
            if (unit.factor() != 1.) {
                throw CanteraError("UnitSystem::setDefaults", "Temperature scales "
                    "with non-unity conversion factor from Kelvin are not supported.");
            }
        } else if (unit.convertible(knownUnits.at("A"))) {
            // Do nothing -- no other scales are supported for current
            if (unit.factor() != 1.) {
                throw CanteraError("UnitSystem::setDefaults", "Current scales "
                    "with non-unity conversion factor from Ampere are not supported.");
            }
        } else {
            throw CanteraError("UnitSystem::setDefaults",
                "Unable to match unit '{}' to a basic dimension", name);
        }
    }
    if (!m_explicit_activation_energy) {
        m_activation_energy_factor = m_energy_factor / m_quantity_factor;
    }
}

void UnitSystem::setDefaults(const map<string, string>& units)
{
    for (const auto& [dimension, name] : units) {
        Units unit(name);
        if (dimension == "mass" && unit.convertible(knownUnits.at("kg"))) {
            m_mass_factor = unit.factor();
            m_defaults["mass"] = name;
        } else if (dimension == "length" && unit.convertible(knownUnits.at("m"))) {
            m_length_factor = unit.factor();
            m_defaults["length"] = name;
        } else if (dimension == "time" && unit.convertible(knownUnits.at("s"))) {
            m_time_factor = unit.factor();
            m_defaults["time"] = name;
        } else if (dimension == "temperature" && unit.convertible(knownUnits.at("K"))) {
            // do nothing - no other temperature scales are supported
            if (unit.factor() != 1.) {
                throw CanteraError("UnitSystem::setDefaults", "Temperature scales "
                    "with non-unity conversion factor from Kelvin are not supported.");
            }
        } else if (dimension == "current" && unit.convertible(knownUnits.at("A"))) {
            // do nothing - no other current scales are supported
            if (unit.factor() != 1.) {
                throw CanteraError("UnitSystem::setDefaults", "Current scales "
                    "with non-unity conversion factor from Ampere are not supported.");
            }
        } else if (dimension == "quantity" && unit.convertible(knownUnits.at("kmol"))) {
            m_quantity_factor = unit.factor();
            m_defaults["quantity"] = name;
        } else if (dimension == "pressure" && unit.convertible(knownUnits.at("Pa"))) {
            m_pressure_factor = unit.factor();
            m_defaults["pressure"] = name;
        } else if (dimension == "energy" && unit.convertible(knownUnits.at("J"))) {
            m_energy_factor = unit.factor();
            m_defaults["energy"] = name;
        } else if (dimension == "activation-energy") {
            // handled separately to allow override
        } else {
            throw CanteraError("UnitSystem::setDefaults",
                "Unable to set default unit for '{}' to '{}' ({}).",
                dimension, name, unit.str());
        }
    }
    if (units.find("activation-energy") != units.end()) {
        setDefaultActivationEnergy(units.at("activation-energy"));
    } else if (!m_explicit_activation_energy) {
        m_activation_energy_factor = m_energy_factor / m_quantity_factor;
    }
}

void UnitSystem::setDefaultActivationEnergy(const string& e_units)
{
    Units u(e_units);
    m_defaults["activation-energy"] = e_units;
    if (u.convertible(Units("J/kmol"))) {
        m_activation_energy_factor = u.factor();
    } else if (u.convertible(knownUnits.at("K"))) {
        m_activation_energy_factor = GasConstant;
    } else if (u.convertible(knownUnits.at("eV"))) {
        m_activation_energy_factor = u.factor() * Avogadro;
    } else {
        throw CanteraError("Units::setDefaultActivationEnergy",
            "Unable to match unit '{}' to a unit of activation energy", e_units);
    }
    m_explicit_activation_energy = true;
}

double UnitSystem::convert(double value, const string& src, const string& dest) const
{
    return convert(value, Units(src), Units(dest));
}

double UnitSystem::convert(double value, const Units& src,
                           const Units& dest) const
{
    if (!src.convertible(dest)) {
        throw CanteraError("UnitSystem::convert",
            "Incompatible units:\n    Units({}) and\n    Units({})",
            src.str(), dest.str());
    }
    return value * src.factor() / dest.factor();
}

double UnitSystem::convertTo(double value, const string& dest) const
{
    return convertTo(value, Units(dest));
}

double UnitSystem::convertTo(double value, const Units& dest) const
{
    return value / dest.factor()
        * pow(m_mass_factor, dest.m_mass_dim - dest.m_pressure_dim - dest.m_energy_dim)
        * pow(m_length_factor, dest.m_length_dim + dest.m_pressure_dim - 2*dest.m_energy_dim)
        * pow(m_time_factor, dest.m_time_dim + 2*dest.m_pressure_dim + 2*dest.m_energy_dim)
        * pow(m_quantity_factor, dest.m_quantity_dim)
        * pow(m_pressure_factor, dest.m_pressure_dim)
        * pow(m_energy_factor, dest.m_energy_dim);
}

double UnitSystem::convertFrom(double value, const string& dest) const
{
    return convertFrom(value, Units(dest));
}

double UnitSystem::convertFrom(double value, const Units& src) const
{
    return value * src.factor()
        * pow(m_mass_factor, -src.m_mass_dim + src.m_pressure_dim + src.m_energy_dim)
        * pow(m_length_factor, -src.m_length_dim - src.m_pressure_dim + 2*src.m_energy_dim)
        * pow(m_time_factor, -src.m_time_dim - 2*src.m_pressure_dim - 2*src.m_energy_dim)
        * pow(m_quantity_factor, -src.m_quantity_dim)
        * pow(m_pressure_factor, -src.m_pressure_dim)
        * pow(m_energy_factor, -src.m_energy_dim);
}

static pair<double, string> split_unit(const AnyValue& v) {
    if (v.is<string>()) {
        // Should be a value and units, separated by a space, for example '2e4 J/kmol'
        string val_units = v.asString();
        size_t space = val_units.find(" ");
        if (space == npos) {
            throw CanteraError("split_unit (UnitSystem)",
                "Couldn't parse '{}' as a space-separated value/unit pair\n",
                val_units);
        }
        return {fpValueCheck(val_units.substr(0, space)),
                val_units.substr(space+1)};
    } else {
        // Just a value
        return {v.asDouble(), ""};
    }
}

double UnitSystem::convert(const AnyValue& v, const string& dest) const
{
    try {
        return convert(v, Units(dest));
    } catch (InputFileError&) {
        throw; // already have input file context from convert(AnyValue, Units)
    } catch (CanteraError& err) {
        throw InputFileError("UnitSystem::convert", v, err.getMessage());
    }
}

double UnitSystem::convert(const AnyValue& v, const Units& dest) const
{
    try {
        auto [value, units] = split_unit(v);
        if (units.empty()) {
            // Just a value, so convert using default units
            return convertTo(value, dest);
        } else {
            // Both source and destination units are explicit
            return convert(value, Units(units), dest);
        }
    } catch (CanteraError& err) {
        throw InputFileError("UnitSystem::convert", v, err.getMessage());
    }
}

double UnitSystem::convertRateCoeff(const AnyValue& v, const Units& dest) const
{
    if (dest.factor() != 0) {
        // If the destination units are defined, no special handling is required
        return convert(v, dest);
    }

    auto [value, units] = split_unit(v);
    if (units.empty()) {
        if (m_length_factor == 1.0 && m_quantity_factor == 1.0) {
            // Input is a number in the default mks+kmol system, so no conversion is
            // required
            return value;
        }
    } else {
        Units sourceUnits(units);
        if (fabs(sourceUnits.factor() - 1.0) < 1e-14) {
            // Input is explicitly in the mks+kmol system, so no conversion is required
            return value;
        }
    }
    throw InputFileError("UnitSystem::convertRateCoeff", v,
        "Unable to convert value with non-default units to undefined units,\n"
        "likely while creating a standalone ReactionRate object.");
}

vector<double> UnitSystem::convert(const vector<AnyValue>& vals,
                                   const string& dest) const
{
    return convert(vals, Units(dest));
}

vector<double> UnitSystem::convert(const vector<AnyValue>& vals,
                                   const Units& dest) const
{
    vector<double> out;
    for (const auto& val : vals) {
        out.emplace_back(convert(val, dest));
    }
    return out;
}

double UnitSystem::convertActivationEnergy(double value, const string& src,
                                           const string& dest) const
{
    // Convert to J/kmol
    Units usrc(src);
    if (usrc.convertible(Units("J/kmol"))) {
        value *= usrc.factor();
    } else if (usrc.convertible(Units("K"))) {
        value *= GasConstant * usrc.factor();
    } else if (usrc.convertible(Units("eV"))) {
        value *= Avogadro * usrc.factor();
    } else {
        throw CanteraError("UnitSystem::convertActivationEnergy",
            "Don't understand units '{}' as an activation energy", src);
    }

    // Convert from J/kmol
    Units udest(dest);
    if (udest.convertible(Units("J/kmol"))) {
        value /= udest.factor();
    } else if (udest.convertible(Units("K"))) {
        value /= GasConstant * udest.factor();
    } else if (udest.convertible(Units("eV"))) {
        value /= Avogadro * udest.factor();
    } else {
        throw CanteraError("UnitSystem::convertActivationEnergy",
            "Don't understand units '{}' as an activation energy", dest);
    }

    return value;
}

double UnitSystem::convertActivationEnergyTo(double value, const string& dest) const
{
    return convertActivationEnergyTo(value, Units(dest));
}

double UnitSystem::convertActivationEnergyTo(double value,
                                             const Units& dest) const
{
    if (dest.convertible(Units("J/kmol"))) {
        return value * m_activation_energy_factor / dest.factor();
    } else if (dest.convertible(knownUnits.at("K"))) {
        return value * m_activation_energy_factor / GasConstant;
    } else if (dest.convertible(knownUnits.at("eV"))) {
        return value * m_activation_energy_factor / (Avogadro * dest.factor());
    } else {
        throw CanteraError("UnitSystem::convertActivationEnergyTo",
            "'{}' is not a unit of activation energy", dest.str());
    }
}

double UnitSystem::convertActivationEnergyFrom(double value, const string& src) const
{
    Units usrc(src);
    if (usrc.convertible(Units("J/kmol"))) {
        return value * usrc.factor() / m_activation_energy_factor;
    } else if (usrc.convertible(knownUnits.at("K"))) {
        return value * GasConstant / m_activation_energy_factor;
    } else if (usrc.convertible(knownUnits.at("eV"))) {
        return value * Avogadro * usrc.factor() / m_activation_energy_factor;
    } else {
        throw CanteraError("UnitSystem::convertActivationEnergyFrom",
            "'{}' is not a unit of activation energy", src);
    }
}

double UnitSystem::convertActivationEnergy(const AnyValue& v, const string& dest) const
{
    try {
        auto [value, units] = split_unit(v);
        if (units.empty()) {
            // Just a value, so convert using default units
            return convertActivationEnergyTo(value, dest);
        } else {
            // Both source and destination units are explicit
            return convertActivationEnergy(value, units, dest);
        }
    } catch (CanteraError& err) {
        throw InputFileError(
            "UnitSystem::convertActivationEnergy", v, err.getMessage());
    }
}

AnyMap UnitSystem::getDelta(const UnitSystem& other) const
{
    AnyMap delta;
    // Create a local alias because the template specialization can't be deduced
    // automatically
    const auto& get = getValue<string, string>;
    if (m_mass_factor != other.m_mass_factor) {
        delta["mass"] = get(m_defaults, "mass", "kg");
    }
    if (m_length_factor != other.m_length_factor) {
        delta["length"] = get(m_defaults, "length", "m");
    }
    if (m_time_factor != other.m_time_factor) {
        delta["time"] = get(m_defaults, "time", "s");
    }
    if (m_pressure_factor != other.m_pressure_factor) {
        delta["pressure"] = get(m_defaults, "pressure", "Pa");
    }
    if (m_energy_factor != other.m_energy_factor) {
        delta["energy"] = get(m_defaults, "energy", "J");
    }
    if (m_quantity_factor != other.m_quantity_factor) {
        delta["quantity"] = get(m_defaults, "quantity", "kmol");
    }
    if (m_explicit_activation_energy
        || (other.m_explicit_activation_energy
            && m_activation_energy_factor != m_energy_factor / m_quantity_factor))
    {
        delta["activation-energy"] = get(m_defaults, "activation-energy", "J/kmol");
    }
    return delta;
}

}
