sum.hpp

#ifndef ABLATELIBRARY_RADIATIONPROPERTIESADD_HPP
#define ABLATELIBRARY_RADIATIONPROPERTIESADD_HPP
 
#include <map>
#include <memory>
#include "radiationProperties.hpp"
 
namespace ablate::eos::radiationProperties {
 
class Sum : public RadiationModel {
   private:
    const std::vector<std::shared_ptr<ablate::eos::radiationProperties::RadiationModel>> models;
 
    static PetscErrorCode EmissionTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
    static PetscErrorCode SumTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
   public:
    explicit Sum(std::vector<std::shared_ptr<ablate::eos::radiationProperties::RadiationModel>> models);
    explicit Sum(const Sum&) = delete;
    void operator=(const Sum&) = delete;
 
    [[nodiscard]] ThermodynamicTemperatureFunction GetRadiationPropertiesTemperatureFunction(RadiationProperty property, const std::vector<domain::Field>& fields) const override;
};
 
}  // namespace ablate::eos::radiationProperties
 
#endif  // ABLATELIBRARY_RADIATIONPROPERTIESZIMMER_H