From: Martin Kronbichler Date: Fri, 21 Jan 2022 16:24:14 +0000 (+0100) Subject: New test case X-Git-Tag: v9.4.0-rc1~574^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13278%2Fhead;p=dealii.git New test case --- diff --git a/tests/mappings/mapping_q_points_real_to_unit_02.cc b/tests/mappings/mapping_q_points_real_to_unit_02.cc new file mode 100644 index 0000000000..3438764f7f --- /dev/null +++ b/tests/mappings/mapping_q_points_real_to_unit_02.cc @@ -0,0 +1,75 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2020 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.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + + +// check MappingQ::transform_real_to_unit_points on MappingQ1 with a point +// that would at some point return an invalid answer due to a wrongly computed +// affine approximation, using a point extracted from +// tests/particles/particle_handler_21. + +#include + +#include +#include + +#include "../tests.h" + + +template +void +test() +{ + Triangulation tria; + GridGenerator::hyper_ball(tria, Point(), 1.0); + tria.refine_global(3); + + // pick out a particular cell + const typename Triangulation::cell_iterator coarse_cell(&tria, 0, 5); + const auto cell = coarse_cell->child(0)->child(5)->child(7); + + std::vector> real_points{ + Point<3>(-0.0489699, -0.859053, -0.0489699)}; + std::vector> unit_points(real_points.size()); + + for (unsigned int degree = 1; degree < 4; ++degree) + { + MappingQ mapping(degree); + mapping.transform_points_real_to_unit_cell(cell, + real_points, + unit_points); + + deallog << "Transform on cell with center: " << cell->center(true) + << " with mapping degree " << degree << std::endl; + deallog << "Combined transform " << real_points[0] << " gives " + << unit_points[0] << std::endl; + + deallog << "Classic transform " << real_points[0] << " gives " + << mapping.transform_real_to_unit_cell(cell, real_points[0]) + << std::endl; + } + + deallog << "OK" << std::endl; +} + + + +int +main() +{ + initlog(); + deallog << std::setprecision(10); + + test<3>(); +} diff --git a/tests/mappings/mapping_q_points_real_to_unit_02.output b/tests/mappings/mapping_q_points_real_to_unit_02.output new file mode 100644 index 0000000000..d61adada59 --- /dev/null +++ b/tests/mappings/mapping_q_points_real_to_unit_02.output @@ -0,0 +1,11 @@ + +DEAL::Transform on cell with center: -0.08086235035 -0.8200448543 -0.08086235035 with mapping degree 1 +DEAL::Combined transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::Classic transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::Transform on cell with center: -0.08086235035 -0.8200448543 -0.08086235035 with mapping degree 2 +DEAL::Combined transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::Classic transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::Transform on cell with center: -0.08086235035 -0.8200448543 -0.08086235035 with mapping degree 3 +DEAL::Combined transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::Classic transform -0.04896990000 -0.8590530000 -0.04896990000 gives 0.7082205181 0.2510034894 0.7082205181 +DEAL::OK