diff --git a/src/db/unit_tests/dbDeepRegionTests.cc b/src/db/unit_tests/dbDeepRegionTests.cc
index e8e012b42..cbe9a558c 100644
--- a/src/db/unit_tests/dbDeepRegionTests.cc
+++ b/src/db/unit_tests/dbDeepRegionTests.cc
@@ -1490,59 +1490,112 @@ TEST(26_BreakoutCells)
   db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au26.gds");
 }
 
-TEST(27_snap)
+TEST(27a_snap)
 {
+  db::Layout ly;
   {
-    db::Layout ly;
-    {
-      std::string fn (tl::testsrc ());
-      fn += "/testdata/algo/scale_and_snap.gds";
-      tl::InputStream stream (fn);
-      db::Reader reader (stream);
-      reader.read (ly);
-    }
-
-    db::cell_index_type top_cell_index = *ly.begin_top_down ();
-    db::Cell &top_cell = ly.cell (top_cell_index);
-
-    db::DeepShapeStore dss;
-
-    unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
-    db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
-    r1.set_merged_semantics (false);
-    db::Region r2 = r1.snapped (19, 19);
-
-    r2.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
-
-    CHECKPOINT();
-    db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au27.gds");
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
   }
 
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  db::DeepShapeStore dss;
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
+  r1.set_merged_semantics (false);
+  db::Region r2 = r1.snapped (19, 19);
+
+  r2.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au27.gds");
+}
+
+TEST(27b_snap)
+{
+  db::Layout ly;
   {
-    db::Layout ly;
-    {
-      std::string fn (tl::testsrc ());
-      fn += "/testdata/algo/scale_and_snap.gds";
-      tl::InputStream stream (fn);
-      db::Reader reader (stream);
-      reader.read (ly);
-    }
-
-    db::cell_index_type top_cell_index = *ly.begin_top_down ();
-    db::Cell &top_cell = ly.cell (top_cell_index);
-
-    db::DeepShapeStore dss;
-
-    unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
-    db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
-    r1.set_merged_semantics (false);
-    r1.snap (19, 19);
-
-    r1.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
-
-    CHECKPOINT();
-    db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au27.gds");
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
   }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  db::DeepShapeStore dss;
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
+  r1.set_merged_semantics (false);
+  r1.snap (19, 19);
+
+  r1.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au27.gds");
+}
+
+TEST(28a_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  db::DeepShapeStore dss;
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
+  r1.set_merged_semantics (false);
+  db::Region r2 = r1.scaled_and_snapped (19, 2, 10, 19, 2, 10);
+
+  r2.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au28.gds");
+}
+
+TEST(28b_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  db::DeepShapeStore dss;
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1), dss);
+  r1.set_merged_semantics (false);
+  r1.scale_and_snap (19, 2, 10, 19, 2, 10);
+
+  r1.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/deep_region_au28.gds");
 }
 
 TEST(100_Integration)
diff --git a/src/db/unit_tests/dbRegion.cc b/src/db/unit_tests/dbRegion.cc
index 9f1a4b40f..21ba05b60 100644
--- a/src/db/unit_tests/dbRegion.cc
+++ b/src/db/unit_tests/dbRegion.cc
@@ -28,6 +28,10 @@
 #include "dbRegionProcessors.h"
 #include "dbEdgesUtils.h"
 #include "dbBoxScanner.h"
+#include "dbReader.h"
+#include "dbTestSupport.h"
+
+#include "tlStream.h"
 
 #include <cstdio>
 
@@ -1391,6 +1395,106 @@ TEST(31)
   EXPECT_EQ (db::Region (db::Box (db::Point (0, 3999), db::Point (1001, 6000))).pull_overlapping (r).to_string (), "(1000,0;1000,4000;6000,4000;6000,0)");
 }
 
+TEST(32a_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1));
+  r1.set_merged_semantics (false);
+  db::Region r2 = r1.snapped (19, 19);
+
+  r2.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/region_au32.gds");
+}
+
+TEST(32b_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1));
+  r1.set_merged_semantics (false);
+  r1.snap (19, 19);
+
+  r1.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/region_au32.gds");
+}
+
+TEST(33a_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1));
+  r1.set_merged_semantics (false);
+  db::Region r2 = r1.scaled_and_snapped (19, 2, 10, 19, 2, 10);
+
+  r2.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/region_au33.gds");
+}
+
+TEST(33b_snap)
+{
+  db::Layout ly;
+  {
+    std::string fn (tl::testsrc ());
+    fn += "/testdata/algo/scale_and_snap.gds";
+    tl::InputStream stream (fn);
+    db::Reader reader (stream);
+    reader.read (ly);
+  }
+
+  db::cell_index_type top_cell_index = *ly.begin_top_down ();
+  db::Cell &top_cell = ly.cell (top_cell_index);
+
+  unsigned int l1 = ly.get_layer (db::LayerProperties (1, 0));
+  db::Region r1 (db::RecursiveShapeIterator (ly, top_cell, l1));
+  r1.set_merged_semantics (false);
+  r1.scale_and_snap (19, 2, 10, 19, 2, 10);
+
+  r1.insert_into (&ly, top_cell_index, ly.get_layer (db::LayerProperties (100, 0)));
+
+  CHECKPOINT();
+  db::compare_layouts (_this, ly, tl::testsrc () + "/testdata/algo/region_au33.gds");
+}
+
 TEST(100_Processors)
 {
   db::Region r;
diff --git a/testdata/algo/deep_region_au28.gds b/testdata/algo/deep_region_au28.gds
new file mode 100644
index 000000000..b7d0b6ef1
Binary files /dev/null and b/testdata/algo/deep_region_au28.gds differ
diff --git a/testdata/algo/region_au32.gds b/testdata/algo/region_au32.gds
new file mode 100644
index 000000000..f510c8dca
Binary files /dev/null and b/testdata/algo/region_au32.gds differ
diff --git a/testdata/algo/region_au33.gds b/testdata/algo/region_au33.gds
new file mode 100644
index 000000000..bae4cd78e
Binary files /dev/null and b/testdata/algo/region_au33.gds differ