fe_values_views_internal.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_fe_values_views_internal_h
#define dealii_fe_values_views_internal_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/array_view.h>
#include <deal.II/base/symmetric_tensor.h>
#include <deal.II/base/table.h>
#include <deal.II/base/tensor.h>
 
#include <deal.II/fe/fe_values_views.h>
 
#include <vector>
 
DEAL_II_NAMESPACE_OPEN
 
namespace FEValuesViews
{
  namespace internal
  {
    // Given values of degrees of freedom, evaluate the
    // values/gradients/... at quadrature points
 
    // ------------------------- scalar functions --------------------------
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_values(
      const ArrayView<const Number> &dof_values,
      const Table<2, double>        &shape_values,
      const std::vector<typename Scalar<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename ProductType<Number, double>::type> &values);
 
    template <int order, int dim, int spacedim, typename Number>
    void
    do_function_derivatives(
      const ArrayView<const Number>                   &dof_values,
      const Table<2, ::Tensor<order, spacedim>> &shape_derivatives,
      const std::vector<typename Scalar<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number, ::Tensor<order, spacedim>>::type>
        &derivatives);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_laplacians(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<2, spacedim>> &shape_hessians,
      const std::vector<typename Scalar<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename Scalar<dim, spacedim>::
                    template solution_laplacian_type<Number>> &laplacians);
 
    // ----------------------------- vector part ---------------------------
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_values(
      const ArrayView<const Number> &dof_values,
      const Table<2, double>        &shape_values,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number, ::Tensor<1, spacedim>>::type>
        &values);
 
    template <int order, int dim, int spacedim, typename Number>
    void
    do_function_derivatives(
      const ArrayView<const Number>                   &dof_values,
      const Table<2, ::Tensor<order, spacedim>> &shape_derivatives,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number, ::Tensor<order + 1, spacedim>>::type>
        &derivatives);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_symmetric_gradients(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number,
                             ::SymmetricTensor<2, spacedim>>::type>
        &symmetric_gradients);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_divergences(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename Vector<dim, spacedim>::
                    template solution_divergence_type<Number>> &divergences);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_curls(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename ProductType<
        Number,
        typename ::internal::CurlType<spacedim>::type>::type> &curls);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_laplacians(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<2, spacedim>> &shape_hessians,
      const std::vector<typename Vector<dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename Vector<dim, spacedim>::
                    template solution_laplacian_type<Number>> &laplacians);
 
    // ---------------------- symmetric tensor part ------------------------
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_values(
      const ArrayView<const Number>  &dof_values,
      const ::Table<2, double> &shape_values,
      const std::vector<
        typename SymmetricTensor<2, dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number,
                             ::SymmetricTensor<2, spacedim>>::type>
        &values);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_divergences(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<
        typename SymmetricTensor<2, dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename SymmetricTensor<2, dim, spacedim>::
                    template solution_divergence_type<Number>> &divergences);
 
    // ---------------------- non-symmetric tensor part ------------------------
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_values(
      const ArrayView<const Number>  &dof_values,
      const ::Table<2, double> &shape_values,
      const std::vector<typename Tensor<2, dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<
        typename ProductType<Number, ::Tensor<2, spacedim>>::type>
        &values);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_divergences(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<typename Tensor<2, dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename Tensor<2, dim, spacedim>::
                    template solution_divergence_type<Number>> &divergences);
 
    template <int dim, int spacedim, typename Number>
    void
    do_function_gradients(
      const ArrayView<const Number>               &dof_values,
      const Table<2, ::Tensor<1, spacedim>> &shape_gradients,
      const std::vector<typename Tensor<2, dim, spacedim>::ShapeFunctionData>
        &shape_function_data,
      std::vector<typename Tensor<2, dim, spacedim>::
                    template solution_gradient_type<Number>> &gradients);
  } // end of namespace internal
} // namespace FEValuesViews
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif