From e25314fbef62d6514eb8d32745a82ae36822dd4a Mon Sep 17 00:00:00 2001 From: Daniel Arndt Date: Fri, 26 May 2023 08:50:45 -0400 Subject: [PATCH] Remove unused file --- tests/cuda/create_mesh.h | 144 --------------------------------------- 1 file changed, 144 deletions(-) delete mode 100644 tests/cuda/create_mesh.h diff --git a/tests/cuda/create_mesh.h b/tests/cuda/create_mesh.h deleted file mode 100644 index e503ee441f..0000000000 --- a/tests/cuda/create_mesh.h +++ /dev/null @@ -1,144 +0,0 @@ -// --------------------------------------------------------------------- -// -// Copyright (C) 2008 - 2021 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. -// -// --------------------------------------------------------------------- - - -// this creates a mesh that contains cells of all different kinds detected in -// the MatrixFree class: Use a mesh that consists of a square cell, then a -// triangular cell, a parallelepiped cell and two trapezoidal cells, where one -// is very close to being Cartesian (by 1e-8). - -#include - -#include - -#include -#include - - -void -create_mesh(Triangulation<2> &tria, const double scale_grid = 1.) -{ - const unsigned int dim = 2; - std::vector> points(12); - - // build the mesh layer by layer from points - - // 1. cube cell - points[0] = Point(0, 0); - points[1] = Point(0, 1); - points[2] = Point(1, 0); - points[3] = Point(1, 1); - - // 2. rectangular cell - points[4] = Point(3., 0); - points[5] = Point(3., 1); - - // 3. parallelogram cell - points[6] = Point(5., 1.); - points[7] = Point(5., 2.); - - // almost square cell (but trapezoidal by - // 1e-8) - points[8] = Point(6., 1.); - points[9] = Point(6., 2. + 1e-8); - - // apparently trapezoidal cell - points[10] = Point(7., 1.4); - points[11] = Point(7.5, numbers::PI); - - if (scale_grid != 1.) - for (unsigned int i = 0; i < points.size(); ++i) - points[i] *= scale_grid; - - - // connect the points to cells - std::vector> cells(5); - for (unsigned int i = 0; i < 5; ++i) - { - cells[i].vertices[0] = 0 + 2 * i; - cells[i].vertices[1] = 2 + 2 * i; - cells[i].vertices[2] = 1 + 2 * i; - cells[i].vertices[3] = 3 + 2 * i; - cells[i].material_id = 0; - } - - tria.create_triangulation(points, cells, SubCellData()); -} - - - -void -create_mesh(Triangulation<3> &tria, const double scale_grid = 1.) -{ - const unsigned int dim = 3; - std::vector> points(24); - - // build the mesh layer by layer from points - - // 1. cube cell - points[0] = Point(0, 0, 0); - points[1] = Point(0, 1., 0); - points[2] = Point(0, 0, 1); - points[3] = Point(0, 1., 1); - points[4] = Point(1., 0, 0); - points[5] = Point(1., 1., 0); - points[6] = Point(1., 0, 1); - points[7] = Point(1., 1., 1); - - // 2. rectangular cell - points[8] = Point(3., 0, 0); - points[9] = Point(3., 1, 0); - points[10] = Point(3., 0, 1); - points[11] = Point(3., 1, 1); - - // 3. parallelogram cell - points[12] = Point(5., 1., 1.); - points[13] = Point(5., 2., 1.); - points[14] = Point(5., 1., 2.); - points[15] = Point(5., 2., 2.); - - // almost square cell (but trapezoidal by - // 1e-8 in y-direction) - points[16] = Point(6., 1., 1.); - points[17] = Point(6., 2. + 1e-8, 1.); - points[18] = Point(6., 1., 2.); - points[19] = Point(6., 2., 2.); - - // apparently trapezoidal cell - points[20] = Point(7., 1.4, 1.2231); - points[21] = Point(7.5, numbers::PI, 1.334); - points[22] = Point(7., 1.5, 7.1); - points[23] = Point(7.5, 3.8, 2.99); - - if (scale_grid != 1.) - for (unsigned int i = 0; i < points.size(); ++i) - points[i] *= scale_grid; - - // connect the points to cells - std::vector> cells(5); - for (unsigned int i = 0; i < 5; ++i) - { - cells[i].vertices[0] = 0 + 4 * i; - cells[i].vertices[1] = 4 + 4 * i; - cells[i].vertices[2] = 1 + 4 * i; - cells[i].vertices[3] = 5 + 4 * i; - cells[i].vertices[4] = 2 + 4 * i; - cells[i].vertices[5] = 6 + 4 * i; - cells[i].vertices[6] = 3 + 4 * i; - cells[i].vertices[7] = 7 + 4 * i; - cells[i].material_id = 0; - } - tria.create_triangulation(points, cells, SubCellData()); -} -- 2.39.5