#pragma once

#include "./parallelization/mpi_structure.hpp"

#include <iostream>

#include <sstream>

#include <string>

#include <vector>

#include <array>

#include <map>

#include <cstdlib>

#include <algorithm>

#include <cassert>

template <typename T, typename U>

struct CEmptyMap {

CEmptyMap(std::initializer_list<std::pair<const T, U> >) {}

};

#ifdef ENABLE_MAPS

template<class T, class U>

using MapType = std::map<T,U>;

#define MakePair(a,b) {a,b},

#else

template<class T, class U>

using MapType = CEmptyMap<T,U>;

#define MakePair(a,b)

#endif

enum class SU2_COMPONENT {

};

const unsigned int EXIT_DIVERGENCE = 2; /*!< \brief Exit code (divergence). */

const unsigned int MAX_PARAMETERS = 10; /*!< \brief Maximum number of parameters for a design variable definition. */

const unsigned int MAX_NUMBER_PERIODIC = 10; /*!< \brief Maximum number of periodic boundary conditions. */

const unsigned int MAX_STRING_SIZE = 400; /*!< \brief Maximum size of a generic string. */

const unsigned int MAX_NUMBER_FFD = 15; /*!< \brief Maximum number of FFDBoxes for the FFD. */

enum: unsigned int{MAX_SOLS = 13}; /*!< \brief Maximum number of solutions at the same time (dimension of solution container array). */

const unsigned int MAX_TERMS = 7; /*!< \brief Maximum number of terms in the numerical equations (dimension of solver container array). */

const unsigned int MAX_ZONES = 3; /*!< \brief Maximum number of zones. */

const unsigned int MAX_FE_KINDS = 7; /*!< \brief Maximum number of Finite Elements. */

const unsigned int NO_RK_ITER = 0; /*!< \brief No Runge-Kutta iteration. */

const unsigned int OVERHEAD = 4; /*!< \brief Overhead space above nMarker when allocating space for boundary elems (MPI + periodic). */

const unsigned int MESH_0 = 0; /*!< \brief Definition of the finest grid level. */

const unsigned int MESH_1 = 1; /*!< \brief Definition of the finest grid level. */

const unsigned int ZONE_0 = 0; /*!< \brief Definition of the first grid domain. */

const unsigned int ZONE_1 = 1; /*!< \brief Definition of the second grid domain. */

const unsigned int INST_0 = 0; /*!< \brief Definition of the first instance per grid level. */

constexpr passivedouble STANDARD_GRAVITY = 9.80665; /*!< \brief Acceleration due to gravity at surface of earth. */

constexpr passivedouble UNIVERSAL_GAS_CONSTANT = 8.3144598; /*!< \brief Universal gas constant in J/(mol*K) */

constexpr passivedouble BOLTZMANN_CONSTANT = 1.3806503E-23; /*!< \brief Boltzmann's constant [J K^-1] */

constexpr passivedouble AVOGAD_CONSTANT = 6.0221415E26; /*!< \brief Avogadro's constant, number of particles in one kmole. */

constexpr passivedouble FUND_ELEC_CHARGE_CGS = 4.8032047E-10; /*!< \brief Fundamental electric charge in CGS units, cm^(3/2) g^(1/2) s^(-1). */

constexpr passivedouble STD_REF_TEMP = 298.15; /*!< \brief Standard reference temperature for enthalpy in Kelvin. */

constexpr passivedouble EPS = 1.0E-16; /*!< \brief Error scale. */

constexpr passivedouble TURB_EPS = 1.0E-16; /*!< \brief Turbulent Error scale. */

constexpr passivedouble ONE2 = 0.5; /*!< \brief One divided by two. */

constexpr passivedouble ONE3 = 1.0 / 3.0; /*!< \brief One divided by three. */

constexpr passivedouble TWO3 = 2.0 / 3.0; /*!< \brief Two divided by three. */

constexpr passivedouble FOUR3 = 4.0 / 3.0; /*!< \brief Four divided by three. */

