constant.hpp

#ifndef ABLATELIBRARY_TRANSPORT_MODEL_CONSTANT_HPP
#define ABLATELIBRARY_TRANSPORT_MODEL_CONSTANT_HPP
#include "transportModel.hpp"
namespace ablate::eos::transport {
 
class Constant : public TransportModel {
   private:
    const bool active;
    const PetscReal k;
    const PetscReal mu;
    const std::vector<PetscReal> diff;
 
    static PetscErrorCode ConstantFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
 
    static PetscErrorCode ConstantVectorFunction(const PetscReal conserved[], PetscReal* property, void* ctx);
 
    static PetscErrorCode ConstantTemperatureFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
    static PetscErrorCode ConstantTemperatureVectorFunction(const PetscReal conserved[], PetscReal temperature, PetscReal* property, void* ctx);
 
   public:
    explicit Constant(double k = 0, double mu = 0, double diff = 0);
    explicit Constant(double k, double mu, std::vector<double> difft);
    explicit Constant(const Constant&) = delete;
    void operator=(const Constant&) = delete;
 
    [[nodiscard]] ThermodynamicFunction GetTransportFunction(TransportProperty property, const std::vector<domain::Field>& fields) const override;
 
    [[nodiscard]] ThermodynamicTemperatureFunction GetTransportTemperatureFunction(TransportProperty property, const std::vector<domain::Field>& fields) const override;
};
}  // namespace ablate::eos::transport
#endif  // ABLATELIBRARY_CONSTANT_HPP