CNumericsSIMD.hpp

Go to the documentation of this file.

/*!
 * \file CNumericsSIMD.hpp
 * \brief Vectorized (SIMD) numerics classes.
 * \author P. Gomes
 * \version 8.5.0 "Harrier"
 *
 * SU2 Project Website: https://su2code.github.io
 *
 * The SU2 Project is maintained by the SU2 Foundation
 * (http://su2foundation.org)
 *
 * Copyright 2012-2026, SU2 Contributors (cf. AUTHORS.md)
 *
 * SU2 is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * SU2 is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with SU2. If not, see <http://www.gnu.org/licenses/>.
 */
 
#pragma once
 
#include "../../../Common/include/parallelization/vectorization.hpp"
#include "../numerics/util.hpp"
 
/*!
 * \brief Define Double and Int SIMD types.
 */
using Double = simd::Array<su2double>;
using Int = simd::Array<unsigned long, Double::Size>;
 
/*--- Forward declare a few classes used in name only by the interface. ---*/
class CConfig;
class CGeometry;
class CVariable;
 
/*!
 * \class CNumericsSIMD
 * \ingroup ConvDiscr
 * \brief Base class to define the interface.
 * \note See CNumericsEmptyDecorator.
 */
class CNumericsSIMD {
public:
  /*!
   * \brief Interface for edge flux computation.
   * \param[in] iEdge - The edges for flux computation.
   * \param[in] config - Problem definitions.
   * \param[in] geometry - Problem geometry.
   * \param[in] solution - Solution variables.
   * \param[in] updateType - Type of update done on vector and matrix.
   * \param[in] updateMask - SIMD array of 1's and 0's, the latter prevent the update.
   * \param[in,out] vector - Target for the fluxes.
   * \param[in,out] matrix - Target for the flux Jacobians.
   * \param[out] edgeMassFluxes - Optional, per-edge mass flux (density-equation flux),
   *             used by "bounded scalar" discretization of transported scalars (turbulence, species).
   * \note The update mask is used to handle "remainder" edges (nEdge mod simdSize).
   */
  virtual void ComputeFlux(Int iEdge,
                           const CConfig& config,
                           const CGeometry& geometry,
                           const CVariable& solution,
                           UpdateType updateType,
                           Double updateMask,
                           CSysVector<su2double>& vector,
                           SparseMatrixType& matrix,
                           su2activevector* edgeMassFluxes) const = 0;
 
  /*! \brief Destructor of the class. */
  virtual ~CNumericsSIMD(void) = default;
 
  /*!
   * \brief Factory method.
   * \param[in] config - Problem definitions.
   * \param[in] nDim - 2D or 3D.
   * \param[in] iMesh - Grid index.
   * \param[in] turbVars - Turbulence variables.
   */
  static CNumericsSIMD* CreateNumerics(const CConfig& config, int nDim, int iMesh, const CVariable* turbVars = nullptr);
 
};