34template <
int dim,
int spacedim>
38 Polynomials::generate_complete_Lagrange_basis_on_subdivisions(
48 ExcMessage(
"This element can only be used with a positive number of "
57 const std::vector<unsigned int> R =
58 FETools::lexicographic_to_hierarchic_numbering<dim>(subdivisions);
67 const int N1d = points.size();
68 for (
unsigned int i = 0; i < N; ++i)
69 for (
unsigned int j = 0; j < N; ++j)
76 for (
unsigned int d = 0; d < dim; ++d)
78 int current_distance =
std::abs((xi % N1d) - (xj % N1d));
81 distance =
std::max(distance, current_distance);
95template <
int dim,
int spacedim>
97 const std::vector<
Point<1>> &support_points)
100 Polynomials::generate_complete_linear_basis_on_subdivisions(
104 support_points.
size() - 1,
108 Assert(support_points.size() > 1,
109 ExcMessage(
"This element can only be used with a positive number of "
117template <
int dim,
int spacedim>
125 std::ostringstream namebuf;
127 << this->degree <<
")";
128 return namebuf.str();
133template <
int dim,
int spacedim>
138 std::vector<double> &nodal_values)
const
141 this->get_unit_support_points().size());
145 for (
unsigned int i = 0; i < this->n_dofs_per_cell(); ++i)
149 nodal_values[i] = support_point_values[i](0);
155template <
int dim,
int spacedim>
156std::unique_ptr<FiniteElement<dim, spacedim>>
159 return std::make_unique<FE_Q_iso_Q1<dim, spacedim>>(*this);
164template <
int dim,
int spacedim>
168 const unsigned int codim)
const
191 if (this->degree < fe_q_iso_q1_other->degree &&
192 fe_q_iso_q1_other->degree % this->degree == 0)
194 else if (this->degree == fe_q_iso_q1_other->degree)
196 else if (this->degree > fe_q_iso_q1_other->degree &&
197 this->degree % fe_q_iso_q1_other->degree == 0)
205 if (fe_nothing->is_dominating())
221#include "fe_q_iso_q1.inst"
void initialize(const std::vector< Point< 1 > > &support_points_1d)
virtual FiniteElementDomination::Domination compare_for_domination(const FiniteElement< dim, spacedim > &fe_other, const unsigned int codim=0) const override final
virtual std::string get_name() const override
virtual void convert_generalized_support_point_values_to_dof_values(const std::vector< Vector< double > > &support_point_values, std::vector< double > &nodal_values) const override
virtual std::unique_ptr< FiniteElement< dim, spacedim > > clone() const override
FE_Q_iso_Q1(const unsigned int n_subdivisions)
unsigned int n_dofs_per_cell() const
Table< 2, bool > local_dof_sparsity_pattern
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
static ::ExceptionBase & ExcImpossibleInDim(int arg1)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define DEAL_II_NOT_IMPLEMENTED()
@ either_element_can_dominate
@ other_element_dominates
@ neither_element_dominates
std::string dim_string(const int dim, const int spacedim)
constexpr T pow(const T base, const int iexp)
::VectorizedArray< Number, width > max(const ::VectorizedArray< Number, width > &, const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > abs(const ::VectorizedArray< Number, width > &)