Bump version, update readme/changelog, add tests

Author: FranzBangar

CHANGELOG.md

@@ -4,6 +4,10 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+# [1.7.1]
+### Changed
+- Bugfix: wrong grading propagation on inverted wires
+
 # [1.7.0]
 ### Added
 - Automatic grading:

examples/bug.py

@@ -0,0 +1,23 @@
+import numpy as np
+
+import classy_blocks as cb
+
+box = cb.Box([0, 0, 0], [1, 1, 1])
+        
+base_face = box.bottom_face
+neighbour_face = base_face.copy().translate([4, 0, 0])
+common_face = base_face.copy().rotate(np.pi/2, [0, 1, 0]).translate([2, 0, 2])
+
+left_loft = cb.Loft(base_face, common_face)
+right_loft = cb.Loft(neighbour_face, common_face.copy().shift(2).invert())
+
+for axis in range(3):
+    left_loft.chop(axis, start_size=0.05, total_expansion=5)
+
+right_loft.chop(2, count=10)
+
+mesh = cb.Mesh()
+mesh.add(left_loft)
+mesh.add(right_loft)
+
+mesh.write("case/system/blockMeshDict", debug_path="debug.vtk")

pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "classy_blocks"
-version = "1.7.0"
+version = "1.7.1"
 description = "Python classes for easier creation of openFoam's blockMesh dictionaries."
 readme = "README.md"
 license = { file = "LICENSE" }

tests/test_bugs/__init__.py

@@ -0,0 +1,43 @@
+import unittest
+from typing import get_args
+
+import numpy as np
+
+import classy_blocks as cb
+from classy_blocks.cbtyping import DirectionType
+
+
+class GradingBugTests(unittest.TestCase):
+    def test_invert_grading(self):
+        # Bug case; two blocks with separate bottom faces share the same top face.
+        # Grading in one direction must be inverted
+        #          /|\
+        #        / / \ \
+        #      /  /   \  \
+        #    /   /     \   \
+        #  /____/   ^   \____\
+        #  base  common  neighbour
+
+        box = cb.Box([0, 0, 0], [1, 1, 1])
+
+        base_face = box.bottom_face
+        neighbour_face = base_face.copy().translate([4, 0, 0])
+        common_face = base_face.copy().rotate(np.pi / 2, [0, 1, 0]).translate([2, 0, 2])
+
+        left_loft = cb.Loft(base_face, common_face)
+        right_loft = cb.Loft(neighbour_face, common_face.copy().shift(2).invert())
+
+        for axis in get_args(DirectionType):
+            left_loft.chop(axis, start_size=0.05, total_expansion=5)
+
+        right_loft.chop(2, count=10)
+
+        mesh = cb.Mesh()
+        mesh.add(left_loft)
+        mesh.add(right_loft)
+
+        mesh.assemble()
+        mesh.block_list.assemble()
+
+        self.assertEqual(mesh.blocks[0].format_grading(), "simpleGrading ( 5 5 5 )")
+        self.assertEqual(mesh.blocks[1].format_grading(), "simpleGrading ( 0.2 5 1 )")