constexpr passivedouble PI_NUMBER = 3.14159265358979323846; /*!< \brief Pi number (not using M_PI to avoid Windows issues). */

constexpr passivedouble STEFAN_BOLTZMANN = 5.670367E-08; /*!< \brief Stefan-Boltzmann constant in W/(m^2*K^4). */

const int MASTER_NODE = 0; /*!< \brief Master node for MPI parallelization. */

const int SINGLE_NODE = 1; /*!< \brief There is only a node in the MPI parallelization. */

const int SINGLE_ZONE = 1; /*!< \brief There is only a zone. */

const unsigned short COMM_TYPE_UNSIGNED_LONG = 1; /*!< \brief Communication type for unsigned long. */

const unsigned short COMM_TYPE_LONG = 2; /*!< \brief Communication type for long. */

const unsigned short COMM_TYPE_UNSIGNED_SHORT = 3; /*!< \brief Communication type for unsigned short. */

const unsigned short COMM_TYPE_DOUBLE = 4; /*!< \brief Communication type for double. */

const unsigned short COMM_TYPE_CHAR = 5; /*!< \brief Communication type for char. */

const unsigned short COMM_TYPE_SHORT = 6; /*!< \brief Communication type for short. */

const unsigned short COMM_TYPE_INT = 7; /*!< \brief Communication type for int. */

enum GEO_TYPE {

};

constexpr unsigned short N_ELEM_TYPES = 7; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_LINE = 2; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_TRIANGLE = 3; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_QUADRILATERAL = 4; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_TETRAHEDRON = 4; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_HEXAHEDRON = 8; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_PYRAMID = 5; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_PRISM = 6; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_POINTS_MAXIMUM = 8; /*!< \brief Max. out of the above, used for static arrays, keep it up to date. */

constexpr unsigned short N_FACES_LINE = 1; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_TRIANGLE = 3; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_QUADRILATERAL = 4; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_TETRAHEDRON = 4; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_PYRAMID = 5; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_PRISM = 5; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_HEXAHEDRON = 6; /*!< \brief General output & CGNS defines. */

constexpr unsigned short N_FACES_MAXIMUM = 6; /*!< \brief Max. out of the above, used for static arrays, keep it up to date. */

inline unsigned short nFacesOfElementType(unsigned short elementType) {

}

inline unsigned short nPointsOfElementType(unsigned short elementType) {

}

const int CGNS_STRING_SIZE = 33; /*!< \brief Length of strings used in the CGNS format. */

const int SU2_CONN_SIZE = 10; /*!< \brief Size of the connectivity array that is allocated for each element

const int SU2_CONN_SKIP = 2; /*!< \brief Offset to skip the globalID and VTK type at the start of the element connectivity list for each CGNS element. */

constexpr passivedouble COLORING_EFF_THRESH = 0.875; /*!< \brief Below this value fallback strategies are used instead. */

constexpr int N_POLY_COEFFS = 5; /*!< \brief Number of coefficients in temperature polynomial fits for fluid models. */

enum ANSWER {

};

enum AVERAGE_TYPE {

AVERAGE_AREA = 1, /*!< \brief Area-weighted average. */

AVERAGE_MASSFLUX = 2 /*!< \brief Mass-flux weighted average. */

};

static const MapType<std::string, AVERAGE_TYPE> Average_Map = {

MakePair("AREA", AVERAGE_AREA)

MakePair("MASSFLUX", AVERAGE_MASSFLUX)

};

enum class MAIN_SOLVER {

};

