|
| 1 | +// ------------------------------------------------------------------------ |
| 2 | +// |
| 3 | +// SPDX-License-Identifier: LGPL-2.1-or-later |
| 4 | +// Copyright (C) 2007 - 2020 by the deal.II authors |
| 5 | +// |
| 6 | +// This file is part of the deal.II library. |
| 7 | +// |
| 8 | +// Part of the source code is dual licensed under Apache-2.0 WITH |
| 9 | +// LLVM-exception OR LGPL-2.1-or-later. Detailed license information |
| 10 | +// governing the source code and code contributions can be found in |
| 11 | +// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II. |
| 12 | +// |
| 13 | +// ------------------------------------------------------------------------ |
| 14 | + |
| 15 | +// Check gradients in support points for FE_Q_iso_Q1 |
| 16 | + |
| 17 | +#include "../tests.h" |
| 18 | + |
| 19 | +#include <deal.II/base/logstream.h> |
| 20 | +#include <deal.II/base/qprojector.h> |
| 21 | +#include <deal.II/base/quadrature_lib.h> |
| 22 | + |
| 23 | +#include <deal.II/dofs/dof_handler.h> |
| 24 | +#include <deal.II/grid/grid_generator.h> |
| 25 | +#include <deal.II/grid/reference_cell.h> |
| 26 | +#include <deal.II/grid/tria.h> |
| 27 | +#include <deal.II/grid/tria_iterator.h> |
| 28 | + |
| 29 | +#include <deal.II/fe/fe_q_iso_q1.h> |
| 30 | +#include <deal.II/fe/fe_values.h> |
| 31 | +#include <deal.II/fe/mapping_q.h> |
| 32 | + |
| 33 | +template <int dim> |
| 34 | +void check(const FiniteElement<dim> &fe, const Quadrature<dim>& quad, const char* name) |
| 35 | +{ |
| 36 | + deallog << name << ":" << std::endl; |
| 37 | + |
| 38 | + Triangulation<dim> tr; |
| 39 | + MappingQ<dim> mapping(1); |
| 40 | + DoFHandler<dim> dof(tr); |
| 41 | + GridGenerator::hyper_cube(tr, 0., 1.); |
| 42 | + dof.distribute_dofs(fe); |
| 43 | + typename DoFHandler<dim>::cell_iterator c = dof.begin(); |
| 44 | + |
| 45 | + FEValues<dim> fev(mapping, |
| 46 | + fe, |
| 47 | + quad, |
| 48 | + UpdateFlags(update_gradients)); |
| 49 | + |
| 50 | + fev.reinit(c); |
| 51 | + const unsigned int n_q_points = quad.size(); |
| 52 | + |
| 53 | + for (unsigned int q=0;q<n_q_points;++q) |
| 54 | + { |
| 55 | + deallog << quad.point(q); |
| 56 | + for (unsigned int i = 0; i < fe.dofs_per_cell; ++i) { |
| 57 | + deallog << " " << fev.shape_grad(i, q).norm(); |
| 58 | + } |
| 59 | + deallog << std::endl; |
| 60 | + } |
| 61 | + |
| 62 | +} |
| 63 | + |
| 64 | +int |
| 65 | +main() |
| 66 | +{ |
| 67 | + initlog(); |
| 68 | + |
| 69 | + FE_Q_iso_Q1<2> f1(2); |
| 70 | + FE_Q_iso_Q1<2> f2(QGaussLobatto<1>(2 + 1).get_points()); |
| 71 | + QGauss<2> quadrature(3); |
| 72 | + |
| 73 | + check<>(f1, quadrature, "iso_q1(2)-subd"); |
| 74 | + check<>(f2, quadrature, "iso_q1(2)-GL"); |
| 75 | +} |
0 commit comments