From fea3aa3cc4cfba751496946f647c81368515e033 Mon Sep 17 00:00:00 2001 From: heltai Date: Tue, 17 May 2011 13:34:50 +0000 Subject: [PATCH] Added failing test for transform_real_to_unit in codim one. git-svn-id: https://svn.dealii.org/trunk@23711 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/codim_one/mapping_03.cc | 67 ++++++++++++++++++++++++++ tests/codim_one/mapping_03/cmp/generic | 0 2 files changed, 67 insertions(+) create mode 100644 tests/codim_one/mapping_03.cc create mode 100644 tests/codim_one/mapping_03/cmp/generic diff --git a/tests/codim_one/mapping_03.cc b/tests/codim_one/mapping_03.cc new file mode 100644 index 0000000000..a4090f45e9 --- /dev/null +++ b/tests/codim_one/mapping_03.cc @@ -0,0 +1,67 @@ +//---------------------------- mapping_02.cc --------------------------- +// $Id: bem.cc 22693 2010-11-11 20:11:47Z kanschat $ +// Version: $Name$ +// +// Copyright (C) 2010, 2011 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. +// +//---------------------------- mapping_02.cc --------------------------- + + +// like _01, but use a quadratic mapping. since we now map line segments to +// curves, the normal vectors at different quadrature points should no longer +// be parallel + +#include "../tests.h" + +#include +#include +#include +#include +#include +#include +#include +#include +#include + + +template +void test () +{ + Triangulation mesh; + GridGenerator::hyper_cube(mesh); + + QGauss quadrature(dim == 2 ? 3 : 2); + MappingQ mapping(1); + Point p; + + // Try to project a point on the + // surface + for(unsigned int i=0; i q = + mapping.transform_real_to_unit_cell(mesh.begin_active(), p); + + deallog << "P: " << p + << ", on unit: " << q << endl; + +} + + + +int main () +{ + std::ofstream logfile("mapping_03/output"); + deallog.attach(logfile); + deallog.depth_console(0); + + test<2> (); + test<3> (); + + return 0; +} diff --git a/tests/codim_one/mapping_03/cmp/generic b/tests/codim_one/mapping_03/cmp/generic new file mode 100644 index 0000000000..e69de29bb2 -- 2.39.5