static const MapType<std::string, MAIN_SOLVER> Solver_Map = {

MakePair("NONE", MAIN_SOLVER::NONE)

MakePair("EULER", MAIN_SOLVER::EULER)

MakePair("NAVIER_STOKES", MAIN_SOLVER::NAVIER_STOKES)

MakePair("RANS", MAIN_SOLVER::RANS)

MakePair("INC_EULER", MAIN_SOLVER::INC_EULER)

MakePair("INC_NAVIER_STOKES", MAIN_SOLVER::INC_NAVIER_STOKES)

MakePair("INC_RANS", MAIN_SOLVER::INC_RANS)

MakePair("FEM_EULER", MAIN_SOLVER::FEM_EULER)

MakePair("FEM_NAVIER_STOKES", MAIN_SOLVER::FEM_NAVIER_STOKES)

MakePair("FEM_RANS", MAIN_SOLVER::FEM_RANS)

MakePair("FEM_LES", MAIN_SOLVER::FEM_LES)

MakePair("NEMO_EULER",MAIN_SOLVER::NEMO_EULER)

MakePair("NEMO_NAVIER_STOKES",MAIN_SOLVER::NEMO_NAVIER_STOKES)

MakePair("HEAT_EQUATION", MAIN_SOLVER::HEAT_EQUATION)

MakePair("ELASTICITY", MAIN_SOLVER::FEM_ELASTICITY)

MakePair("TEMPLATE_SOLVER", MAIN_SOLVER::TEMPLATE_SOLVER)

MakePair("MULTIPHYSICS", MAIN_SOLVER::MULTIPHYSICS)

};

enum class ENUM_MULTIZONE {

MZ_BLOCK_GAUSS_SEIDEL, /*!< \brief Definition of a Block-Gauss-Seidel multizone solver. */

MZ_BLOCK_JACOBI, /*!< \brief Definition of a Block-Jacobi solver. */

};

static const MapType<std::string, ENUM_MULTIZONE> Multizone_Map = {

MakePair("BLOCK_GAUSS_SEIDEL", ENUM_MULTIZONE::MZ_BLOCK_GAUSS_SEIDEL)

MakePair("BLOCK_JACOBI", ENUM_MULTIZONE::MZ_BLOCK_JACOBI)

};

enum class STRUCT_DEFORMATION {

SMALL, /*!< \brief Definition of linear elastic material. */

LARGE, /*!< \brief Definition of Neo-Hookean material. */

};

static const MapType<std::string, STRUCT_DEFORMATION> Struct_Map = {

MakePair("SMALL_DEFORMATIONS", STRUCT_DEFORMATION::SMALL)

MakePair("LARGE_DEFORMATIONS", STRUCT_DEFORMATION::LARGE)

};

enum class STRUCT_MODEL {

};

static const MapType<std::string, STRUCT_MODEL> Material_Map = {

MakePair("LINEAR_ELASTIC", STRUCT_MODEL::LINEAR_ELASTIC)

MakePair("NEO_HOOKEAN", STRUCT_MODEL::NEO_HOOKEAN)

MakePair("KNOWLES", STRUCT_MODEL::KNOWLES)

MakePair("IDEAL_DE", STRUCT_MODEL::IDEAL_DE)

};

enum class STRUCT_COMPRESS {

COMPRESSIBLE, /*!< \brief Definition of compressible material. */

NEARLY_INCOMP, /*!< \brief Definition of nearly incompressible material. */

};

static const MapType<std::string, STRUCT_COMPRESS> MatComp_Map = {

MakePair("COMPRESSIBLE", STRUCT_COMPRESS::COMPRESSIBLE)

MakePair("NEARLY_INCOMPRESSIBLE", STRUCT_COMPRESS::NEARLY_INCOMP)

};

enum class INTERFACE_INTERPOLATOR {

};

static const MapType<std::string, INTERFACE_INTERPOLATOR> Interpolator_Map = {

MakePair("NEAREST_NEIGHBOR", INTERFACE_INTERPOLATOR::NEAREST_NEIGHBOR)

MakePair("ISOPARAMETRIC", INTERFACE_INTERPOLATOR::ISOPARAMETRIC)

MakePair("WEIGHTED_AVERAGE", INTERFACE_INTERPOLATOR::WEIGHTED_AVERAGE)

MakePair("RADIAL_BASIS_FUNCTION", INTERFACE_INTERPOLATOR::RADIAL_BASIS_FUNCTION)

};

enum class RADIAL_BASIS {

};

