Dear all,
I like to create a three dimensional Lagrange finite element that uses
ansatz functions of the following form: N_i (x,y,z) = q_i(x,y) p_i (z),
where q_i (x,y) are the basis functions of a two dimensional Lagrange
finite element of arbitrary degree, basically FE_Q<2> (degree) and p_i are
the basis functions of a linear Lagrange finite element FE<1> (1).
The class AnisotropicPolynomials<dim> provides the functionality to exactly
create the desired ansatz functions. However, I'm not sure how to use this
class to create a respective Finite Element class.
Therefore, I checked the code in fe_q.cc and the constructor is implemented
as follows.
52 template <int dim, int spacedim>
53 FE_Q<dim,spacedim>::FE_Q (const unsigned int degree)
54 :
55 FE_Q_Base<TensorProductPolynomials<dim>, dim, spacedim>
56 (TensorProductPolynomials<dim>(Polynomials::
generate_complete_Lagrange_basis(internal::FE_Q::get_QGaussLobatto_points(
degree))),
57 FiniteElementData<dim>(this->get_dpo_vector(degree),
58 1, degree,
59 FiniteElementData<dim>::H1),
60 std::vector<bool> (1, false))
61 {
62 this->initialize(internal::FE_Q::get_QGaussLobatto_points(degree));
63 }
Would it be possible to replace the TensorProductPolynomial class by the
AnisotropicPolynomials class?
Best regards,
Dustin
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
