From: David Wells <drwells@email.unc.edu>
Date: Wed, 14 Apr 2021 04:03:20 +0000 (-0400)
Subject: Get rid of create_triangulation_compatibility in the DB Mesh reader.
X-Git-Tag: v9.3.0-rc1~221^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d2f1977aa675e8a0675e718a3106fdc1040ee688;p=dealii.git

Get rid of create_triangulation_compatibility in the DB Mesh reader.
---

diff --git a/source/grid/grid_in.cc b/source/grid/grid_in.cc
index 5f82626d63..8ac1ef56b1 100644
--- a/source/grid/grid_in.cc
+++ b/source/grid/grid_in.cc
@@ -1319,7 +1319,7 @@ GridIn<dim, spacedim>::read_dbmesh(std::istream &in)
       cells.emplace_back();
       for (const unsigned int i : GeometryInfo<dim>::vertex_indices())
         {
-          in >> cells.back().vertices[i];
+          in >> cells.back().vertices[GeometryInfo<dim>::ucd_to_deal[i]];
 
           AssertThrow((cells.back().vertices[i] >= 1) &&
                         (static_cast<unsigned int>(cells.back().vertices[i]) <=
@@ -1357,9 +1357,10 @@ GridIn<dim, spacedim>::read_dbmesh(std::istream &in)
   GridTools::delete_unused_vertices(vertices, cells, subcelldata);
   // ...and cells
   GridReordering<dim, spacedim>::invert_all_cells_of_negative_grid(vertices,
-                                                                   cells);
-  GridReordering<dim, spacedim>::reorder_cells(cells);
-  tria->create_triangulation_compatibility(vertices, cells, subcelldata);
+                                                                   cells,
+                                                                   true);
+  GridReordering<dim, spacedim>::reorder_cells(cells, true);
+  tria->create_triangulation(vertices, cells, subcelldata);
 }