static const MapType<std::string, RADIAL_BASIS> RadialBasisFunction_Map = {

MakePair("WENDLAND_C2", RADIAL_BASIS::WENDLAND_C2)

MakePair("INV_MULTI_QUADRIC", RADIAL_BASIS::INV_MULTI_QUADRIC)

MakePair("GAUSSIAN", RADIAL_BASIS::GAUSSIAN)

MakePair("THIN_PLATE_SPLINE", RADIAL_BASIS::THIN_PLATE_SPLINE)

MakePair("MULTI_QUADRIC", RADIAL_BASIS::MULTI_QUADRIC)

};

enum class INLET_SPANWISE_INTERP {

};

static const MapType<std::string, INLET_SPANWISE_INTERP> Inlet_SpanwiseInterpolation_Map = {

MakePair("NONE", INLET_SPANWISE_INTERP::NONE)

MakePair("LINEAR_1D", INLET_SPANWISE_INTERP::LINEAR_1D)

MakePair("AKIMA_1D", INLET_SPANWISE_INTERP::AKIMA_1D)

MakePair("CUBIC_1D", INLET_SPANWISE_INTERP::CUBIC_1D)

};

enum class INLET_INTERP_TYPE {

VR_VTHETA,

ALPHA_PHI,

};

static const MapType<std::string, INLET_INTERP_TYPE> Inlet_SpanwiseInterpolationType_Map = {

MakePair("VR_VTHETA", INLET_INTERP_TYPE::VR_VTHETA)

MakePair("ALPHA_PHI", INLET_INTERP_TYPE::ALPHA_PHI)

};

enum ENUM_TRANSFER {

};

enum class ENUM_REGIME {

};

enum ENUM_KIND_NONDIM {

};

static const MapType<std::string, ENUM_KIND_NONDIM> NonDim_Map = {

MakePair("DIMENSIONAL", DIMENSIONAL)

MakePair("FREESTREAM_PRESS_EQ_ONE", FREESTREAM_PRESS_EQ_ONE)

MakePair("FREESTREAM_VEL_EQ_MACH", FREESTREAM_VEL_EQ_MACH)

MakePair("FREESTREAM_VEL_EQ_ONE", FREESTREAM_VEL_EQ_ONE)

MakePair("INITIAL_VALUES", INITIAL_VALUES)

MakePair("REFERENCE_VALUES", REFERENCE_VALUES)

};

enum ENUM_MEASUREMENTS {

SI = 0, /*!< \brief Definition of compressible solver. */

US = 1 /*!< \brief Definition of incompressible solver. */

};

static const MapType<std::string, ENUM_MEASUREMENTS> Measurements_Map = {

MakePair("SI", SI)

MakePair("US", US)

};

enum RUNTIME_TYPE {

};

enum SOLVER_TYPE : const int {

FLOW_SOL=0, /*!< \brief Position of the mean flow solution in the solver container array. */

ADJFLOW_SOL=1, /*!< \brief Position of the continuous adjoint flow solution in the solver container array. */

TURB_SOL=2, /*!< \brief Position of the turbulence model solution in the solver container array. */

ADJTURB_SOL=3, /*!< \brief Position of the continuous adjoint turbulence solution in the solver container array. */

TRANS_SOL=4, /*!< \brief Position of the transition model solution in the solver container array. */

HEAT_SOL=5, /*!< \brief Position of the heat equation in the solution solver array. */

ADJHEAT_SOL=6, /*!< \brief Position of the adjoint heat equation in the solution solver array. */

RAD_SOL=7, /*!< \brief Position of the radiation equation in the solution solver array. */

ADJRAD_SOL=8, /*!< \brief Position of the continuous adjoint turbulence solution in the solver container array. */

MESH_SOL=9, /*!< \brief Position of the mesh solver. */

ADJMESH_SOL=10, /*!< \brief Position of the adjoint of the mesh solver. */

SPECIES_SOL=11, /*!< \brief Position of the species solver. */

ADJSPECIES_SOL=12, /*!< \brief Position of the adjoint of the species solver. */

FEA_SOL=0, /*!< \brief Position of the Finite Element flow solution in the solver container array. */

ADJFEA_SOL=1, /*!< \brief Position of the continuous adjoint Finite Element flow solution in the solver container array. */

TEMPLATE_SOL=0, /*!< \brief Position of the template solution. */

};

