Reference documentation for deal.II version Git 840703ecf3 2019-10-14 17:15:20 -0400
\(\newcommand{\dealcoloneq}{\mathrel{\vcenter{:}}=}\)
polynomials_bernardi_raugel.h
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 2004 - 2018 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE.md at
12 // the top level directory of deal.II.
13 //
14 // ---------------------------------------------------------------------
15 
16 
17 #ifndef dealii_polynomials_bernardi_raugel_h
18 #define dealii_polynomials_bernardi_raugel_h
19 
20 #include <deal.II/base/geometry_info.h>
21 #include <deal.II/base/point.h>
22 #include <deal.II/base/polynomial.h>
23 #include <deal.II/base/polynomial_space.h>
24 #include <deal.II/base/tensor.h>
25 #include <deal.II/base/tensor_polynomials_base.h>
26 #include <deal.II/base/tensor_product_polynomials.h>
27 
28 #include <vector>
29 
30 DEAL_II_NAMESPACE_OPEN
31 
32 
84 template <int dim>
86 {
87 public:
95  PolynomialsBernardiRaugel(const unsigned int k);
96 
100  std::string
101  name() const override;
102 
115  void
116  evaluate(const Point<dim> & unit_point,
117  std::vector<Tensor<1, dim>> &values,
118  std::vector<Tensor<2, dim>> &grads,
119  std::vector<Tensor<3, dim>> &grad_grads,
120  std::vector<Tensor<4, dim>> &third_derivatives,
121  std::vector<Tensor<5, dim>> &fourth_derivatives) const override;
122 
128  static unsigned int
129  n_polynomials(const unsigned int k);
130 
134  virtual std::unique_ptr<TensorPolynomialsBase<dim>>
135  clone() const override;
136 
137 private:
145 
152 
157  static std::vector<std::vector<Polynomials::Polynomial<double>>>
159 
164  static std::vector<std::vector<Polynomials::Polynomial<double>>>
166 };
167 
168 
169 template <int dim>
170 inline std::string
172 {
173  return "BernardiRaugel";
174 }
175 
176 
177 DEAL_II_NAMESPACE_CLOSE
178 
179 #endif
void evaluate(const Point< dim > &unit_point, std::vector< Tensor< 1, dim >> &values, std::vector< Tensor< 2, dim >> &grads, std::vector< Tensor< 3, dim >> &grad_grads, std::vector< Tensor< 4, dim >> &third_derivatives, std::vector< Tensor< 5, dim >> &fourth_derivatives) const override
std::string name() const override
const AnisotropicPolynomials< dim > polynomial_space_Q
static unsigned int n_polynomials(const unsigned int k)
virtual std::unique_ptr< TensorPolynomialsBase< dim > > clone() const override
static std::vector< std::vector< Polynomials::Polynomial< double > > > create_polynomials_Q()
static std::vector< std::vector< Polynomials::Polynomial< double > > > create_polynomials_bubble()
PolynomialsBernardiRaugel(const unsigned int k)
const AnisotropicPolynomials< dim > polynomial_space_bubble