composition.h

Go to the documentation of this file.

/* ***********************************************************************
//
//   Copyright (C) 2025 -- The 4D-STAR Collaboration
//   File Author: Emily Boudreaux
//   Last Modified: March 26, 2025
//
//   4DSSE is free software; you can use it and/or modify
//   it under the terms and restrictions the GNU General Library Public
//   License version 3 (GPLv3) as published by the Free Software Foundation.
//
//   4DSSE is distributed in the hope that it will be useful,
//   but WITHOUT ANY WARRANTY; without even the implied warranty of
//   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
//   See the GNU Library General Public License for more details.
//
//   You should have received a copy of the GNU Library General Public License
//   along with this software; if not, write to the Free Software
//   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
// *********************************************************************** */
#pragma once
 
#include <iostream>
#include <string>
#include <unordered_map>
#include <set>
 
#include <utility>
 
#include "probe.h"
#include "config.h"
 
#include "atomicSpecies.h"
 
namespace serif::composition {
    struct CanonicalComposition {
        double X = 0.0; 
        double Y = 0.0; 
        double Z = 0.0; 
 
        friend std::ostream& operator<<(std::ostream& os, const CanonicalComposition& composition) {
            os << "<CanonicalComposition: "
               << "X = " << composition.X << ", "
               << "Y = " << composition.Y << ", "
               << "Z = " << composition.Z << ">";
            return os;
        }
    };

    struct GlobalComposition {
        double specificNumberDensity; 
        double meanParticleMass; 
 
        // Overload the output stream operator for GlobalComposition
        friend std::ostream& operator<<(std::ostream& os, const GlobalComposition& comp);
    };

    struct CompositionEntry {
        std::string m_symbol; 
        serif::atomic::Species m_isotope; 
        bool m_massFracMode = true; 
 
        double m_massFraction = 0.0; 
        double m_numberFraction = 0.0; 
        double m_relAbundance = 0.0; 
 
        bool m_initialized = false; 

        CompositionEntry();

        CompositionEntry(const std::string& symbol, bool massFracMode=true);

        CompositionEntry(const CompositionEntry& entry);

        void setSpecies(const std::string& symbol);

        std::string symbol() const;

        double mass_fraction() const;

        double mass_fraction(double meanMolarMass) const;

        double number_fraction() const;

        double number_fraction(double totalMoles) const;

        double rel_abundance() const;

        serif::atomic::Species isotope() const;

        bool getMassFracMode() const;

        void setMassFraction(double mass_fraction);

        void setNumberFraction(double number_fraction);

        bool setMassFracMode(double meanMolarMass);

        bool setNumberFracMode(double totalMoles);

        friend std::ostream& operator<<(std::ostream& os, const CompositionEntry& entry);
    };

    class Composition {
    private:
        serif::config::Config& m_config = serif::config::Config::getInstance();
        serif::probe::LogManager& m_logManager = serif::probe::LogManager::getInstance();
        quill::Logger* m_logger = m_logManager.getLogger("log");
 
        bool m_finalized = false; 
        double m_specificNumberDensity = 0.0; 
        double m_meanParticleMass = 0.0; 
        bool m_massFracMode = true; 
 
        std::set<std::string> m_registeredSymbols; 
        std::unordered_map<std::string, CompositionEntry> m_compositions; 

        static bool isValidSymbol(const std::string& symbol);

        bool isValidComposition(const std::vector<double>& fractions) const;

        void validateComposition(const std::vector<double>& fractions) const;

        bool finalizeMassFracMode(bool norm);

        bool finalizeNumberFracMode(bool norm);
 
    public:
        Composition() = default;

        ~Composition() = default;

        bool finalize(bool norm=false);

        explicit Composition(const std::vector<std::string>& symbols);

        explicit Composition(const std::set<std::string>& symbols);

        Composition(const std::vector<std::string>& symbols, const std::vector<double>& mass_fractions, bool massFracMode=true);

        Composition(const Composition& composition);
 
        Composition& operator=(Composition const& other);

        void registerSymbol(const std::string& symbol, bool massFracMode=true);

        void registerSymbol(const std::vector<std::string>& symbols, bool massFracMode=true);

        [[nodiscard]] std::set<std::string> getRegisteredSymbols() const;

        double setMassFraction(const std::string& symbol, const double& mass_fraction);

        std::vector<double> setMassFraction(const std::vector<std::string>& symbols, const std::vector<double>& mass_fractions);

        double setNumberFraction(const std::string& symbol, const double& number_fraction);

        std::vector<double> setNumberFraction(const std::vector<std::string>& symbols, const std::vector<double>& number_fractions);

        Composition mix(const Composition& other, double fraction) const;

        [[nodiscard]] std::unordered_map<std::string, double> getMassFraction() const;

        [[nodiscard]] double getMassFraction(const std::string& symbol) const;

        [[nodiscard]] double getNumberFraction(const std::string& symbol) const;

        [[nodiscard]] std::unordered_map<std::string, double> getNumberFraction() const;

        [[nodiscard]] std::pair<CompositionEntry, GlobalComposition> getComposition(const std::string& symbol) const;

        [[nodiscard]] std::pair<std::unordered_map<std::string, CompositionEntry>, GlobalComposition> getComposition() const;

        [[nodiscard]] double getMeanParticleMass() const;

        [[nodiscard]] double getMeanAtomicNumber() const;

        Composition subset(const std::vector<std::string>& symbols, std::string method="norm") const;

        bool hasSymbol(const std::string& symbol) const;
 
        bool contains(const serif::atomic::Species& isotope) const;

        void setCompositionMode(bool massFracMode);

        [[nodiscard]] CanonicalComposition getCanonicalComposition(bool harsh=false) const;

        friend std::ostream& operator<<(std::ostream& os, const Composition& composition);
 
        // Overload the + operator to call mix with a fraction of 0.5
        Composition operator+(const Composition& other) const;
 
    };
}; // namespace serif::composition