const int CONV_TERM = 0; /*!< \brief Position of the convective terms in the numerics container array. */

const int VISC_TERM = 1; /*!< \brief Position of the viscous terms in the numerics container array. */

const int SOURCE_FIRST_TERM = 2; /*!< \brief Position of the first source term in the numerics container array. */

const int SOURCE_SECOND_TERM = 3; /*!< \brief Position of the second source term in the numerics container array. */

const int CONV_BOUND_TERM = 4; /*!< \brief Position of the convective boundary terms in the numerics container array. */

const int VISC_BOUND_TERM = 5; /*!< \brief Position of the viscous boundary terms in the numerics container array. */

const int GRAD_TERM = 6; /*!< \brief Position of the gradient smoothing terms in the numerics container array. */

const int FEA_TERM = 0; /*!< \brief Position of the finite element analysis terms in the numerics container array. */

const int DE_TERM = 1; /*!< \brief Position of the dielectric terms in the numerics container array. */

const int MAT_NHCOMP = 2; /*!< \brief Position of the Neo-Hookean compressible material model. */

const int MAT_IDEALDE = 3; /*!< \brief Position of the Ideal-DE material model. */

const int MAT_KNOWLES = 4; /*!< \brief Position of the Knowles material model. */

const int EL_LINE = 6; /*!< \brief Elements of two nodes, with second order gauss quadrature (1D). */

const int EL_TRIA = 0; /*!< \brief Elements of three nodes (2D). */

const int EL_QUAD = 1; /*!< \brief Elements of four nodes (2D). */

const int EL_TRIA2 = 2; /*!< \brief Elements of three nodes (2D), with second order gauss quadrature. */

const int EL_TETRA = 0; /*!< \brief Elements of four nodes (3D). */

const int EL_HEXA = 1; /*!< \brief Elements of eight nodes (3D). */

const int EL_PYRAM = 2; /*!< \brief Elements of five nodes (3D). */

const int EL_PRISM = 3; /*!< \brief Elements of six nodes (3D). */

const int EL_TETRA2 = 4; /*!< \brief Elements of four nodes, with second order gauss quadrature (3D). */

const int EL_PYRAM2 = 5; /*!< \brief Elements of five nodes, with third order gauss quadrature (3D). */

enum ENUM_SPACE {

};

enum ENUM_FLUIDMODEL {

};

static const MapType<std::string, ENUM_FLUIDMODEL> FluidModel_Map = {

MakePair("STANDARD_AIR", STANDARD_AIR)

MakePair("IDEAL_GAS", IDEAL_GAS)

MakePair("VW_GAS", VW_GAS)

MakePair("PR_GAS", PR_GAS)

MakePair("CONSTANT_DENSITY", CONSTANT_DENSITY)

MakePair("INC_IDEAL_GAS", INC_IDEAL_GAS)

MakePair("INC_IDEAL_GAS_POLY", INC_IDEAL_GAS_POLY)

MakePair("MUTATIONPP", MUTATIONPP)

MakePair("SU2_NONEQ", SU2_NONEQ)

MakePair("FLUID_MIXTURE", FLUID_MIXTURE)

MakePair("COOLPROP", COOLPROP)

MakePair("DATADRIVEN_FLUID", DATADRIVEN_FLUID)

MakePair("FLUID_FLAMELET", FLUID_FLAMELET)

};

enum ENUM_GASMODEL {

};

static const MapType<std::string, ENUM_GASMODEL> GasModel_Map = {

MakePair("NONE", NO_MODEL)

MakePair("ARGON", ARGON)

MakePair("AIR-7", AIR7)

MakePair("AIR-21", AIR21)

MakePair("O2", O2)

MakePair("N2", N2)

MakePair("AIR-5", AIR5)

MakePair("ARGON-SID",ARGON_SID)

MakePair("ONESPECIES", ONESPECIES)

};

