curl_compute.h

/*  ______________________________________________________________________
 *
 *  ExaDG - High-Order Discontinuous Galerkin for the Exa-Scale
 *
 *  Copyright (C) 2021 by the ExaDG authors
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  This program 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 General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 *  ______________________________________________________________________
 */
 
#ifndef INCLUDE_EXADG_INCOMPRESSIBLE_NAVIER_STOKES_SPATIAL_DISCRETIZATION_CURL_COMPUTE_H_
#define INCLUDE_EXADG_INCOMPRESSIBLE_NAVIER_STOKES_SPATIAL_DISCRETIZATION_CURL_COMPUTE_H_
 
#include <deal.II/base/tensor.h>
#include <deal.II/base/vectorization.h>
 
namespace ExaDG
{
namespace IncNS
{
template<int dim, typename FEEval>
struct CurlCompute
{
  static dealii::Tensor<1, dim, dealii::VectorizedArray<typename FEEval::number_type>>
  compute(FEEval & fe_eval, unsigned int const q_point)
  {
    return fe_eval.get_curl(q_point);
  }
};
 
// use partial specialization of templates to handle the 2-dimensional case
template<typename FEEval>
struct CurlCompute<2, FEEval>
{
  static dealii::Tensor<1, 2, dealii::VectorizedArray<typename FEEval::number_type>>
  compute(FEEval & fe_eval, unsigned int const q_point)
  {
    dealii::Tensor<1, 2, dealii::VectorizedArray<typename FEEval::number_type>> curl;
    /*
     * fe_eval = / phi \   _____\   grad(fe_eval) = / d(phi)/dx1   d(phi)/dx2 \
     *           \  0  /        /                   \     0             0     /
     */
    dealii::Tensor<2, 2, dealii::VectorizedArray<typename FEEval::number_type>> temp =
      fe_eval.get_gradient(q_point);
    /*
     *         __    /  0  \     /   d(phi)/dx2 \
     *  curl = \/ X  |  0  |  =  | - d(phi)/dx1 |
     *               \ phi /     \      0       /
     */
    curl[0] = temp[0][1];  //   d(phi)/dx2
    curl[1] = -temp[0][0]; // - d(phi)/dx1
    return curl;
  }
};
 
} // namespace IncNS
} // namespace ExaDG
 
 
#endif /* INCLUDE_EXADG_INCOMPRESSIBLE_NAVIER_STOKES_SPATIAL_DISCRETIZATION_CURL_COMPUTE_H_ */