GetTemperature.H

Go to the documentation of this file.

/* Copyright 2021 Hannah Klion
 *
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
 
#ifndef WARPX_GET_TEMPERATURE_H_
#define WARPX_GET_TEMPERATURE_H_
 
#include "TemperatureProperties.H"

struct GetTemperature
{
    /* Type of temperature initialization */
    TemperatureInitType m_type;
 
    /* Constant temperature value, if m_type == TempConstantValue */
    amrex::Real m_temperature;
    /* Temperature parser function, if m_type == TempParserFunction */
    amrex::ParserExecutor<3> m_temperature_parser;

    explicit GetTemperature (TemperatureProperties const& temp) noexcept;

    AMREX_GPU_HOST_DEVICE
    amrex::Real operator() (amrex::Real const x, amrex::Real const y, amrex::Real const z) const noexcept
    {
        switch (m_type)
        {
            case (TempConstantValue):
            {
                return m_temperature;
            }
            case (TempParserFunction):
            {
                return m_temperature_parser(x,y,z);
            }
            default:
            {
                amrex::Abort("Get initial temperature: unknown type");
                return 0.0;
            }
        }
    }
};
struct GetTemperatureVector
{
    /* Type of temperature initialization */
    TemperatureInitType m_type;
 
    amrex::Real m_ux_std{0}, m_uy_std{0}, m_uz_std{0};
    amrex::ParserExecutor<3> m_ux_std_parser, m_uy_std_parser, m_uz_std_parser;
    explicit GetTemperatureVector (TemperatureProperties const& temp) noexcept;

    AMREX_GPU_HOST_DEVICE
    amrex::XDim3 operator() (amrex::Real const x, amrex::Real const y, amrex::Real const z) const noexcept
    {
        switch (m_type)
        {
            case (TempConstantVector):
            {
                return amrex::XDim3{m_ux_std, m_uy_std, m_uz_std};
            }
            case (TempParserFunctionVector):
            {
                return amrex::XDim3{m_ux_std_parser(x,y,z), m_uy_std_parser(x,y,z), m_uz_std_parser(x,y,z)};
            }
            default:
            {
                amrex::Abort("Get initial temperature: unknown type");
                return amrex::XDim3{0.0, 0.0, 0.0};
            }
        }
    }
};
#endif //WARPX_GET_TEMPERATURE_H_