fe_p1nc.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2016 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_fe_p1nc_h
#define dealii_fe_p1nc_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/qprojector.h>
#include <deal.II/base/quadrature.h>
 
#include <deal.II/fe/fe.h>
 
 
DEAL_II_NAMESPACE_OPEN
 
 
class FE_P1NC : public FiniteElement<2, 2>
{
public:
  FE_P1NC();
 
  virtual std::string
  get_name() const override;
 
  virtual UpdateFlags
  requires_update_flags(const UpdateFlags flags) const override;
 
  virtual std::unique_ptr<FiniteElement<2, 2>>
  clone() const override;
 
  virtual ~FE_P1NC() override = default;
 
 
 
private:
  static std::vector<unsigned int>
  get_dpo_vector();
 
  static ndarray<double, 4, 3>
  get_linear_shape_coefficients(const Triangulation<2, 2>::cell_iterator &cell);
 
  virtual std::unique_ptr<FiniteElement<2, 2>::InternalDataBase>
  get_data(
    const UpdateFlags update_flags,
    const Mapping<2, 2> &,
    const Quadrature<2> &quadrature,
    ::internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  using FiniteElement<2, 2>::get_face_data;
 
  virtual std::unique_ptr<FiniteElement<2, 2>::InternalDataBase>
  get_face_data(
    const UpdateFlags update_flags,
    const Mapping<2, 2> &,
    const hp::QCollection<1> &quadrature,
    ::internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  virtual std::unique_ptr<FiniteElement<2, 2>::InternalDataBase>
  get_subface_data(
    const UpdateFlags update_flags,
    const Mapping<2, 2> &,
    const Quadrature<1> &quadrature,
    ::internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  virtual void
  fill_fe_values(
    const Triangulation<2, 2>::cell_iterator &cell,
    const CellSimilarity::Similarity          cell_similarity,
    const Quadrature<2>                      &quadrature,
    const Mapping<2, 2>                      &mapping,
    const Mapping<2, 2>::InternalDataBase    &mapping_internal,
    const internal::FEValuesImplementation::MappingRelatedData<2, 2>
                                                &mapping_data,
    const FiniteElement<2, 2>::InternalDataBase &fe_internal,
    internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  using FiniteElement<2, 2>::fill_fe_face_values;
 
  virtual void
  fill_fe_face_values(
    const Triangulation<2, 2>::cell_iterator &cell,
    const unsigned int                        face_no,
    const hp::QCollection<1>                 &quadrature,
    const Mapping<2, 2>                      &mapping,
    const Mapping<2, 2>::InternalDataBase    &mapping_internal,
    const internal::FEValuesImplementation::MappingRelatedData<2, 2>
                           &mapping_data,
    const InternalDataBase &fe_internal,
    ::internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  virtual void
  fill_fe_subface_values(
    const Triangulation<2, 2>::cell_iterator &cell,
    const unsigned int                        face_no,
    const unsigned int                        sub_no,
    const Quadrature<1>                      &quadrature,
    const Mapping<2, 2>                      &mapping,
    const Mapping<2, 2>::InternalDataBase    &mapping_internal,
    const internal::FEValuesImplementation::MappingRelatedData<2, 2>
                           &mapping_data,
    const InternalDataBase &fe_internal,
    ::internal::FEValuesImplementation::FiniteElementRelatedData<2, 2>
      &output_data) const override;
 
  void
  initialize_constraints();
};
 
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif