From 3c2df429b3451fe9c350f9a2a9d72b0652e14b7d Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 14 Sep 2012 12:32:17 +0000 Subject: [PATCH] Step one in breaking apart our monster testcase 'internals'. That testcase had an output of over 6MB with more than 12000 lines, a good fraction of which is susceptible to round-off errors. It is time to split this testcase by finite element to be able to tell more easily what is going on. These new testcases take over the part where we test support points. git-svn-id: https://svn.dealii.org/trunk@26376 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/fe/fe_support_points_bdm.cc | 29 +++ tests/fe/fe_support_points_bdm/cmp/generic | 11 + tests/fe/fe_support_points_common.h | 82 +++++++ tests/fe/fe_support_points_nedelec.cc | 31 +++ .../fe/fe_support_points_nedelec/cmp/generic | 25 +++ tests/fe/fe_support_points_q.cc | 33 +++ tests/fe/fe_support_points_q/cmp/generic | 208 ++++++++++++++++++ tests/fe/fe_support_points_q_hierarchical.cc | 32 +++ .../cmp/generic | 208 ++++++++++++++++++ tests/fe/fe_support_points_rt.cc | 31 +++ tests/fe/fe_support_points_rt/cmp/generic | 25 +++ tests/fe/fe_support_points_sys_01.cc | 28 +++ tests/fe/fe_support_points_sys_01/cmp/generic | 42 ++++ tests/fe/fe_support_points_sys_02.cc | 28 +++ tests/fe/fe_support_points_sys_02/cmp/generic | 24 ++ tests/fe/fe_support_points_sys_03.cc | 28 +++ tests/fe/fe_support_points_sys_03/cmp/generic | 60 +++++ tests/fe/fe_support_points_sys_04.cc | 28 +++ tests/fe/fe_support_points_sys_04/cmp/generic | 138 ++++++++++++ tests/fe/fe_support_points_sys_05.cc | 31 +++ tests/fe/fe_support_points_sys_05/cmp/generic | 166 ++++++++++++++ tests/fe/fe_support_points_sys_06.cc | 32 +++ tests/fe/fe_support_points_sys_06/cmp/generic | 198 +++++++++++++++++ tests/fe/fe_support_points_sys_07.cc | 30 +++ tests/fe/fe_support_points_sys_07/cmp/generic | 6 + tests/fe/fe_support_points_sys_08.cc | 32 +++ tests/fe/fe_support_points_sys_08/cmp/generic | 6 + tests/fe/fe_support_points_sys_09.cc | 30 +++ tests/fe/fe_support_points_sys_09/cmp/generic | 6 + tests/fe/fe_support_points_sys_10.cc | 31 +++ tests/fe/fe_support_points_sys_10/cmp/generic | 6 + 31 files changed, 1665 insertions(+) create mode 100644 tests/fe/fe_support_points_bdm.cc create mode 100644 tests/fe/fe_support_points_bdm/cmp/generic create mode 100644 tests/fe/fe_support_points_common.h create mode 100644 tests/fe/fe_support_points_nedelec.cc create mode 100644 tests/fe/fe_support_points_nedelec/cmp/generic create mode 100644 tests/fe/fe_support_points_q.cc create mode 100644 tests/fe/fe_support_points_q/cmp/generic create mode 100644 tests/fe/fe_support_points_q_hierarchical.cc create mode 100644 tests/fe/fe_support_points_q_hierarchical/cmp/generic create mode 100644 tests/fe/fe_support_points_rt.cc create mode 100644 tests/fe/fe_support_points_rt/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_01.cc create mode 100644 tests/fe/fe_support_points_sys_01/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_02.cc create mode 100644 tests/fe/fe_support_points_sys_02/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_03.cc create mode 100644 tests/fe/fe_support_points_sys_03/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_04.cc create mode 100644 tests/fe/fe_support_points_sys_04/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_05.cc create mode 100644 tests/fe/fe_support_points_sys_05/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_06.cc create mode 100644 tests/fe/fe_support_points_sys_06/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_07.cc create mode 100644 tests/fe/fe_support_points_sys_07/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_08.cc create mode 100644 tests/fe/fe_support_points_sys_08/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_09.cc create mode 100644 tests/fe/fe_support_points_sys_09/cmp/generic create mode 100644 tests/fe/fe_support_points_sys_10.cc create mode 100644 tests/fe/fe_support_points_sys_10/cmp/generic diff --git a/tests/fe/fe_support_points_bdm.cc b/tests/fe/fe_support_points_bdm.cc new file mode 100644 index 0000000000..ec505ba1ff --- /dev/null +++ b/tests/fe/fe_support_points_bdm.cc @@ -0,0 +1,29 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_ALL(BDM, 1, 2); + CHECK_ALL(BDM, 2, 2); +} diff --git a/tests/fe/fe_support_points_bdm/cmp/generic b/tests/fe/fe_support_points_bdm/cmp/generic new file mode 100644 index 0000000000..65f9349cf0 --- /dev/null +++ b/tests/fe/fe_support_points_bdm/cmp/generic @@ -0,0 +1,11 @@ + +DEAL::BDM1<2> cell support points +DEAL::BDM1<2> face 0 support points +DEAL::BDM1<2> face 1 support points +DEAL::BDM1<2> face 2 support points +DEAL::BDM1<2> face 3 support points +DEAL::BDM2<2> cell support points +DEAL::BDM2<2> face 0 support points +DEAL::BDM2<2> face 1 support points +DEAL::BDM2<2> face 2 support points +DEAL::BDM2<2> face 3 support points diff --git a/tests/fe/fe_support_points_common.h b/tests/fe/fe_support_points_common.h new file mode 100644 index 0000000000..e5ff06123b --- /dev/null +++ b/tests/fe/fe_support_points_common.h @@ -0,0 +1,82 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + + + +template +inline void +check_support (const FiniteElement& finel, const char* name) +{ + Triangulation tr; + GridGenerator::hyper_cube(tr, 0., 1.); + DoFHandler dof (tr); + dof.distribute_dofs (finel); + + deallog << name << '<' << dim << '>' << " cell support points" << std::endl; + + const std::vector > &cell_points = finel.get_unit_support_points (); + + for (unsigned int k=0;k > &face_points = finel.get_unit_face_support_points (); + const std::vector dummy_weights (face_points.size()); + + Quadrature q(face_points, dummy_weights); + + for (unsigned int i=0;i::faces_per_cell;++i) + { + std::vector > q_points (q.get_points().size()); + QProjector::project_to_face (q, i, q_points); + Quadrature qp(q_points); + deallog << name << '<' << dim << '>' << " face " << i + << " support points" << std::endl; + + for (unsigned int k=0; k EL(deg); check_support(EL, #EL #deg); } +#define CHECK_SYS1(sub1,N1,dim) { FESystem q(sub1, N1); check_support(q, #sub1 #N1); } +#define CHECK_SYS2(sub1,N1,sub2,N2,dim) { FESystem q(sub1, N1, sub2, N2); \ + check_support(q, #sub1 #N1 #sub2 #N2); } +#define CHECK_SYS3(sub1,N1,sub2,N2,sub3,N3,dim) { FESystem q(sub1, N1, sub2, N2, sub3, N3); \ + check_support(q, #sub1 #N1 #sub2 #N2 #sub3 #N3); } + + diff --git a/tests/fe/fe_support_points_nedelec.cc b/tests/fe/fe_support_points_nedelec.cc new file mode 100644 index 0000000000..e7eef71d71 --- /dev/null +++ b/tests/fe/fe_support_points_nedelec.cc @@ -0,0 +1,31 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_ALL(Nedelec, 0, 2); + CHECK_ALL(Nedelec, 0, 3); + CHECK_ALL(Nedelec, 1, 2); + CHECK_ALL(Nedelec, 1, 3); +} diff --git a/tests/fe/fe_support_points_nedelec/cmp/generic b/tests/fe/fe_support_points_nedelec/cmp/generic new file mode 100644 index 0000000000..05ead30500 --- /dev/null +++ b/tests/fe/fe_support_points_nedelec/cmp/generic @@ -0,0 +1,25 @@ + +DEAL::Nedelec0<2> cell support points +DEAL::Nedelec0<2> face 0 support points +DEAL::Nedelec0<2> face 1 support points +DEAL::Nedelec0<2> face 2 support points +DEAL::Nedelec0<2> face 3 support points +DEAL::Nedelec0<3> cell support points +DEAL::Nedelec0<3> face 0 support points +DEAL::Nedelec0<3> face 1 support points +DEAL::Nedelec0<3> face 2 support points +DEAL::Nedelec0<3> face 3 support points +DEAL::Nedelec0<3> face 4 support points +DEAL::Nedelec0<3> face 5 support points +DEAL::Nedelec1<2> cell support points +DEAL::Nedelec1<2> face 0 support points +DEAL::Nedelec1<2> face 1 support points +DEAL::Nedelec1<2> face 2 support points +DEAL::Nedelec1<2> face 3 support points +DEAL::Nedelec1<3> cell support points +DEAL::Nedelec1<3> face 0 support points +DEAL::Nedelec1<3> face 1 support points +DEAL::Nedelec1<3> face 2 support points +DEAL::Nedelec1<3> face 3 support points +DEAL::Nedelec1<3> face 4 support points +DEAL::Nedelec1<3> face 5 support points diff --git a/tests/fe/fe_support_points_q.cc b/tests/fe/fe_support_points_q.cc new file mode 100644 index 0000000000..6818322071 --- /dev/null +++ b/tests/fe/fe_support_points_q.cc @@ -0,0 +1,33 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_ALL(Q,1,2); + CHECK_ALL(Q,2,2); + CHECK_ALL(Q,3,2); + + CHECK_ALL(Q,1,3); + CHECK_ALL(Q,2,3); +} diff --git a/tests/fe/fe_support_points_q/cmp/generic b/tests/fe/fe_support_points_q/cmp/generic new file mode 100644 index 0000000000..ef1829c25b --- /dev/null +++ b/tests/fe/fe_support_points_q/cmp/generic @@ -0,0 +1,208 @@ + +DEAL::Q1<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::Q1<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::Q1<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::Q1<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::Q1<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::Q2<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.500 +DEAL::1.00 0.500 +DEAL::0.500 0.00 +DEAL::0.500 1.00 +DEAL::0.500 0.500 +DEAL::Q2<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 0.500 +DEAL::Q2<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 0.500 +DEAL::Q2<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.500 0.00 +DEAL::Q2<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.500 1.00 +DEAL::Q3<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::Q3<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::Q3<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::Q3<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::Q3<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::Q1<3> cell support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q1<3> face 0 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::Q1<3> face 1 support points +DEAL::1.00 0.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q1<3> face 2 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::Q1<3> face 3 support points +DEAL::0.00 1.00 0.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 1.00 1.00 +DEAL::Q1<3> face 4 support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::Q1<3> face 5 support points +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q2<3> cell support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::0.00 0.500 0.00 +DEAL::1.00 0.500 0.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 1.00 0.00 +DEAL::0.00 0.500 1.00 +DEAL::1.00 0.500 1.00 +DEAL::0.500 0.00 1.00 +DEAL::0.500 1.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::1.00 0.00 0.500 +DEAL::0.00 1.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::0.00 0.500 0.500 +DEAL::1.00 0.500 0.500 +DEAL::0.500 0.00 0.500 +DEAL::0.500 1.00 0.500 +DEAL::0.500 0.500 0.00 +DEAL::0.500 0.500 1.00 +DEAL::0.500 0.500 0.500 +DEAL::Q2<3> face 0 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::0.00 1.00 0.500 +DEAL::0.00 0.500 0.00 +DEAL::0.00 0.500 1.00 +DEAL::0.00 0.500 0.500 +DEAL::Q2<3> face 1 support points +DEAL::1.00 0.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::1.00 0.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::1.00 0.500 0.00 +DEAL::1.00 0.500 1.00 +DEAL::1.00 0.500 0.500 +DEAL::Q2<3> face 2 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 0.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::1.00 0.00 0.500 +DEAL::0.500 0.00 0.500 +DEAL::Q2<3> face 3 support points +DEAL::0.00 1.00 0.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 1.00 1.00 +DEAL::0.500 1.00 0.00 +DEAL::0.500 1.00 1.00 +DEAL::0.00 1.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::0.500 1.00 0.500 +DEAL::Q2<3> face 4 support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.500 0.00 +DEAL::1.00 0.500 0.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 1.00 0.00 +DEAL::0.500 0.500 0.00 +DEAL::Q2<3> face 5 support points +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::0.00 0.500 1.00 +DEAL::1.00 0.500 1.00 +DEAL::0.500 0.00 1.00 +DEAL::0.500 1.00 1.00 +DEAL::0.500 0.500 1.00 diff --git a/tests/fe/fe_support_points_q_hierarchical.cc b/tests/fe/fe_support_points_q_hierarchical.cc new file mode 100644 index 0000000000..cc26fffc6f --- /dev/null +++ b/tests/fe/fe_support_points_q_hierarchical.cc @@ -0,0 +1,32 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_ALL(Q_Hierarchical,1,2); + CHECK_ALL(Q_Hierarchical,2,2); + CHECK_ALL(Q_Hierarchical,3,2); + CHECK_ALL(Q_Hierarchical,1,3); + CHECK_ALL(Q_Hierarchical,2,3); +} diff --git a/tests/fe/fe_support_points_q_hierarchical/cmp/generic b/tests/fe/fe_support_points_q_hierarchical/cmp/generic new file mode 100644 index 0000000000..a240ce5ddc --- /dev/null +++ b/tests/fe/fe_support_points_q_hierarchical/cmp/generic @@ -0,0 +1,208 @@ + +DEAL::Q_Hierarchical1<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::Q_Hierarchical1<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::Q_Hierarchical1<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::Q_Hierarchical1<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::Q_Hierarchical1<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::Q_Hierarchical2<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.500 +DEAL::1.00 0.500 +DEAL::0.500 0.00 +DEAL::0.500 1.00 +DEAL::0.500 0.500 +DEAL::Q_Hierarchical2<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 0.500 +DEAL::Q_Hierarchical2<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 0.500 +DEAL::Q_Hierarchical2<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.500 0.00 +DEAL::Q_Hierarchical2<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.500 1.00 +DEAL::Q_Hierarchical3<2> cell support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::Q_Hierarchical3<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::Q_Hierarchical3<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::Q_Hierarchical3<2> face 2 support points +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::Q_Hierarchical3<2> face 3 support points +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::Q_Hierarchical1<3> cell support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q_Hierarchical1<3> face 0 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::Q_Hierarchical1<3> face 1 support points +DEAL::1.00 0.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q_Hierarchical1<3> face 2 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::Q_Hierarchical1<3> face 3 support points +DEAL::0.00 1.00 0.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 1.00 1.00 +DEAL::Q_Hierarchical1<3> face 4 support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::Q_Hierarchical1<3> face 5 support points +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::Q_Hierarchical2<3> cell support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::0.00 0.500 0.00 +DEAL::1.00 0.500 0.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 1.00 0.00 +DEAL::0.00 0.500 1.00 +DEAL::1.00 0.500 1.00 +DEAL::0.500 0.00 1.00 +DEAL::0.500 1.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::1.00 0.00 0.500 +DEAL::0.00 1.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::0.00 0.500 0.500 +DEAL::1.00 0.500 0.500 +DEAL::0.500 0.00 0.500 +DEAL::0.500 1.00 0.500 +DEAL::0.500 0.500 0.00 +DEAL::0.500 0.500 1.00 +DEAL::0.500 0.500 0.500 +DEAL::Q_Hierarchical2<3> face 0 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::0.00 1.00 0.500 +DEAL::0.00 0.500 0.00 +DEAL::0.00 0.500 1.00 +DEAL::0.00 0.500 0.500 +DEAL::Q_Hierarchical2<3> face 1 support points +DEAL::1.00 0.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::1.00 0.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::1.00 0.500 0.00 +DEAL::1.00 0.500 1.00 +DEAL::1.00 0.500 0.500 +DEAL::Q_Hierarchical2<3> face 2 support points +DEAL::0.00 0.00 0.00 +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 0.00 +DEAL::1.00 0.00 1.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 0.00 1.00 +DEAL::0.00 0.00 0.500 +DEAL::1.00 0.00 0.500 +DEAL::0.500 0.00 0.500 +DEAL::Q_Hierarchical2<3> face 3 support points +DEAL::0.00 1.00 0.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 0.00 +DEAL::1.00 1.00 1.00 +DEAL::0.500 1.00 0.00 +DEAL::0.500 1.00 1.00 +DEAL::0.00 1.00 0.500 +DEAL::1.00 1.00 0.500 +DEAL::0.500 1.00 0.500 +DEAL::Q_Hierarchical2<3> face 4 support points +DEAL::0.00 0.00 0.00 +DEAL::1.00 0.00 0.00 +DEAL::0.00 1.00 0.00 +DEAL::1.00 1.00 0.00 +DEAL::0.00 0.500 0.00 +DEAL::1.00 0.500 0.00 +DEAL::0.500 0.00 0.00 +DEAL::0.500 1.00 0.00 +DEAL::0.500 0.500 0.00 +DEAL::Q_Hierarchical2<3> face 5 support points +DEAL::0.00 0.00 1.00 +DEAL::1.00 0.00 1.00 +DEAL::0.00 1.00 1.00 +DEAL::1.00 1.00 1.00 +DEAL::0.00 0.500 1.00 +DEAL::1.00 0.500 1.00 +DEAL::0.500 0.00 1.00 +DEAL::0.500 1.00 1.00 +DEAL::0.500 0.500 1.00 diff --git a/tests/fe/fe_support_points_rt.cc b/tests/fe/fe_support_points_rt.cc new file mode 100644 index 0000000000..dbc6807c1e --- /dev/null +++ b/tests/fe/fe_support_points_rt.cc @@ -0,0 +1,31 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_ALL(RaviartThomas, 0, 2); + CHECK_ALL(RaviartThomas, 0, 3); + CHECK_ALL(RaviartThomas, 1, 2); + CHECK_ALL(RaviartThomas, 1, 3); +} diff --git a/tests/fe/fe_support_points_rt/cmp/generic b/tests/fe/fe_support_points_rt/cmp/generic new file mode 100644 index 0000000000..a60ccf56fa --- /dev/null +++ b/tests/fe/fe_support_points_rt/cmp/generic @@ -0,0 +1,25 @@ + +DEAL::RaviartThomas0<2> cell support points +DEAL::RaviartThomas0<2> face 0 support points +DEAL::RaviartThomas0<2> face 1 support points +DEAL::RaviartThomas0<2> face 2 support points +DEAL::RaviartThomas0<2> face 3 support points +DEAL::RaviartThomas0<3> cell support points +DEAL::RaviartThomas0<3> face 0 support points +DEAL::RaviartThomas0<3> face 1 support points +DEAL::RaviartThomas0<3> face 2 support points +DEAL::RaviartThomas0<3> face 3 support points +DEAL::RaviartThomas0<3> face 4 support points +DEAL::RaviartThomas0<3> face 5 support points +DEAL::RaviartThomas1<2> cell support points +DEAL::RaviartThomas1<2> face 0 support points +DEAL::RaviartThomas1<2> face 1 support points +DEAL::RaviartThomas1<2> face 2 support points +DEAL::RaviartThomas1<2> face 3 support points +DEAL::RaviartThomas1<3> cell support points +DEAL::RaviartThomas1<3> face 0 support points +DEAL::RaviartThomas1<3> face 1 support points +DEAL::RaviartThomas1<3> face 2 support points +DEAL::RaviartThomas1<3> face 3 support points +DEAL::RaviartThomas1<3> face 4 support points +DEAL::RaviartThomas1<3> face 5 support points diff --git a/tests/fe/fe_support_points_sys_01.cc b/tests/fe/fe_support_points_sys_01.cc new file mode 100644 index 0000000000..cec4ebcfb0 --- /dev/null +++ b/tests/fe/fe_support_points_sys_01.cc @@ -0,0 +1,28 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS1(FE_Q<2>(1),3,2); +} diff --git a/tests/fe/fe_support_points_sys_01/cmp/generic b/tests/fe/fe_support_points_sys_01/cmp/generic new file mode 100644 index 0000000000..514de0a89f --- /dev/null +++ b/tests/fe/fe_support_points_sys_01/cmp/generic @@ -0,0 +1,42 @@ + +DEAL::FE_Q<2>(1)3<2> cell support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::FE_Q<2>(1)3<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::FE_Q<2>(1)3<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::FE_Q<2>(1)3<2> face 2 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::FE_Q<2>(1)3<2> face 3 support points +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 diff --git a/tests/fe/fe_support_points_sys_02.cc b/tests/fe/fe_support_points_sys_02.cc new file mode 100644 index 0000000000..ba270a88df --- /dev/null +++ b/tests/fe/fe_support_points_sys_02.cc @@ -0,0 +1,28 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS1(FE_DGQ<2>(2),2,2); +} diff --git a/tests/fe/fe_support_points_sys_02/cmp/generic b/tests/fe/fe_support_points_sys_02/cmp/generic new file mode 100644 index 0000000000..c5ff501a95 --- /dev/null +++ b/tests/fe/fe_support_points_sys_02/cmp/generic @@ -0,0 +1,24 @@ + +DEAL::FE_DGQ<2>(2)2<2> cell support points +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::FE_DGQ<2>(2)2<2> face 0 support points +DEAL::FE_DGQ<2>(2)2<2> face 1 support points +DEAL::FE_DGQ<2>(2)2<2> face 2 support points +DEAL::FE_DGQ<2>(2)2<2> face 3 support points diff --git a/tests/fe/fe_support_points_sys_03.cc b/tests/fe/fe_support_points_sys_03.cc new file mode 100644 index 0000000000..716ddc716c --- /dev/null +++ b/tests/fe/fe_support_points_sys_03.cc @@ -0,0 +1,28 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS2(FE_Q<2>(1),3,FE_DGQ<2>(2),2,2); +} diff --git a/tests/fe/fe_support_points_sys_03/cmp/generic b/tests/fe/fe_support_points_sys_03/cmp/generic new file mode 100644 index 0000000000..ec8e668f21 --- /dev/null +++ b/tests/fe/fe_support_points_sys_03/cmp/generic @@ -0,0 +1,60 @@ + +DEAL::FE_Q<2>(1)3FE_DGQ<2>(2)2<2> cell support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::FE_Q<2>(1)3FE_DGQ<2>(2)2<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::FE_Q<2>(1)3FE_DGQ<2>(2)2<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::FE_Q<2>(1)3FE_DGQ<2>(2)2<2> face 2 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::FE_Q<2>(1)3FE_DGQ<2>(2)2<2> face 3 support points +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 diff --git a/tests/fe/fe_support_points_sys_04.cc b/tests/fe/fe_support_points_sys_04.cc new file mode 100644 index 0000000000..46858df5dd --- /dev/null +++ b/tests/fe/fe_support_points_sys_04.cc @@ -0,0 +1,28 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS3(FE_DGQ<2>(2),2,FE_DGQ<2>(2),2,FE_Q<2>(3),3,2); +} diff --git a/tests/fe/fe_support_points_sys_04/cmp/generic b/tests/fe/fe_support_points_sys_04/cmp/generic new file mode 100644 index 0000000000..1f19a2ce0c --- /dev/null +++ b/tests/fe/fe_support_points_sys_04/cmp/generic @@ -0,0 +1,138 @@ + +DEAL::FE_DGQ<2>(2)2FE_DGQ<2>(2)2FE_Q<2>(3)3<2> cell support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.500 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.500 +DEAL::0.500 0.500 +DEAL::1.00 0.500 +DEAL::0.00 1.00 +DEAL::0.500 1.00 +DEAL::1.00 1.00 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::FE_DGQ<2>(2)2FE_DGQ<2>(2)2FE_Q<2>(3)3<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::0.00 0.333 +DEAL::0.00 0.667 +DEAL::FE_DGQ<2>(2)2FE_DGQ<2>(2)2FE_Q<2>(3)3<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::1.00 0.333 +DEAL::1.00 0.667 +DEAL::FE_DGQ<2>(2)2FE_DGQ<2>(2)2FE_Q<2>(3)3<2> face 2 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::FE_DGQ<2>(2)2FE_DGQ<2>(2)2FE_Q<2>(3)3<2> face 3 support points +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 diff --git a/tests/fe/fe_support_points_sys_05.cc b/tests/fe/fe_support_points_sys_05.cc new file mode 100644 index 0000000000..36f173ebdd --- /dev/null +++ b/tests/fe/fe_support_points_sys_05.cc @@ -0,0 +1,31 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS3((FESystem<2>(FE_Q<2>(1),3)), 3, + FE_DGQ<2>(3), 1, + FE_Q<2>(1), 3, + 2); +} diff --git a/tests/fe/fe_support_points_sys_05/cmp/generic b/tests/fe/fe_support_points_sys_05/cmp/generic new file mode 100644 index 0000000000..53fbc315e6 --- /dev/null +++ b/tests/fe/fe_support_points_sys_05/cmp/generic @@ -0,0 +1,166 @@ + +DEAL::(FESystem<2>(FE_Q<2>(1),3))3FE_DGQ<2>(3)1FE_Q<2>(1)3<2> cell support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.333 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::1.00 0.333 +DEAL::0.00 0.667 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::1.00 0.667 +DEAL::0.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::1.00 1.00 +DEAL::(FESystem<2>(FE_Q<2>(1),3))3FE_DGQ<2>(3)1FE_Q<2>(1)3<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::(FESystem<2>(FE_Q<2>(1),3))3FE_DGQ<2>(3)1FE_Q<2>(1)3<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::(FESystem<2>(FE_Q<2>(1),3))3FE_DGQ<2>(3)1FE_Q<2>(1)3<2> face 2 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::(FESystem<2>(FE_Q<2>(1),3))3FE_DGQ<2>(3)1FE_Q<2>(1)3<2> face 3 support points +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 diff --git a/tests/fe/fe_support_points_sys_06.cc b/tests/fe/fe_support_points_sys_06.cc new file mode 100644 index 0000000000..38f12055c0 --- /dev/null +++ b/tests/fe/fe_support_points_sys_06.cc @@ -0,0 +1,32 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS3(FE_DGQ<2>(3), 1, + FESystem<2>(FE_DGQ<2>(3),3), 1, + FESystem<2>(FE_Q<2>(2),3, + FE_DGQ<2>(0),1),2, + 2); +} diff --git a/tests/fe/fe_support_points_sys_06/cmp/generic b/tests/fe/fe_support_points_sys_06/cmp/generic new file mode 100644 index 0000000000..18ae3cd530 --- /dev/null +++ b/tests/fe/fe_support_points_sys_06/cmp/generic @@ -0,0 +1,198 @@ + +DEAL::FE_DGQ<2>(3)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_DGQ<2>(0),1)2<2> cell support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.333 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::1.00 0.333 +DEAL::0.00 0.667 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::1.00 0.667 +DEAL::0.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.333 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::1.00 0.333 +DEAL::0.00 0.667 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::1.00 0.667 +DEAL::0.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.333 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::1.00 0.333 +DEAL::0.00 0.667 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::1.00 0.667 +DEAL::0.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::1.00 1.00 +DEAL::0.00 0.00 +DEAL::0.333 0.00 +DEAL::0.667 0.00 +DEAL::1.00 0.00 +DEAL::0.00 0.333 +DEAL::0.333 0.333 +DEAL::0.667 0.333 +DEAL::1.00 0.333 +DEAL::0.00 0.667 +DEAL::0.333 0.667 +DEAL::0.667 0.667 +DEAL::1.00 0.667 +DEAL::0.00 1.00 +DEAL::0.333 1.00 +DEAL::0.667 1.00 +DEAL::1.00 1.00 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::0.500 0.500 +DEAL::FE_DGQ<2>(3)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_DGQ<2>(0),1)2<2> face 0 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::0.00 0.500 +DEAL::FE_DGQ<2>(3)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_DGQ<2>(0),1)2<2> face 1 support points +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::1.00 0.500 +DEAL::FE_DGQ<2>(3)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_DGQ<2>(0),1)2<2> face 2 support points +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::0.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::1.00 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::0.500 0.00 +DEAL::FE_DGQ<2>(3)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_DGQ<2>(0),1)2<2> face 3 support points +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::0.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::1.00 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 +DEAL::0.500 1.00 diff --git a/tests/fe/fe_support_points_sys_07.cc b/tests/fe/fe_support_points_sys_07.cc new file mode 100644 index 0000000000..efad599057 --- /dev/null +++ b/tests/fe/fe_support_points_sys_07.cc @@ -0,0 +1,30 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS2 (FE_DGQ<2>(3), 1, + FE_Nedelec<2>(0), 2, + 2); +} diff --git a/tests/fe/fe_support_points_sys_07/cmp/generic b/tests/fe/fe_support_points_sys_07/cmp/generic new file mode 100644 index 0000000000..e6140c4398 --- /dev/null +++ b/tests/fe/fe_support_points_sys_07/cmp/generic @@ -0,0 +1,6 @@ + +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(0)2<2> cell support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(0)2<2> face 0 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(0)2<2> face 1 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(0)2<2> face 2 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(0)2<2> face 3 support points diff --git a/tests/fe/fe_support_points_sys_08.cc b/tests/fe/fe_support_points_sys_08.cc new file mode 100644 index 0000000000..6f0b7f64b1 --- /dev/null +++ b/tests/fe/fe_support_points_sys_08.cc @@ -0,0 +1,32 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS3(FE_Nedelec<2>(0), 1, + FESystem<2>(FE_DGQ<2>(3),3), 1, + FESystem<2>(FE_Q<2>(2),3, + FE_Nedelec<2>(0),2),2, + 2); +} diff --git a/tests/fe/fe_support_points_sys_08/cmp/generic b/tests/fe/fe_support_points_sys_08/cmp/generic new file mode 100644 index 0000000000..8ef502f9aa --- /dev/null +++ b/tests/fe/fe_support_points_sys_08/cmp/generic @@ -0,0 +1,6 @@ + +DEAL::FE_Nedelec<2>(0)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(0),2)2<2> cell support points +DEAL::FE_Nedelec<2>(0)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(0),2)2<2> face 0 support points +DEAL::FE_Nedelec<2>(0)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(0),2)2<2> face 1 support points +DEAL::FE_Nedelec<2>(0)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(0),2)2<2> face 2 support points +DEAL::FE_Nedelec<2>(0)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(0),2)2<2> face 3 support points diff --git a/tests/fe/fe_support_points_sys_09.cc b/tests/fe/fe_support_points_sys_09.cc new file mode 100644 index 0000000000..5bb2e9c175 --- /dev/null +++ b/tests/fe/fe_support_points_sys_09.cc @@ -0,0 +1,30 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS2 (FE_DGQ<2>(3), 1, + FE_Nedelec<2>(1), 2, + 2); +} diff --git a/tests/fe/fe_support_points_sys_09/cmp/generic b/tests/fe/fe_support_points_sys_09/cmp/generic new file mode 100644 index 0000000000..8518f8d62f --- /dev/null +++ b/tests/fe/fe_support_points_sys_09/cmp/generic @@ -0,0 +1,6 @@ + +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(1)2<2> cell support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(1)2<2> face 0 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(1)2<2> face 1 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(1)2<2> face 2 support points +DEAL::FE_DGQ<2>(3)1FE_Nedelec<2>(1)2<2> face 3 support points diff --git a/tests/fe/fe_support_points_sys_10.cc b/tests/fe/fe_support_points_sys_10.cc new file mode 100644 index 0000000000..4ff239f555 --- /dev/null +++ b/tests/fe/fe_support_points_sys_10.cc @@ -0,0 +1,31 @@ +//---------------------------------------------------------------------- +// $Id$ +// Version: $Name$ +// +// Copyright (C) 2007, 2008, 2010, 2012 by the deal.II authors +// +// This file is subject to QPL and may not be distributed +// without copyright and license information. Please refer +// to the file deal.II/doc/license.html for the text and +// further information on this license. +// +//---------------------------------------------------------------------- +// +// Compute support points + +#include "../tests.h" +#include "fe_support_points_common.h" + + +int +main() +{ + initlog(__FILE__); + deallog.threshold_double(1.e-10); + + CHECK_SYS3(FE_Nedelec<2>(1), 1, + FESystem<2>(FE_DGQ<2>(3),3), 1, + FESystem<2>(FE_Q<2>(2),3, + FE_Nedelec<2>(1),2),2, + 2); +} diff --git a/tests/fe/fe_support_points_sys_10/cmp/generic b/tests/fe/fe_support_points_sys_10/cmp/generic new file mode 100644 index 0000000000..05fbb0adec --- /dev/null +++ b/tests/fe/fe_support_points_sys_10/cmp/generic @@ -0,0 +1,6 @@ + +DEAL::FE_Nedelec<2>(1)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(1),2)2<2> cell support points +DEAL::FE_Nedelec<2>(1)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(1),2)2<2> face 0 support points +DEAL::FE_Nedelec<2>(1)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(1),2)2<2> face 1 support points +DEAL::FE_Nedelec<2>(1)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(1),2)2<2> face 2 support points +DEAL::FE_Nedelec<2>(1)1FESystem<2>(FE_DGQ<2>(3),3)1FESystem<2>(FE_Q<2>(2),3, FE_Nedelec<2>(1),2)2<2> face 3 support points -- 2.39.5