fe_data.cc

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

#include <deal.II/base/geometry_info.h>

#include <deal.II/fe/fe.h>

DEAL_II_NAMESPACE_OPEN

template <int dim>
FiniteElementData<dim>::FiniteElementData(
  const std::vector<unsigned int> &dofs_per_object,
  const unsigned int               n_components,
  const unsigned int               degree,
  const Conformity                 conformity,
  const BlockIndices &             block_indices)
  : dofs_per_vertex(dofs_per_object[0])
  , dofs_per_line(dofs_per_object[1])
  , dofs_per_quad(dim > 1 ? dofs_per_object[2] : 0)
  , dofs_per_hex(dim > 2 ? dofs_per_object[3] : 0)
  , first_line_index(GeometryInfo<dim>::vertices_per_cell * dofs_per_vertex)
  , first_quad_index(first_line_index +
                     GeometryInfo<dim>::lines_per_cell * dofs_per_line)
  , first_hex_index(first_quad_index +
                    GeometryInfo<dim>::quads_per_cell * dofs_per_quad)
  , first_face_line_index(GeometryInfo<dim - 1>::vertices_per_cell *
                          dofs_per_vertex)
  , first_face_quad_index(
      (dim == 3 ? GeometryInfo<dim - 1>::vertices_per_cell * dofs_per_vertex :
                  GeometryInfo<dim>::vertices_per_cell * dofs_per_vertex) +
      GeometryInfo<dim - 1>::lines_per_cell * dofs_per_line)
  , dofs_per_face(GeometryInfo<dim>::vertices_per_face * dofs_per_vertex +
                  GeometryInfo<dim>::lines_per_face * dofs_per_line +
                  GeometryInfo<dim>::quads_per_face * dofs_per_quad)
  , dofs_per_cell(GeometryInfo<dim>::vertices_per_cell * dofs_per_vertex +
                  GeometryInfo<dim>::lines_per_cell * dofs_per_line +
                  GeometryInfo<dim>::quads_per_cell * dofs_per_quad +
                  GeometryInfo<dim>::hexes_per_cell * dofs_per_hex)
  , components(n_components)
  , degree(degree)
  , conforming_space(conformity)
  , block_indices_data(block_indices.size() == 0 ?
                         BlockIndices(1, dofs_per_cell) :
                         block_indices)
{
  Assert(dofs_per_object.size() == dim + 1,
         ExcDimensionMismatch(dofs_per_object.size() - 1, dim));
}



template <int dim>
bool
FiniteElementData<dim>::operator==(const FiniteElementData<dim> &f) const
{
  return ((dofs_per_vertex == f.dofs_per_vertex) &&
          (dofs_per_line == f.dofs_per_line) &&
          (dofs_per_quad == f.dofs_per_quad) &&
          (dofs_per_hex == f.dofs_per_hex) && (components == f.components) &&
          (degree == f.degree) && (conforming_space == f.conforming_space));
}


template class FiniteElementData<1>;
template class FiniteElementData<2>;
template class FiniteElementData<3>;

DEAL_II_NAMESPACE_CLOSE