From: bangerth Date: Fri, 14 Sep 2012 12:32:17 +0000 (+0000) Subject: Step one in breaking apart our monster testcase 'internals'. That testcase had an... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=3c2df429b3451fe9c350f9a2a9d72b0652e14b7d;p=dealii-svn.git 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 --- 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