enum class ENUM_DATADRIVEN_METHOD {

LUT = 0,

MLP = 1

};

static const MapType<std::string, ENUM_DATADRIVEN_METHOD> DataDrivenMethod_Map = {

MakePair("LUT", ENUM_DATADRIVEN_METHOD::LUT)

MakePair("MLP", ENUM_DATADRIVEN_METHOD::MLP)

};

enum class TRANSCOEFFMODEL {

};

static const MapType<std::string, TRANSCOEFFMODEL> TransCoeffModel_Map = {

MakePair("SUTHERLAND", TRANSCOEFFMODEL::SUTHERLAND)

MakePair("WILKE", TRANSCOEFFMODEL::WILKE)

MakePair("GUPTA-YOS", TRANSCOEFFMODEL::GUPTAYOS)

MakePair("CHAPMANN-ENSKOG", TRANSCOEFFMODEL::CHAPMANN_ENSKOG)

};

enum class INC_DENSITYMODEL {

};

static const MapType<std::string, INC_DENSITYMODEL> DensityModel_Map = {

MakePair("CONSTANT", INC_DENSITYMODEL::CONSTANT)

MakePair("BOUSSINESQ", INC_DENSITYMODEL::BOUSSINESQ)

MakePair("VARIABLE", INC_DENSITYMODEL::VARIABLE)

MakePair("FLAMELET", INC_DENSITYMODEL::FLAMELET)

};

enum ENUM_INIT_OPTION {

REYNOLDS = 0, /*!< \brief Reynold's number initalization. */

TD_CONDITIONS = 1 /*!< \brief Total conditions initalization. */

};

static const MapType<std::string, ENUM_INIT_OPTION> InitOption_Map = {

MakePair("REYNOLDS", REYNOLDS)

MakePair("TD_CONDITIONS", TD_CONDITIONS)

};

enum class FREESTREAM_OPTION {

TEMPERATURE_FS, /*!< \brief Temperature initialization. */

DENSITY_FS, /*!< \brief Density initalization. */

};

static const MapType<std::string, FREESTREAM_OPTION> FreeStreamOption_Map = {

MakePair("TEMPERATURE_FS", FREESTREAM_OPTION::TEMPERATURE_FS)

MakePair("DENSITY_FS", FREESTREAM_OPTION::DENSITY_FS)

};

enum class VISCOSITYMODEL {

};

static const MapType<std::string, VISCOSITYMODEL> ViscosityModel_Map = {

MakePair("CONSTANT_VISCOSITY", VISCOSITYMODEL::CONSTANT)

MakePair("SUTHERLAND", VISCOSITYMODEL::SUTHERLAND)

MakePair("POLYNOMIAL_VISCOSITY", VISCOSITYMODEL::POLYNOMIAL)

MakePair("FLAMELET", VISCOSITYMODEL::FLAMELET)

MakePair("COOLPROP", VISCOSITYMODEL::COOLPROP)

};

enum class MIXINGVISCOSITYMODEL {

WILKE, /*!< \brief Wilke mixing viscosity model. */

DAVIDSON, /*!< \brief Davidson mixing viscosity model. */

};

static const MapType<std::string, MIXINGVISCOSITYMODEL> MixingViscosityModel_Map = {

MakePair("WILKE", MIXINGVISCOSITYMODEL::WILKE)

MakePair("DAVIDSON", MIXINGVISCOSITYMODEL::DAVIDSON)

};

enum class CONDUCTIVITYMODEL {

};

static const MapType<std::string, CONDUCTIVITYMODEL> ConductivityModel_Map = {

MakePair("CONSTANT_CONDUCTIVITY", CONDUCTIVITYMODEL::CONSTANT)

MakePair("CONSTANT_PRANDTL", CONDUCTIVITYMODEL::CONSTANT_PRANDTL)

MakePair("POLYNOMIAL_CONDUCTIVITY", CONDUCTIVITYMODEL::POLYNOMIAL)

MakePair("FLAMELET", CONDUCTIVITYMODEL::FLAMELET)

MakePair("COOLPROP", CONDUCTIVITYMODEL::COOLPROP)

};

