/usr/include/deal.II/fe/fe_dg_vector.templates.h is in libdeal.ii-dev 8.4.2-2+b1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 | // ---------------------------------------------------------------------
//
// Copyright (C) 2006 - 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.
//
// ---------------------------------------------------------------------
#include <deal.II/fe/fe_dg_vector.h>
#include <deal.II/fe/fe_tools.h>
#include <deal.II/base/quadrature_lib.h>
DEAL_II_NAMESPACE_OPEN
//TODO:[GK] deg+1 is wrong here and should be fixed after FiniteElementData was cleaned up
template <class PolynomialType, int dim, int spacedim>
FE_DGVector<PolynomialType,dim,spacedim>::FE_DGVector (
const unsigned int deg, MappingType map)
:
FE_PolyTensor<PolynomialType, dim, spacedim>(
deg,
FiniteElementData<dim>(
get_dpo_vector(deg), dim, deg+1, FiniteElementData<dim>::L2),
std::vector<bool>(PolynomialType::compute_n_pols(deg), true),
std::vector<ComponentMask>(PolynomialType::compute_n_pols(deg),
ComponentMask(dim,true)))
{
this->mapping_type = map;
const unsigned int polynomial_degree = this->tensor_degree();
QGauss<dim> quadrature(polynomial_degree+1);
this->generalized_support_points = quadrature.get_points();
this->reinit_restriction_and_prolongation_matrices(true, true);
FETools::compute_projection_matrices (*this, this->restriction, true);
FETools::compute_embedding_matrices (*this, this->prolongation, true);
}
template <class PolynomialType, int dim, int spacedim>
FiniteElement<dim, spacedim> *
FE_DGVector<PolynomialType,dim,spacedim>::clone() const
{
return new FE_DGVector<PolynomialType, dim, spacedim>(*this);
}
template <class PolynomialType, int dim, int spacedim>
std::string
FE_DGVector<PolynomialType,dim,spacedim>::get_name() const
{
std::ostringstream namebuf;
namebuf << "FE_DGVector_" << this->poly_space.name()
<< "<" << dim << ">(" << this->degree-1 << ")";
return namebuf.str();
}
template <class PolynomialType, int dim, int spacedim>
std::vector<unsigned int>
FE_DGVector<PolynomialType,dim,spacedim>::get_dpo_vector (const unsigned int deg)
{
std::vector<unsigned int> dpo(dim+1);
dpo[dim] = PolynomialType::compute_n_pols(deg);
return dpo;
}
template <class PolynomialType, int dim, int spacedim>
bool
FE_DGVector<PolynomialType,dim,spacedim>::has_support_on_face
(const unsigned int,
const unsigned int) const
{
return true;
}
template <class PolynomialType, int dim, int spacedim>
void
FE_DGVector<PolynomialType,dim,spacedim>::interpolate
(std::vector<double> &,
const std::vector<double> &) const
{
Assert(false, ExcNotImplemented());
}
template <class PolynomialType, int dim, int spacedim>
void
FE_DGVector<PolynomialType,dim,spacedim>::interpolate
(std::vector<double> & /*local_dofs*/,
const std::vector<Vector<double> > & /*values*/,
unsigned int /*offset*/) const
{
Assert(false, ExcNotImplemented());
}
template <class PolynomialType, int dim, int spacedim>
void
FE_DGVector<PolynomialType,dim,spacedim>::interpolate
(std::vector<double> & /*local_dofs*/,
const VectorSlice<const std::vector<std::vector<double> > > & /*values*/) const
{
Assert(false, ExcNotImplemented());
}
template <class PolynomialType, int dim, int spacedim>
std::size_t
FE_DGVector<PolynomialType,dim,spacedim>::memory_consumption() const
{
Assert(false, ExcNotImplemented());
return 0;
}
DEAL_II_NAMESPACE_CLOSE
|