fix for node ordering in tensor grid constructor such that no negativ…

…e cell volumes occur
WIAS-PDELib · Mar 22, 2022 · f9d47ae · f9d47ae · j-fu · Mar 22, 2022
1 parent 78761c2
commit f9d47ae
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "ExtendableGrids"
 uuid = "cfc395e8-590f-11e8-1f13-43a2532b2fa8"
 authors = ["Juergen Fuhrmann <[email protected]>, Christian Merdon  <[email protected]>"]
-version = "0.9.1"
+version = "0.9.2"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"

diff --git a/src/simplexgrid.jl b/src/simplexgrid.jl
@@ -591,7 +591,22 @@ function simplexgrid(grid2::ExtendableGrid, coordZ; bot_offset=0,cell_offset=0,t
         end
         ishift+=nnodes2
     end
-    simplexgrid(coord3,cells3,cellregions3,bfaces3,bfaceregions3)
+
+    xgrid = simplexgrid(coord3,cells3,cellregions3,bfaces3,bfaceregions3)
+
+    ## quick and dirty fix of local orderings to avoid negative CellVolumes
+    ## and problems with CellFinder
+    cellvolumes = xgrid[CellVolumes]
+    cellnodes = xgrid[CellNodes]
+    for cell = 1 : num_cells(xgrid)
+        if cellvolumes[cell] < 0
+            ishift = cellnodes[1,cell]
+            cellnodes[1,cell] = cellnodes[2,cell]
+            cellnodes[2,cell] = ishift
+            cellvolumes[cell] *= -1
+        end
+    end
+    return xgrid
 end