From: bangerth Date: Tue, 19 Sep 2006 02:59:35 +0000 (+0000) Subject: New test that presently fails. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=33cf9fc3cf3301889bef860cbcc81c2bcf54ba3f;p=dealii-svn.git New test that presently fails. git-svn-id: https://svn.dealii.org/trunk@13927 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/hp/crash_09.cc b/tests/hp/crash_09.cc new file mode 100644 index 0000000000..0fc2f56c96 --- /dev/null +++ b/tests/hp/crash_09.cc @@ -0,0 +1,96 @@ +//---------------------------- crash_09.cc --------------------------- +// $Id: crash_09.cc 12732 2006-03-28 23:15:45Z wolf $ +// Version: $Name$ +// +// Copyright (C) 2006 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. +// +//---------------------------- crash_09.cc --------------------------- + + +// a test where a degree of freedom was constrained multiple times, +// but with different weights. see the hp paper for more on this + +char logname[] = "crash_09/output"; + + +#include "../tests.h" +#include +#include +#include +#include + +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include + + + +int main () +{ + std::ofstream logfile(logname); + logfile.precision (3); + + deallog.attach(logfile); + deallog.depth_console(0); + deallog.threshold_double(1.e-10); + + + // create a mesh like this: + // + // *---*---*---* + // | 6 | 7 | 8 | + // *---*---*---* + // | 3 | 4 | 5 | + // *---*---*---* + // | 0 | 1 | 2 | + // *---*---*---* + Triangulation<3> triangulation; + std::vector subdivisions (3, 2); + subdivisions[2] = 1; + GridGenerator::subdivided_hyper_rectangle (triangulation, subdivisions, + Point<3>(), Point<3>(2,2,1)); + + hp::FECollection<3> fe; + fe.push_back (FE_Q<3>(1)); + fe.push_back (FE_Q<3>(2)); + fe.push_back (FE_Q<3>(3)); + + hp::DoFHandler<3> dof_handler(triangulation); + + hp::DoFHandler<3>::active_cell_iterator + cell = dof_handler.begin_active(); + cell->set_active_fe_index (0); + ++cell; + cell->set_active_fe_index (1); + ++cell; + cell->set_active_fe_index (0); + ++cell; + cell->set_active_fe_index (2); + + dof_handler.distribute_dofs (fe); + + ConstraintMatrix constraints; + DoFTools::make_hanging_node_constraints (dof_handler, + constraints); + constraints.close (); + + constraints.print (deallog.get_file_stream()); +} + diff --git a/tests/hp/crash_09/cmp/generic b/tests/hp/crash_09/cmp/generic new file mode 100644 index 0000000000..656b9cc7c5 --- /dev/null +++ b/tests/hp/crash_09/cmp/generic @@ -0,0 +1,113 @@ + + 8 3: 0.500 + 8 14: 0.500 + 9 6: 0.500 + 9 17: 0.500 + 10 3: 0.500 + 10 6: 0.500 + 11 14: 0.500 + 11 17: 0.500 + 22 1: 0.500 + 22 3: 0.500 + 24 4: 0.500 + 24 6: 0.500 + 25 3: 0.500 + 25 6: 0.500 + 26 3: 0.167 + 26 6: 0.167 + 26 23: 0.667 + 27 3: 0.333 + 27 6: 0.333 + 27 23: 0.333 + 28 4: 0.333 + 28 6: 0.667 + 29 4: 0.167 + 29 6: 0.833 + 30 4: 0.167 + 30 6: 0.167 + 30 23: 0.667 + 31 4: 0.333 + 31 6: 0.333 + 31 23: 0.333 + 32 3: 0.333 + 32 6: 0.667 + 33 3: 0.167 + 33 6: 0.833 + 38 2: 0.500 + 38 3: 0.500 + 41 2: 0.333 + 41 3: 0.333 + 41 40: 0.333 + 42 2: 0.167 + 42 3: 0.167 + 42 40: 0.667 + 43 3: 0.833 + 43 14: 0.167 + 44 3: 0.667 + 44 14: 0.333 + 45 2: 0.333 + 45 3: 0.667 + 46 2: 0.167 + 46 3: 0.833 + 47 3: 0.167 + 47 14: 0.167 + 47 40: 0.667 + 48 3: 0.333 + 48 14: 0.333 + 48 40: 0.333 + 53 13: 0.500 + 53 14: 0.500 + 54 3: 0.500 + 54 14: 0.500 + 55 6: 0.500 + 55 7: 0.500 + 56 6: 0.500 + 56 17: 0.500 + 59 6: 0.333 + 59 17: 0.667 + 60 6: 0.167 + 60 17: 0.833 + 61 17: 0.167 + 61 18: 0.167 + 61 57: 0.667 + 62 17: 0.333 + 62 18: 0.333 + 62 57: 0.333 + 63 6: 0.333 + 63 17: 0.333 + 63 57: 0.333 + 64 6: 0.167 + 64 17: 0.167 + 64 57: 0.667 + 65 17: 0.833 + 65 18: 0.167 + 66 17: 0.667 + 66 18: 0.333 + 71 17: 0.500 + 71 18: 0.500 + 72 14: 0.833 + 72 16: 0.167 + 73 14: 0.667 + 73 16: 0.333 + 74 14: 0.333 + 74 17: 0.333 + 74 85: 0.333 + 75 14: 0.167 + 75 17: 0.167 + 75 85: 0.667 + 76 14: 0.833 + 76 17: 0.167 + 77 14: 0.667 + 77 17: 0.333 + 78 14: 0.333 + 78 16: 0.333 + 78 85: 0.333 + 79 14: 0.167 + 79 16: 0.167 + 79 85: 0.667 + 84 14: 0.500 + 84 17: 0.500 + 86 17: 0.500 + 86 19: 0.500 + 87 14: 0.500 + 87 16: 0.500