fe_p1nc.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2015 - 2018 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.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------

#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 std::array<std::array<double, 3>, 4>
  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;

  virtual std::unique_ptr<FiniteElement<2, 2>::InternalDataBase>
  get_face_data(
    const UpdateFlags update_flags,
    const Mapping<2, 2> &,
    const Quadrature<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;

  virtual void
  fill_fe_face_values(
    const Triangulation<2, 2>::cell_iterator &cell,
    const unsigned int                        face_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;

  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