fe_p1nc.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2015 - 2016 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__fe_p1nc_h
#define dealii__fe_p1nc_h

#include <deal.II/base/config.h>
#include <deal.II/fe/fe.h>
#include <deal.II/base/quadrature.h>
#include <deal.II/base/qprojector.h>


DEAL_II_NAMESPACE_OPEN



class FE_P1NC : public FiniteElement<2,2>
{

public:
  FE_P1NC();

  virtual std::string get_name () const;

  virtual UpdateFlags requires_update_flags (const UpdateFlags flags) const;

  virtual FiniteElement<2,2> *clone () const;

  virtual ~FE_P1NC ();



private:

  static std::vector<unsigned int> get_dpo_vector ();

  static std_cxx11::array<std_cxx11::array<double,3>,4>
  get_linear_shape_coefficients (const Triangulation<2,2>::cell_iterator &cell);

  virtual FiniteElement<2,2>::InternalDataBase *
  get_data (const UpdateFlags update_flags,
            const Mapping<2,2> &,
            const Quadrature<2> &quadrature,
            ::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  virtual FiniteElement<2,2>::InternalDataBase *
  get_face_data (const UpdateFlags update_flags,
                 const Mapping<2,2> &,
                 const Quadrature<1> &quadrature,
                 ::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  virtual FiniteElement<2,2>::InternalDataBase *
  get_subface_data (const UpdateFlags update_flags,
                    const Mapping<2,2> &,
                    const Quadrature<1> &quadrature,
                    ::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  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::FEValues::MappingRelatedData<2,2> &mapping_data,
                  const FiniteElement<2,2>::InternalDataBase        &fe_internal,
                  internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  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::FEValues::MappingRelatedData<2,2> &mapping_data,
                       const InternalDataBase                                    &fe_internal,
                       ::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  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::FEValues::MappingRelatedData<2,2> &mapping_data,
                          const InternalDataBase                                    &fe_internal,
                          ::internal::FEValues::FiniteElementRelatedData<2,2> &output_data) const;

  void initialize_constraints ();

};




DEAL_II_NAMESPACE_CLOSE

#endif