From: Martin Kronbichler Date: Sat, 27 Aug 2016 13:51:16 +0000 (+0200) Subject: Add similar test to matrix_vector_14b for the FE_Q case. X-Git-Tag: v8.5.0-rc1~729^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F3024%2Fhead;p=dealii.git Add similar test to matrix_vector_14b for the FE_Q case. --- diff --git a/tests/matrix_free/matrix_vector_04b.cc b/tests/matrix_free/matrix_vector_04b.cc new file mode 100644 index 0000000000..a6cc1f0a82 --- /dev/null +++ b/tests/matrix_free/matrix_vector_04b.cc @@ -0,0 +1,48 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + + + +// Same as matrix_vector_04, but using more quadrature points than the FE +// degree would suggest + +#include "../tests.h" + +std::ofstream logfile("output"); + +#include "matrix_vector_common.h" + + +template +void test () +{ + const SphericalManifold manifold; + Triangulation tria; + GridGenerator::hyper_shell (tria, Point(), + 0.5, 1., 96, true); + tria.set_all_manifold_ids(0); + tria.set_manifold (0, manifold); + if (dim == 2) + tria.refine_global (2); + + FE_Q fe (fe_degree); + DoFHandler dof (tria); + dof.distribute_dofs(fe); + ConstraintMatrix constraints; + DoFTools::make_hanging_node_constraints(dof, constraints); + constraints.close(); + + do_test (dof, constraints); +} diff --git a/tests/matrix_free/matrix_vector_04b.output b/tests/matrix_free/matrix_vector_04b.output new file mode 100644 index 0000000000..c8d1cd1e9f --- /dev/null +++ b/tests/matrix_free/matrix_vector_04b.output @@ -0,0 +1,13 @@ + +DEAL:2d::Testing FE_Q<2>(1) +DEAL:2d::Norm of difference: 0 +DEAL:2d:: +DEAL:2d::Testing FE_Q<2>(2) +DEAL:2d::Norm of difference: 0 +DEAL:2d:: +DEAL:3d::Testing FE_Q<3>(1) +DEAL:3d::Norm of difference: 0 +DEAL:3d:: +DEAL:3d::Testing FE_Q<3>(2) +DEAL:3d::Norm of difference: 0 +DEAL:3d::