From 17ea5cf57aa48cd09b65e62c70f204eb69e381ff Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Mon, 22 Jul 2024 18:10:15 +0200 Subject: [PATCH] New test case --- tests/distributed_grids/2d_periodicity_01.cc | 73 ++++++++++++++++++++ 1 file changed, 73 insertions(+) create mode 100644 tests/distributed_grids/2d_periodicity_01.cc diff --git a/tests/distributed_grids/2d_periodicity_01.cc b/tests/distributed_grids/2d_periodicity_01.cc new file mode 100644 index 0000000000..0e570c75ca --- /dev/null +++ b/tests/distributed_grids/2d_periodicity_01.cc @@ -0,0 +1,73 @@ +// ------------------------------------------------------------------------ +// +// SPDX-License-Identifier: LGPL-2.1-or-later +// Copyright (C) 2024 by the deal.II authors +// +// This file is part of the deal.II library. +// +// Part of the source code is dual licensed under Apache-2.0 WITH +// LLVM-exception OR LGPL-2.1-or-later. Detailed license information +// governing the source code and code contributions can be found in +// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II. +// +// ------------------------------------------------------------------------ + + + +// Test case with p4est with a generated mesh where the orientation on the +// lines does not match and we would previously run into a bug + + +#include + +#include +#include + +#include "../tests.h" + +int +main(int argc, char *argv[]) +{ + initlog(); + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + + constexpr unsigned int dim = 2; + parallel::distributed::Triangulation tria(MPI_COMM_WORLD); + + GridIn grid_in; + grid_in.attach_triangulation(tria); + std::ifstream in_stream(SOURCE_DIR "/2d_periodicity_01.inp"); + grid_in.read(in_stream, GridIn::abaqus); + + deallog << "Number of mesh elements: " << tria.n_active_cells() << std::endl; + + const BoundingBox bounding_box = GridTools::compute_bounding_box(tria); + for (const auto &cell : tria.cell_iterators()) + for (const auto &face : cell->face_iterators()) + if (face->at_boundary()) + for (unsigned int d = 0; d < dim; ++d) + { + if (std::abs(face->center()[d] - + bounding_box.get_boundary_points().first[d]) < 1e-8) + face->set_boundary_id(1 + 2 * d); + else if (std::abs(face->center()[d] - + bounding_box.get_boundary_points().second[d]) < + 1e-8) + face->set_boundary_id(2 + 2 * d); + } + + std::vector< + dealii::GridTools::PeriodicFacePair::cell_iterator>> + periodic_face_pairs; + + for (unsigned int d = 0; d < dim; ++d) + GridTools::collect_periodic_faces( + tria, 1 + 2 * d, 2 + 2 * d, d, periodic_face_pairs); + + deallog << "Number of periodic face pairs found: " + << periodic_face_pairs.size() << std::endl; + tria.add_periodicity(periodic_face_pairs); + + tria.refine_global(); + deallog << "OK" << std::endl; +} -- 2.39.5