shape_info.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2011 - 2017 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, 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.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------


#ifndef dealii__matrix_free_shape_info_h
#define dealii__matrix_free_shape_info_h


#include <deal.II/base/exceptions.h>
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/base/vectorization.h>
#include <deal.II/base/aligned_vector.h>
#include <deal.II/fe/fe.h>

#include <deal.II/matrix_free/helper_functions.h>


DEAL_II_NAMESPACE_OPEN

namespace internal
{
  namespace MatrixFreeFunctions
  {
    enum ElementType
    {
      tensor_general,
      tensor_symmetric,
      truncated_tensor,
      tensor_symmetric_plus_dg0,
      tensor_gausslobatto
    };

    template <typename Number>
    struct ShapeInfo
    {
      ShapeInfo ();

      template <int dim>
      ShapeInfo (const Quadrature<1> &quad,
                 const FiniteElement<dim> &fe,
                 const unsigned int base_element = 0);

      template <int dim>
      void reinit (const Quadrature<1> &quad,
                   const FiniteElement<dim> &fe_dim,
                   const unsigned int base_element = 0);

      std::size_t memory_consumption () const;

      ElementType element_type;

      AlignedVector<VectorizedArray<Number> > shape_values;

      AlignedVector<VectorizedArray<Number> > shape_gradients;

      AlignedVector<VectorizedArray<Number> > shape_hessians;

      AlignedVector<VectorizedArray<Number> > shape_val_evenodd;

      AlignedVector<VectorizedArray<Number> > shape_gra_evenodd;

      AlignedVector<VectorizedArray<Number> > shape_hes_evenodd;

      ::Table<2,unsigned int>  face_indices;

      std::vector<Number>    face_value[2];

      std::vector<Number>    face_gradient[2];

      std::vector<Number>    subface_value[2];

      std::vector<Number>    shape_values_number;

      std::vector<Number>    shape_gradient_number;

      std::vector<unsigned int> lexicographic_numbering;

      unsigned int fe_degree;

      unsigned int n_q_points_1d;

      unsigned int n_q_points;

      unsigned int dofs_per_cell;

      unsigned int n_q_points_face;

      unsigned int dofs_per_face;

      bool check_1d_shapes_symmetric(const unsigned int n_q_points_1d);

      bool check_1d_shapes_gausslobatto();
    };



    // ------------------------------------------ inline functions

    template <typename Number>
    template <int dim>
    inline
    ShapeInfo<Number>::ShapeInfo (const Quadrature<1> &quad,
                                  const FiniteElement<dim> &fe_in,
                                  const unsigned int base_element_number)
      :
      fe_degree (0),
      n_q_points_1d (0),
      n_q_points (0),
      dofs_per_cell (0),
      n_q_points_face (0),
      dofs_per_face (0)
    {
      reinit (quad, fe_in, base_element_number);
    }



  } // end of namespace MatrixFreeFunctions
} // end of namespace internal

DEAL_II_NAMESPACE_CLOSE

#endif