enum class CONDUCTIVITYMODEL_TURB {

NONE, /*!< \brief No turbulent contribution to the effective thermal conductivity for RANS. */

CONSTANT_PRANDTL, /*!< \brief Include contribution to effective conductivity using constant turbulent Prandtl number for RANS. */

};

static const MapType<std::string, CONDUCTIVITYMODEL_TURB> TurbConductivityModel_Map = {

MakePair("NONE", CONDUCTIVITYMODEL_TURB::NONE)

MakePair("CONSTANT_PRANDTL_TURB", CONDUCTIVITYMODEL_TURB::CONSTANT_PRANDTL)

};

enum class DIFFUSIVITYMODEL {

};

static const MapType<std::string, DIFFUSIVITYMODEL> Diffusivity_Model_Map = {

MakePair("CONSTANT_DIFFUSIVITY", DIFFUSIVITYMODEL::CONSTANT_DIFFUSIVITY)

MakePair("CONSTANT_SCHMIDT", DIFFUSIVITYMODEL::CONSTANT_SCHMIDT)

MakePair("UNITY_LEWIS", DIFFUSIVITYMODEL::UNITY_LEWIS)

MakePair("CONSTANT_LEWIS", DIFFUSIVITYMODEL::CONSTANT_LEWIS)

MakePair("FLAMELET", DIFFUSIVITYMODEL::FLAMELET)

};

enum ENUM_GRIDMOVEMENT {

};

static const MapType<std::string, ENUM_GRIDMOVEMENT> GridMovement_Map = {

MakePair("NONE", NO_MOVEMENT)

MakePair("RIGID_MOTION", RIGID_MOTION)

MakePair("ROTATING_FRAME", ROTATING_FRAME)

MakePair("STEADY_TRANSLATION", STEADY_TRANSLATION)

MakePair("GUST", GUST)

};

enum ENUM_SURFACEMOVEMENT {

};

static const MapType<std::string, ENUM_SURFACEMOVEMENT> SurfaceMovement_Map = {

MakePair("DEFORMING", DEFORMING)

MakePair("MOVING_WALL", MOVING_WALL)

MakePair("AEROELASTIC_RIGID_MOTION", AEROELASTIC_RIGID_MOTION)

MakePair("AEROELASTIC", AEROELASTIC)

MakePair("EXTERNAL", EXTERNAL)

MakePair("EXTERNAL_ROTATION", EXTERNAL_ROTATION)

};

enum ENUM_GUST_TYPE {

};

static const MapType<std::string, ENUM_GUST_TYPE> Gust_Type_Map = {

MakePair("NONE", NO_GUST)

MakePair("TOP_HAT", TOP_HAT)

MakePair("SINE", SINE)

MakePair("ONE_M_COSINE", ONE_M_COSINE)

MakePair("VORTEX", VORTEX)

MakePair("EOG", EOG)

};

enum ENUM_GUST_DIR {

};

static const MapType<std::string, ENUM_GUST_DIR> Gust_Dir_Map = {

MakePair("X_DIR", X_DIR)

MakePair("Y_DIR", Y_DIR)

MakePair("Z_DIR", Z_DIR)

};

enum class CENTERED {

};

static const MapType<std::string, CENTERED> Centered_Map = {

MakePair("NONE", CENTERED::NONE)

MakePair("JST", CENTERED::JST)

MakePair("JST_KE", CENTERED::JST_KE)

MakePair("JST_MAT", CENTERED::JST_MAT)

MakePair("LAX-FRIEDRICH", CENTERED::LAX)

};

enum class UPWIND {

};

