mapping_c1.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2001 - 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 at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#ifndef dealii_mapping_c1_h
#define dealii_mapping_c1_h


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

DEAL_II_NAMESPACE_OPEN


template <int dim, int spacedim=dim>
class MappingC1 : public MappingQ<dim,spacedim>
{
public:
  MappingC1 ();

  virtual
  std::unique_ptr<Mapping<dim,spacedim>> clone () const override;

protected:

  class MappingC1Generic : public MappingQGeneric<dim,spacedim>
  {
  public:

    MappingC1Generic ();

    virtual
    void
    add_line_support_points (const typename Triangulation<dim>::cell_iterator &cell,
                             std::vector<Point<dim> > &a) const override;

    virtual
    void
    add_quad_support_points(const typename Triangulation<dim>::cell_iterator &cell,
                            std::vector<Point<dim> > &a) const override;
  };
};

/* -------------- declaration of explicit specializations ------------- */

#ifndef DOXYGEN

template <> void MappingC1<1>::MappingC1Generic::add_line_support_points (
  const Triangulation<1>::cell_iterator &,
  std::vector<Point<1> > &) const;
template <> void MappingC1<2>::MappingC1Generic::add_line_support_points (
  const Triangulation<2>::cell_iterator &cell,
  std::vector<Point<2> > &a) const;

template <> void MappingC1<1>::MappingC1Generic::add_quad_support_points (
  const Triangulation<1>::cell_iterator &,
  std::vector<Point<1> > &) const;
template <> void MappingC1<2>::MappingC1Generic::add_quad_support_points (
  const Triangulation<2>::cell_iterator &,
  std::vector<Point<2> > &) const;


#endif // DOXYGEN

DEAL_II_NAMESPACE_CLOSE

#endif