increase size of coarsest level

src/Classic/AbsBeamline/FlexibleCollimator.cpp

@@ -232,6 +232,14 @@ void FlexibleCollimator::setDescription(const std::string &desc) {
         urc[1] = std::max(urc[1], bb.center_m[1] + 0.5 * bb.height_m);
     }
 
+    double width = urc[0] - llc[0];
+    double height = urc[1] - llc[1];
+
+    llc[0] -= 1e-3 * width;
+    urc[0] += 1e-3 * width;
+    llc[1] -= 1e-3 * height;
+    urc[1] += 1e-3 * height;
+
     bb_m = mslang::BoundingBox(llc, urc);
 
     tree_m.bb_m = bb_m;

src/Classic/Utilities/MSLang.h

@@ -193,6 +193,12 @@ namespace mslang {
         void writeGnuplot(std::ofstream &out) const {
             out << "# level: " << level_m << ", size: " << objects_m.size() << std::endl;
             bb_m.writeGnuplot(out);
+            out << "# num holes: " << objects_m.size() << std::endl;
+            for (const Base *obj: objects_m) {
+                obj->writeGnuplot(out);
+            }
+            out << std::endl;
+
             if (nodes_m != 0) {
                 for (unsigned int i = 0; i < 4u; ++ i) {
                     nodes_m[i].writeGnuplot(out);

tools/mslang/mslang.cpp

@@ -71,6 +71,13 @@ int main(int argc, char *argv[])
                 urc[1] = std::max(urc[1], bb.center_m[1] + 0.5 * bb.height_m);
             }
 
+            double width = urc[0] - llc[0];
+            double height = urc[1] - llc[1];
+            llc[0] -= 1e-3 * width;
+            urc[0] += 1e-3 * width;
+            llc[1] -= 1e-3 * height;
+            urc[1] += 1e-3 * width;
+
             mslang::QuadTree tree;
             tree.bb_m = mslang::BoundingBox(llc, urc);
             tree.objects_m.insert(tree.objects_m.end(), baseBlocks.begin(), baseBlocks.end());