From: Luca Heltai Date: Tue, 19 Sep 2017 13:22:45 +0000 (+0200) Subject: Added test with MappingQEulerian. X-Git-Tag: v9.0.0-rc1~1018^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c9dd6726723d1698dfbc1d206fcfe72d2b06d7d5;p=dealii.git Added test with MappingQEulerian. --- diff --git a/tests/grid/grid_tools_extract_used_vertices_03.cc b/tests/grid/grid_tools_extract_used_vertices_03.cc new file mode 100644 index 0000000000..f6d4ada182 --- /dev/null +++ b/tests/grid/grid_tools_extract_used_vertices_03.cc @@ -0,0 +1,72 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2001 - 2015 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 at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + +// Check extract used_vertices and find_closest_vertex when a mapping +// is present + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include + + +template +void test (const Point &p) +{ + deallog << "dim: " << dim << ", spacedim: " + << spacedim << std::endl; + + Triangulation tria; + GridGenerator::hyper_cube(tria); + tria.refine_global(1); + + // Vector fe + FESystem fe(FE_Q(1)^spacedim); + DoFHandler dh(tria); + dh.distribute_dofs(fe); + + // Displacement vector. + Vector disp(dh.n_dofs()); + for (unsigned int i=0; i, spacedim> mapping(1, dh, disp); + + auto m = GridTools::extract_used_vertices(tria, mapping); + + for (auto &e: m) + deallog << "Vertex: " << e.first << ": " << e.second << std::endl; + + auto i = GridTools::find_closest_vertex(m,p); + deallog << "Closest vertex to " << p + << ", v["<< i << "] :" + << m[i] << std::endl; +}; + + +int main () +{ + initlog(); + test<2,2> (Point<2>(.2,.2)); + test<2,2> (Point<2>(.6,.9)); + return 0; +} + diff --git a/tests/grid/grid_tools_extract_used_vertices_03.output b/tests/grid/grid_tools_extract_used_vertices_03.output new file mode 100644 index 0000000000..0331bd7a3a --- /dev/null +++ b/tests/grid/grid_tools_extract_used_vertices_03.output @@ -0,0 +1,23 @@ + +DEAL::dim: 2, spacedim: 2 +DEAL::Vertex: 0: 0.500000 0.500000 +DEAL::Vertex: 1: 1.50000 0.500000 +DEAL::Vertex: 2: 0.500000 1.50000 +DEAL::Vertex: 3: 1.50000 1.50000 +DEAL::Vertex: 4: 1.00000 0.500000 +DEAL::Vertex: 5: 0.500000 1.00000 +DEAL::Vertex: 6: 1.50000 1.00000 +DEAL::Vertex: 7: 1.00000 1.50000 +DEAL::Vertex: 8: 1.00000 1.00000 +DEAL::Closest vertex to 0.200000 0.200000, v[0] :0.500000 0.500000 +DEAL::dim: 2, spacedim: 2 +DEAL::Vertex: 0: 0.500000 0.500000 +DEAL::Vertex: 1: 1.50000 0.500000 +DEAL::Vertex: 2: 0.500000 1.50000 +DEAL::Vertex: 3: 1.50000 1.50000 +DEAL::Vertex: 4: 1.00000 0.500000 +DEAL::Vertex: 5: 0.500000 1.00000 +DEAL::Vertex: 6: 1.50000 1.00000 +DEAL::Vertex: 7: 1.00000 1.50000 +DEAL::Vertex: 8: 1.00000 1.00000 +DEAL::Closest vertex to 0.600000 0.900000, v[5] :0.500000 1.00000