From b0c44f0dc2753a717ecc52fe018df316dd80e141 Mon Sep 17 00:00:00 2001 From: Jaeryun Yim Date: Thu, 5 Nov 2015 12:12:50 +0900 Subject: [PATCH] Add a first test for the P1NC element. --- tests/fe/p1nc_parametric_01.cc | 77 ++++++++++++++++++++++++++++++ tests/fe/p1nc_parametric_01.output | 10 ++++ 2 files changed, 87 insertions(+) create mode 100644 tests/fe/p1nc_parametric_01.cc create mode 100644 tests/fe/p1nc_parametric_01.output diff --git a/tests/fe/p1nc_parametric_01.cc b/tests/fe/p1nc_parametric_01.cc new file mode 100644 index 0000000000..436d30f4ff --- /dev/null +++ b/tests/fe/p1nc_parametric_01.cc @@ -0,0 +1,77 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2015 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. +// +// --------------------------------------------------------------------- + + + +// check the P1NC parametric element on a square + +#include "../tests.h" +#include +#include +#include +#include +#include +#include + +#include + +#include +#include + +#define PRECISION 5 + + + +template +void +check() +{ + Triangulation triangulation; + GridGenerator::hyper_cube (triangulation, 0, 5); + + FE_P1NCParametric fe; + DoFHandler dof_handler (triangulation); + dof_handler.distribute_dofs(fe); + + QGauss quadrature(3); + FEValues fe_values (fe, quadrature, update_values | update_q_points); + fe_values.reinit (dof_handler.begin_active()); + + for (unsigned int q=0; q(); +} + + + diff --git a/tests/fe/p1nc_parametric_01.output b/tests/fe/p1nc_parametric_01.output new file mode 100644 index 0000000000..5436e9b991 --- /dev/null +++ b/tests/fe/p1nc_parametric_01.output @@ -0,0 +1,10 @@ + +DEAL::index=0 position=0.56351 0.56351 values=0.63730 0.25000 0.25000 -0.13730 +DEAL::index=1 position=2.50000 0.56351 values=0.44365 0.44365 0.05635 0.05635 +DEAL::index=2 position=4.43649 0.56351 values=0.25000 0.63730 -0.13730 0.25000 +DEAL::index=3 position=0.56351 2.50000 values=0.44365 0.05635 0.44365 0.05635 +DEAL::index=4 position=2.50000 2.50000 values=0.25000 0.25000 0.25000 0.25000 +DEAL::index=5 position=4.43649 2.50000 values=0.05635 0.44365 0.05635 0.44365 +DEAL::index=6 position=0.56351 4.43649 values=0.25000 -0.13730 0.63730 0.25000 +DEAL::index=7 position=2.50000 4.43649 values=0.05635 0.05635 0.44365 0.44365 +DEAL::index=8 position=4.43649 4.43649 values=-0.13730 0.25000 0.25000 0.63730 -- 2.39.5