variables.h

Go to the documentation of this file.

// Copyright (C) 2024 The m-AIA AUTHORS
//
// This file is part of m-AIA (https://git.rwth-aachen.de/aia/m-AIA/m-AIA)
//
// SPDX-License-Identifier: LGPL-3.0-only
 
#ifndef MAIA_VARIABLES_H_
#define MAIA_VARIABLES_H_
#include <array>
#include <limits>
#include "INCLUDE/maiatypes.h"
 
template <MInt nDim>
class MConservativeVariables {
 public:
  explicit MConservativeVariables(MInt noSpecies, MInt noRANSEq = 0) {
    m_noSpecies = noSpecies;
    m_noRansEquations = noRANSEq;
    for(MInt i = 0; i < nDim; i++) {
      RHO_VV[i] = RHO_U + i;
    }
    RHO_Y = new MInt[m_noSpecies];
    for(MInt i = 0; i < m_noSpecies; i++) {
      RHO_Y[i] = RHO_Z + i;
    }
    RANS_FIRST = nDim + 2;
    if(m_noSpecies != 0) {
      RANS_FIRST = nDim + 2 + m_noSpecies;
    }
    if(m_noRansEquations > 0) {
      RANS_VAR = new MInt[m_noRansEquations];
      for(MInt i = 0; i < m_noRansEquations; ++i) {
        RANS_VAR[i] = RANS_FIRST + i;
      }
    } else {
      RANS_VAR = nullptr;
    }
    noVariables = nDim + 2 + m_noSpecies + m_noRansEquations;
  }
 
  ~MConservativeVariables() {
    delete[] RHO_Y;
    if(RANS_VAR != nullptr) {
      delete[] RANS_VAR;
    }
  }
 
  static constexpr MInt RHO_U = 0;
  static constexpr MInt RHO_V = 1;
  static constexpr MInt RHO_W = 2;
  std::array<MInt, nDim> RHO_VV;
  static constexpr MInt RHO_E = nDim;
  static constexpr MInt RHO = nDim + 1;
  static constexpr MInt RHO_Z = nDim + 2;
  static constexpr MInt RHO_C = nDim + 2;
  MInt* RHO_Y = nullptr;
  MInt RANS_FIRST;
  MInt* RANS_VAR = nullptr;
 
 
  MInt noVariables;
  MInt m_noSpecies;
  MInt m_noRansEquations;
  MFloat rhoUInfinity{}, rhoVInfinity{}, rhoWInfinity{}, rhoEInfinity{}, rhoInfinity{};
  std::array<MFloat, nDim> ransInfinity{};
  std::array<MFloat, nDim> rhoVVInfinity{};
};
 
template <MInt nDim>
class MPrimitiveVariables {
 public:
  explicit MPrimitiveVariables(MInt noSpecies, MInt noRANSEq = 0) {
    m_noSpecies = noSpecies;
    m_noRansEquations = noRANSEq;
 
    Y = new MInt[m_noSpecies];
    for(MInt i = 0; i < m_noSpecies; i++) {
      Y[i] = Z + i;
    }
    // G = nDim+1;
 
    RANS_FIRST = nDim + 2;
    if(m_noSpecies != 0) {
      RANS_FIRST = nDim + 2 + m_noSpecies;
    }
    if(m_noRansEquations > 0) {
      RANS_VAR = new MInt[m_noRansEquations];
      for(MInt i = 0; i < m_noRansEquations; ++i) {
        RANS_VAR[i] = RANS_FIRST + i;
      }
    } else {
      RANS_VAR = nullptr;
    }
 
    noVariables = nDim + 2 + m_noSpecies + m_noRansEquations;
  }
  ~MPrimitiveVariables() {
    delete[] Y;
    if(RANS_VAR != nullptr) {
      delete[] RANS_VAR;
    }
  }
 
  static constexpr MInt U = 0;
  static constexpr MInt V = 1;
  static constexpr MInt W = 2;
#if !defined(MAIA_PGI_COMPILER) && defined(NDEBUG) && !defined(MAIA_SANITIZE_ADDRESS) && !defined(MAIA_INTEL_COMPILER)
  static constexpr std::array<MInt, 3> VV = {0, 1, 2};
#else
  std::array<MInt, 3> VV = {0, 1, 2};
#endif
  static constexpr MInt P = nDim + 1;
  static constexpr MInt RHO = nDim;
  static constexpr MInt T = nDim + 1;
  static constexpr MInt Z = nDim + 2;
  static constexpr MInt C = nDim + 2;
  MInt* Y;
  MInt RANS_FIRST;
  MInt* RANS_VAR;
 
  MInt noVariables;
  MInt m_noSpecies;
  MInt m_noRansEquations;
  MFloat UInfinity{}, VInfinity{}, WInfinity{}, PInfinity{}, TInfinity{};
  std::array<MFloat, nDim> VVInfinity;
  std::array<MFloat, nDim> ransInfinity;
  MFloat DthInfinity{}, muInfinity{}, DInfinity{};
};
 
 
namespace maia {
 
namespace fv {
 
namespace variables {
static const MInt Segfault = std::numeric_limits<MInt>::min();
} // namespace variables
 
template <MInt nd>
struct ConservativeVariables {
  static constexpr MInt RHO_U = 0;
  static constexpr MInt RHO_V = 1;
  static constexpr MInt RHO_W = nd == 3 ? 2 : variables::Segfault;
  static constexpr std::array<MInt, 3> RHO_VV = {0, 1, 2};
  static constexpr MInt RHO_E = nd;
  static constexpr MInt RHO = nd + 1;
  const MInt RHO_Z;
  const MInt RHO_C;
  static constexpr MInt RHO_N = nd + 2;
 
  const MInt m_noSpecies;
  const MInt m_noRansEquations;
  const MInt noVariables;
 
  MInt* RHO_NN = nullptr;
  MInt* RHO_Y = nullptr;
 
  explicit ConservativeVariables(const MInt noSpecies, const MInt noRans);
  ~ConservativeVariables();
};
 
template <MInt nd>
struct PrimitiveVariables {
  static constexpr MInt U = 0;
  static constexpr MInt V = 1;
  static constexpr MInt W = nd == 3 ? 2 : variables::Segfault;
  static constexpr std::array<MInt, 3> VV = {0, 1, 2};
  static constexpr MInt RHO = nd;
  static constexpr MInt P = nd + 1;
  static constexpr MInt T = nd + 1;
 
  const MInt Z;
  const MInt C;
  static constexpr MInt N = nd + 2;
 
  const MInt m_noSpecies;
  const MInt m_noRansEquations;
  const MInt noVariables;
 
  MInt* NN = nullptr;
  MInt* Y = nullptr;
 
  explicit PrimitiveVariables(const MInt noSpecies, const MInt noRans);
  ~PrimitiveVariables();
};
 
} // namespace fv
} // namespace maia
#endif