From: Peter Munch Date: Fri, 17 Jan 2020 04:43:22 +0000 (+0100) Subject: Tests to ensure that the order of pft and be the base triangulations are the same X-Git-Tag: v9.2.0-rc1~607^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F9362%2Fhead;p=dealii.git Tests to ensure that the order of pft and be the base triangulations are the same --- diff --git a/tests/fullydistributed_grids/copy_distributed_tria_03.cc b/tests/fullydistributed_grids/copy_distributed_tria_03.cc new file mode 100644 index 0000000000..5faddbff15 --- /dev/null +++ b/tests/fullydistributed_grids/copy_distributed_tria_03.cc @@ -0,0 +1,94 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2019 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. +// +// --------------------------------------------------------------------- + + +// Create a distributed triangulation with multigrid levels, copy it, +// and show that the CellIds have not changed. + +#include + +#include +#include +#include +#include + +#include + +#include +#include +#include +#include + +#include "./tests.h" + +using namespace dealii; + +template +void +test(int n_refinements, const int n_subdivisions, MPI_Comm comm) +{ + // create pdt + parallel::distributed::Triangulation tria_pdt( + comm, + dealii::Triangulation::none, + parallel::distributed::Triangulation::construct_multigrid_hierarchy); + GridGenerator::subdivided_hyper_cube(tria_pdt, n_subdivisions); + tria_pdt.refine_global(n_refinements); + + // create instance of pft + parallel::fullydistributed::Triangulation tria_pft(comm); + + // extract relevant information form pdt + auto construction_data = parallel::fullydistributed::Utilities:: + create_construction_data_from_triangulation(tria_pdt, comm, true); + + // actually create triangulation + tria_pft.create_triangulation(construction_data); + + // copare centers of cells + for (auto &cell : tria_pft.active_cell_iterators()) + { + CellId id = cell->id(); + auto cell_base = id.to_cell(tria_pdt); + for (unsigned int d = 0; d < dim; d++) + Assert(std::abs(cell->center()[d] - cell_base->center()[d]) < 1e-9, + ExcMessage("Cells do not match")); + } +} + +int +main(int argc, char *argv[]) +{ + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + mpi_initlog(); + + const int dim = 2; + const MPI_Comm comm = MPI_COMM_WORLD; + + { + deallog.push("2d"); + const int n_refinements = 3; + const int n_subdivisions = 5; + test<2>(n_refinements, n_subdivisions, comm); + deallog.pop(); + } + { + deallog.push("3d"); + const int n_refinements = 3; + const int n_subdivisions = 4; + test<3>(n_refinements, n_subdivisions, comm); + deallog.pop(); + } +} diff --git a/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=1.output b/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=1.output new file mode 100644 index 0000000000..8b13789179 --- /dev/null +++ b/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=1.output @@ -0,0 +1 @@ + diff --git a/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=5.output b/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=5.output new file mode 100644 index 0000000000..8b13789179 --- /dev/null +++ b/tests/fullydistributed_grids/copy_distributed_tria_03.with_p4est=true.mpirun=5.output @@ -0,0 +1 @@ + diff --git a/tests/fullydistributed_grids/copy_serial_tria_08.cc b/tests/fullydistributed_grids/copy_serial_tria_08.cc new file mode 100644 index 0000000000..ffa76144a3 --- /dev/null +++ b/tests/fullydistributed_grids/copy_serial_tria_08.cc @@ -0,0 +1,94 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2019 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. +// +// --------------------------------------------------------------------- + + +// Create a serial triangulation with multigrid level, copy it, and +// show that the CellIds have not changed. + +#include + +#include +#include +#include +#include + +#include + +#include +#include +#include +#include + +#include "./tests.h" + +using namespace dealii; + +template +void +test(int n_refinements, MPI_Comm comm) +{ + // create serial triangulation + Triangulation basetria( + Triangulation::limit_level_difference_at_vertices); + GridGenerator::hyper_L(basetria); + basetria.refine_global(n_refinements); + + GridTools::partition_triangulation_zorder( + Utilities::MPI::n_mpi_processes(comm), basetria); + GridTools::partition_multigrid_levels(basetria); + + // create instance of pft + parallel::fullydistributed::Triangulation tria_pft(comm); + + // extract relevant information form serial triangulation + auto construction_data = parallel::fullydistributed::Utilities:: + create_construction_data_from_triangulation(basetria, comm, true); + + // actually create triangulation + tria_pft.create_triangulation(construction_data); + + + for (auto &cell : tria_pft.active_cell_iterators()) + { + CellId id = cell->id(); + auto cell_base = id.to_cell(basetria); + // Assert(cell->center() == cell_base->center(), + // ExcMessage("Cells do not match")); + for (unsigned int d = 0; d < dim; d++) + Assert(std::abs(cell->center()[d] - cell_base->center()[d]) < 1e-9, + ExcMessage("Cells do not match")); + } +} + +int +main(int argc, char *argv[]) +{ + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + mpi_initlog(); + + const int n_refinements = 3; + const MPI_Comm comm = MPI_COMM_WORLD; + + { + deallog.push("2d"); + test<2>(n_refinements, comm); + deallog.pop(); + } + { + deallog.push("3d"); + test<3>(n_refinements, comm); + deallog.pop(); + } +} diff --git a/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=1.output b/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=1.output new file mode 100644 index 0000000000..8b13789179 --- /dev/null +++ b/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=1.output @@ -0,0 +1 @@ + diff --git a/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=5.output b/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=5.output new file mode 100644 index 0000000000..8b13789179 --- /dev/null +++ b/tests/fullydistributed_grids/copy_serial_tria_08.mpirun=5.output @@ -0,0 +1 @@ +