static const MapType<std::string, UPWIND> Upwind_Map = {

MakePair("NONE", UPWIND::NONE)

MakePair("ROE", UPWIND::ROE)

MakePair("TURKEL_PREC", UPWIND::TURKEL)

MakePair("AUSM", UPWIND::AUSM)

MakePair("AUSMPLUSUP", UPWIND::AUSMPLUSUP)

MakePair("AUSMPLUSUP2", UPWIND::AUSMPLUSUP2)

MakePair("AUSMPLUSM", UPWIND::AUSMPLUSM)

MakePair("SLAU", UPWIND::SLAU)

MakePair("HLLC", UPWIND::HLLC)

MakePair("SW", UPWIND::SW)

MakePair("MSW", UPWIND::MSW)

MakePair("SCALAR_UPWIND", UPWIND::SCALAR_UPWIND)

MakePair("BOUNDED_SCALAR", UPWIND::BOUNDED_SCALAR)

MakePair("CONVECTIVE_TEMPLATE", UPWIND::CONVECTIVE_TEMPLATE)

MakePair("L2ROE", UPWIND::L2ROE)

MakePair("LMROE", UPWIND::LMROE)

MakePair("SLAU2", UPWIND::SLAU2)

MakePair("FDS", UPWIND::FDS)

MakePair("LAX-FRIEDRICH", UPWIND::LAX_FRIEDRICH)

};

enum ENUM_FEM {

NO_FEM = 0, /*!< \brief No finite element scheme is used. */

DG = 1 /*!< \brief Discontinuous Galerkin numerical method. */

};

static const MapType<std::string, ENUM_FEM> FEM_Map = {

MakePair("NONE", NO_FEM)

MakePair("DG", DG)

};

enum class FEM_SHOCK_CAPTURING_DG {

NONE, /*!< \brief Shock capturing is not used. */

PERSSON /*!< \brief Per-Olof Persson's sub-cell shock capturing method. */

};

static const MapType<std::string, FEM_SHOCK_CAPTURING_DG> ShockCapturingDG_Map = {

MakePair("NONE", FEM_SHOCK_CAPTURING_DG::NONE)

MakePair("PERSSON", FEM_SHOCK_CAPTURING_DG::PERSSON)

};

enum ENUM_MATRIX_COLORING {

GREEDY_COLORING = 0, /*!< \brief Greedy type of algorithm for the coloring. */

NATURAL_COLORING = 1 /*!< \brief One color for every DOF, very slow. Only to be used for debugging. */

};

static const MapType<std::string, ENUM_MATRIX_COLORING> MatrixColoring_Map = {

MakePair("GREEDY_COLORING", GREEDY_COLORING)

MakePair("NATURAL_COLORING", NATURAL_COLORING)

};

enum class LIMITER {

};

static const MapType<std::string, LIMITER> Limiter_Map = {

MakePair("NONE", LIMITER::NONE)

MakePair("VENKATAKRISHNAN", LIMITER::VENKATAKRISHNAN)

MakePair("NISHIKAWA_R3", LIMITER::NISHIKAWA_R3)

MakePair("NISHIKAWA_R4", LIMITER::NISHIKAWA_R4)

MakePair("NISHIKAWA_R5", LIMITER::NISHIKAWA_R5)

MakePair("VENKATAKRISHNAN_WANG", LIMITER::VENKATAKRISHNAN_WANG)

MakePair("BARTH_JESPERSEN", LIMITER::BARTH_JESPERSEN)

MakePair("VAN_ALBADA_EDGE", LIMITER::VAN_ALBADA_EDGE)

MakePair("SHARP_EDGES", LIMITER::SHARP_EDGES)

MakePair("WALL_DISTANCE", LIMITER::WALL_DISTANCE)

};

enum class TURB_MODEL {

};

static const MapType<std::string, TURB_MODEL> Turb_Model_Map = {

MakePair("NONE", TURB_MODEL::NONE)

MakePair("SA", TURB_MODEL::SA)

MakePair("SST", TURB_MODEL::SST)

};

enum class TURB_FAMILY {

};

inline TURB_FAMILY TurbModelFamily(TURB_MODEL model) {

}

enum class SST_OPTIONS {

};

static const MapType<std::string, SST_OPTIONS> SST_Options_Map = {

MakePair("NONE", SST_OPTIONS::NONE)